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LES  PHENOMENES  DE  TRANSPORT  TURBULENT 

par  Parcel  BAR  RE  RE 

Rapport  du  " Working  group  09" 

Office  National  d’ Etudes  et  de  Recherches  Aerospatiales  (ON ERA) 
92320  Cha  til  Ion  (France) 


1  -  INTRODUCTION 


Ce  groupe  de  travail  du  PEP  a  ete  initiA  par  le  Professeur  A.  Ferri  et  nous  voudrions  pour  cadrer 
1' activity  de  ce  groupe  donner  les  commentaires  du  Professeur  A.  Ferri  sur  ce  problAme  de  l'Atude  des 
phAnomAnes  de  transport  en  milieu  turbulent  : 

"La  relation  entre  les  proprietes  de  transport  et  le  caractere  instationnaire  de  l'Acoulement  n'est 
pas  trAs  bien  comprise  et  les  progres  enregistres  dans  ce  domaine  ne  sont  pas  suffisants.  En  general  les 
approches  nouvelles  compliquent  1' analyse  sans  pour  autant  accroitre  la  faculte  de  comprehension.  De 
nombreuses  faqons  d'aborder  ce  probleme  ont  ete  suggerees  pour  predire  les  proprietes  de  l'ecoulement  du 
type  viscosite  tourbillonnaire  de  Prandtl,  A  savoir  : 

a)  l'emploi  en  particulier  de  moyennes  de  quantitAs  fluctuantes  et  de  proprietes  de  transport  en  relation 
avec  l'Anergie  cinetique  ;  quoiqu'il  en  soit  il  n'existe  pas  de  relation  rationnelle  pour  decrire  ces 
proprietes,  bien  que  de  nombreuses  relations  aient  ete  proposees  ;  de  plus  les  relations  ne  tiennent  pas 
compte  de  1' importance  des  larges  remous  (large  eddies)  en  comparaison  avec  les  remous  de  petites  dimen¬ 
sions,  ce  qui  est  trAs  important  et  A  la  base  du  modele  de  Prandtl  ; 

b)  l'utilisation  de  proprietes  qui  dependent  du  temps  (par  exemple  1' intermittence)  est  une  approche  plus 
physique  ;  cependant  les  methodes  pour  evaluer  les  proprietes  de  transport  A  partir  des  proprietes 
physiques  globales  de  l'ecoulement  ne  sont  pas  encore  disponibles.  II  est  bien  connu  que  de  nombreux 
paramAtres  affectent  les  proprietes  de  transport  :  l'Achelle  d'intensite  et  la  distribution  de  la 
turbulence  A  la  position  initiale,  la  generation  de  la  turbulence  due  aux  variations  du  cisaillement  et 
de  l'energie,  1 ' amortissement  de  la  turbulence  et  le  changement  de  frequence  et  d’echelle  de  celle-ci, 
la  generation  de  bruit  due  A  la  turbulence  qui  produit  les  fluctuations  de  pression  et  qui  diffuse  par 
un  mecanisme  d'ondes  et  qui  de  ce  fait  constitue  une  source  d' intermittence. 

Peu  d' informations  existent  actuellement  sur  le  plan  quantitatif  sur  l'effet  de  tels  paramAtres  et 
aucune  information  n'existe  sur  les  Acoulements  turbulents  tridimensionnels .  Le  "panel"  recommande  qu'un 
effort  important  soit  fait  pour  obtenir  une  meilleure  connaissance  de  telles  quantitAs  ou  de  tels  para¬ 
mAtres.  Le  manque  de  connaissances  en  physique  rend  les  resultats  obtenus  par  les  methodes  numeriques 
imprecis  et  souvent  leur  utilite  et  leur  capacite  sont  dangeureusement  surest imees.  Le  problAme  est  encore 
plus  complexe  pour  les  ecoulements  A  deux  phases". 

Ce  groupe  de  travail  etait  done  destine  A  accroitre  notre  comprehension  de  la  structure  de  l'Acoule- 
ment  pour  arriver  A  une  prediction  des  phAnomAnes  de  transport  dans  un  Acoulement  turbulent.  On  est  arrivA 
A  la  conclusion  que  pour  aboutir  A  ce  stade  de  comprehension  il  etait  trAs  important  de  possAder  de 
meilleurs  instruments  d'analyse  et  de  mesures  et  que,  en  particulier,  devaient  etre  placAs  en  prioritA  les 
mesures  de  grandeurs  scalaires  et  les  mesures  spatio-temporelles  des  fluctuations  de  ces  scalaires. 

On  dispose  par  ailleurs  de  techniques  plus  sures  pour  dAterminer  les  fluctuations  du  vecteur  vitesse 
principalement  par  des  techniques  lasers.  Ce  groupe  de  travail  a  selectionnA  quatre  grandeurs  scalaires  : 
concentration,  tempArature,  pression  et  masse  volumique  qui  paraissaient  signif icatives .  Le  souhait  du 
Professeur  A.  FERRI  avait  AtA  de  mettre  1' accent  sur  la  mesure  des  fluctuations  de  pression,  mais  comme 
il  s'avAre  que  cette  grandeur  est  difficile  A  mesurer  et  qu'il  n'existe  pas,  pour  l'instant,  de  mesure 
directe,  le  groupe  de  travail  a  prAfArA  Atendre  cette  suggestion  A  la  mesure  des  scalaires  en  gAnAral.  Ces 
mesures,  plus  gAnArales,  devaient  assurer  une  meilleure  comprAhension  des  Acoulements  turbulents  avec  ou 
sans  rAactions  chimiques  et  dans  un  milieu  oQ  la  masse  volumique  P  pouvait  fluctuer  (Acoulement  compres¬ 
sible  ou  incompressible). 

La  connaissance  des  phAnomAnes  physiques  et  chimiques  en  milieu  turbulent  a  des  applications  directes 
relatives  aux  methodes  de  calcul  quantitatives  des  ecoulements  dans  les  foyers  et  les  rAacteurs  chimiques, 
e'est  le  cas  en  particulier  des  chambres  de  combustion  de  foyers  principaux  et  de  rechauffe  de  turborAac- 
teurs,  des  foyers  de  statorAacteurs ,  des  chambres  de  lasers  chimiques,  des  rAactions  dans  le  sillage  des 
rAacteurs  et  des  fusAes,  dans  les  rAacteurs  homogAnes  ou  tubulaires  rencontres  en  gAnie  chimique,  dans 
devolution  de  la  composition  de  1'atmosphAre  et  de  la  stratosphAre . 

Avant  d'aborder  le  plan  de  ce  rapport  nous  voulons  souligner  deux  difficultAs  qui  ont  impose  l'orien- 
tation  choisie  : 

a)  -  Les  rAactions  chimiques  se  font  s'il  y  a  un  mAlange  A  toutes  petites  Achelles  des  substances  rAagis- 
santes,  la  turbulence  n'intervient  que  par  la  faqon  dont  elle  met  les  substances  en  prAsence,  elle 
peut  avoir  un  effet  bAnefique  si  elle  met  rapidement  en  prAsence  des  gaz  brfllAs  et  chauds  avec  des 
gaz  frais  prAmAlanges,  cet  effet  est  double  et  porte  A  la  fois  sur  le  paramAtre  tempArature  et  le 


parametre  concentration.  On  comprend  pourquoi,  dans  ce  cas,  la  turbulence  accroit  le  taux  de  produc¬ 
tion  chimique.  Mais  cet  effet  de  melange  par  turbulence  peut  dans  certains  cas  diminuer  le  taux  de 
production  :  c'est  la  situation  d'une  flamme  de  diffusion  pour  laquelle  le  comburant  et  le  combus¬ 
tible  sont  injectes  separement  dans  le  foyer,  la  structure  de  la  turbulence  peut  etre  telle  que 
certaines  molecules  comburantes  et  combustibles  ne  se  rencontrent  pas,  ce  qui  limite  le  taux  de 
production  chimique.  On  peut  done  concevoir  par  quel  mecanisme  la  turbulence  modifie  le  taux  des 
reactions  chimiques  mais  on  con^oit  moins  bien  comment  la  reaction  chimique,  par  l'energie  qu'elle 
apporte  localement,  par  1 ' accroissement  tres  rapide  du  volume  (diminution  de  P  )  et  de  la  tempera¬ 
ture  peut  modifier  la  turbulence,  c'est  pourquoi  des  mesures  dans  le  temps  et  l'espace  de  la  tempe¬ 
rature  et  des  concentrations  sont  necessaires. 

b)  -  II  n'en  reste  pas  moins  que  le  probleme  de  la  mesure  de  la  pression  et  de  ses  fluctuations  en  ecoule- 
ment  turbulent  demeure  encore  sans  solution.  Les  mesures  de  pression  moyenne  dans  un  tel  ecoulement 
sont  souvent  incertaines  ;  il  est  en  effet  necessaire  d'eviter,  d'eliminer  ou  d'etalonner  la  sensibi- 
lite  de  la  sonde  aux  fluctuations  de  pression  ainsi  que  la  sensibilite  parasite  de  la  direction  et  de 
l'intensite  des  fluctuations  de  pression  et  en  particulier  il  convient  de  determiner  dans  quelles 
conditions  1' ecoulement  autour  de  la  sonde  peut  etre  considere  comme  quasi-stationnaire .  En  general 
on  admet  comme  utilisable  l'etalonnage  effectue  en  ecoulement  stationnaire .  Des  verifications  dans 
plusieurs  conditions  experimentales  doivent  donner  la  precision  de  la  mesure  des  fluctuations  de 
pression  en  particulier  ces  experiences  de  verification  doivent  comporter  des  mesures  de  ' 
derriere  une  grille  engendrant  une  turbulence  isotropique,  les  mesures  de  peuvent  etre  deduites 
des  correlations  de  vitesse.  Cette  mesure  des  fluctuations  de  pression  comporte,  comme  nous  le 
verrons  par  la  suite,  de  nombreuses  difficultes. 

Apres  ces  remarques  preliminaires ,  1' expose  comprend  quatre  parties  : 

-  une  introduction  generale  concernant  la  structure  de  1' ecoulement  turbulent, 

-  la  mesure  des  temperatures  en  donnant  principalement  les  methodes  qui  nous  apparaissent  les  mieux  adap- 
tees  A  1' etude  de  la  turbulence, 

-  la  mesure  des  concentrations  en  insistant  plus  particulierement  sur  les  nouvelles  techniques  "laser", 

-  la  mesure  de  la  masse  volumique  et  les  procedes  de  visualisation, 

-  la  mesure  de  la  pression  pour  terminer  notre  rapport.  Nous  precisons  dans  la  conclusion  les  orientations 
que  l'on  peut  donner  sur  le  plan  scientifique  et  les  recommandations  &  proposer  au  PEP. 

2  -  STRUCTURE  VE  L’ ECOULEMENT  TURBULENT 


Pour  mieux  modeliser  les  ecoulements  turbulents  il  faut,  a  l'aide  des  mesures  vectorielles  et  scalai- 
res,  etre  capable  de  mieux  les  decrire  et  de  mieux  les  analyser.  La  premiere  difficulty  provient  de  ce 
qu'il  n'y  a  pas  un  ecoulement  turbulent  mais  une  grande  variete  d ' ecoulements  qui  peuvent  comporter  des 
structures  coherentes  et  des  structures  aleatoires.  Un  des  premiers  roles  de  1' analyse  est  de  distinguer 
ces  structures  coherentes  et  une  grande  partie  de  l'activite  mesure  repose  sur  cette  recherche. 

Pour  souligner  le  genre  de  solutions  que  peuvent  apporter  les  methodes  de  diagnostic,  nous  allons 
prendre  quelques  exemples  en  soulignant  tout  d'abord  les  difficultes  qui  decoulent  de  la  mesure  des  para- 
metres  vectoriels  ou  scalaires  : 

-  la  mesure  se  fait  en  un  point  fixe  et  il  est  souvent  difficile  de  suivre  des  traceurs  places  dans 
1' ecoulement, 

-  la  mesure,  si  elle  est  faite  par  un  procede  optique,  integre  sur  tout  le  trajet  du  faisceau  des  pheno- 
menes  et  si  l'on  introduit  ume  sonde  celle-ci  risque  de  perturber  1 ' ecoulement , 

-  il  faut  egalement  que  le  systeme  de  mesure  dans  son  ensemble  possede  un  temps  de  mesure  suffisamment 
court  pour  suivre  la  moindre  fluctuation, 

-  il  est  souvent  difficile  de  placer,  en  un  meme  point,  plusieurs  dispositifs  de  mesure  et  de  reconstituer 
avec  precision  la  structure  de  1 'ecoulement  aussi  bien  dans  l'espace  que  dans  le  temps  car  on  ne  dispose 
pas  de  suffisamment  d' informations . 

On  se  trouve  done  devant  un  champ  scalaire  et  vectoriel  complexe  resultant  de  la  superposition  et  de 
1' interaction  de  champs  elementaires . 

Considerons  par  exemple  la  structure  d'une  flamme  de  diffusion,  et  ce  sera  le  premier  exemple,  cons- 
titue  par  un  jet  de  comburant  froid  qui  est  de  l'air  et  un  jet  parallele  de  combustible  (hydrog^ne)  plac£ 
dans  l'axe  du  precedent.  A  l'aide  de  methodes  optiques,  de  fil  chaud,  on  se  propose  d'etudier  : 

a)  le  champ  de  repartition  spatiale  de  la  vitesse  de  1 'ecoulement  comburant 

-  ■£.(*,  .  **■'(-*,  ^ ^  ■  *3 

comprenant  deux  parties  :  une  valeur  moyenne  ax.  et  une  valeur  fluctuante  ^  . 

b)  le  champ  des  masses  volumiques  qui  fluctue  par  suite  de  la  turbulence  des  deux  ycoulements  :  gaz  combu¬ 
rant  d’une  part  et  gaz  combustible  injecte  dans  le  gaz  comburant  d’autre  part  et  qui  se  mflangent,  ce 
dernier  provoquant  une  variation  de  masse  volumique 

?(J  -  p  ('*.  ♦  (°e  ('w  ^  ) 

c)  le  champ  des  masses  volumiques  variables  dans  l'espace  par  suite  du  melange  des  produits  de  combustion 
avec  les  produits  frais  comburant  et  combustible 

wf  - fr('x  ^  l 
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d)  enfin  le  champ  des  variations  dans  l'espace  de  la  masse  volumique  dues  aux  reactions  chimiques  de 
combustion  qui  liberent  de  l'energie 

?/--  fr<*. 

Ces  champs  elementaires  peuvent  Itre  etudies  separement  et  leur  structure  peut  etre  d^taillee  S  par- 
tir  des  fonctions  de  correlation  et  des  spectres  energetiques  exprimes  en  fonction  du  nombre  d'onde  [1J. 

Ces  fonctions  peuvent  expliquer  le  mecanisme  de  transport  et  de  dissipation  en  consid€rant  la  super¬ 
position  et  1 ' interaction  de  ces  champs  elementaires.  Une  fa$on  d'operer  consiste  &  ytudier  separement  la 
structure  de  chaque  champ  elementaire  et  ensuite  la  structure  complexe  du  champ  resultant.  II  est  tout 
d'abord  important  de  connaitre  le  champ  puisque  les  quantites  scalaires  concentration-temperature 

sont  transportees  par  la  turbulence  due  aux  fluctuations  de  vitesse.  Dans  ce  cas  particulier  cette  analyse 
montre  que  la  turbulence  est  isotrope  et  obeit  &  la  loi  de  dissipation  visqueuse  de  Kolmogorof.  Le  champ 
de  melange  &  froid  du  jet  d' hydrogene  dans  un  ecoulement  d'air  peut  etre  etudie  par  un  resultat 
important  fournU  par  1' analyse  des  spectres  energetiques,  corame  ceux  represent^s  sur  la  figure  1,  ainsi 
que  le  champ  dans  la  zone  de  combustion.  On  note  en  particulier  sur  cette  figure  la  difference  des 

spectres  avec  et  sans  combustion.  Remarquons  tout  d'abord  que  les  petites  structures  correspondant  U  des 
grandes  valeurs  du  nombre  d'onde  Ak  ont  des  pentes  differentes  de  la  loi  de  Kolmogorof  :  (-  3)  au  lieu 
de  C-  5/3),  la  dissipation  visqueuse  est  plus  ^brupte  pour  les  scalaires  que  celle  correspondant  &  la 
vitesse.  Le  champ  du  &  la  reaction  chimique  5^  a  une  structure  tres  diff€rente  des  autres  champs 
(  ,  ?r  ,  ?c  )  .  Le  spectre  presente  un  maximum  autour  de  1  mm  pour  k“l.  Cela  indiquerait  des  paquets 

fonctionnant  comme  des  micro-foyers  et  une  stabilite  accrue  de  zones  de  reaction  chimique  de  cette  dimen¬ 
sion.  Cette  conception  exclut  done  un  mecanisme  de  combustion  suivant  un  front  de  flamme.  L' analyse  du 
bilan  d'energie  turbulente  defini  comme  1' integrate  par  rapport  5  toutes  les  frequences  du  spectre  £ner- 
getique  fournit  des  renseignements  complementaires  sur  les  mecanismes  d'echange  dans  les  flammes.  La  loi 
de  repartition  de  l'energie  dans  la  flamme  d'apres  les  nombres  d'onde  et  la  diminution  de  l'energie  totale 
peut  resulter  de  la  diminution  du  poids  des  coefficients  de  Fourier  pour  des  frequences  les  nlus  basses 
(  K<  0,3  mnfl)  et  ce  sont  les  coefficients  de  Fourier  qui  sont  situes  autour  de  K  =  0,9  mnfl  qui  ont  le 
plus  grand  poids. 

La  figure  2  montre  egalement  que  le  melange  hydrogene-air  &  froid  ^E  (  K))  et  le  melange  gaz  frais 
gaz  brules  ont  des  spectres  identiques  ainsi  que  le  spectre  global  E  c*  C K)  avec  deux  pentes  :  l'une 
relative  &  K  pour  les  grandes  structures  et  une  autre  pente  relative  &  K  pour  les  petites 

structures  montrant  un  mecanisme  identique  pour  le  melange  hydrogene-air  et  le  melange  gaz  frais-produits 
brules . 

Nous  pensons  que  des  resultats  de  ce  type  sont  tres  utiles  pour  etudier  la  structure  de  1 ' ecoulement , 
les  processus  de  melange  et  les  processus  de  combustion. 

Notre  deuxieme  exeraple  choisi,  montrant  tout  l'interet  de  ces  mesures  scalaires,  concerne  1' etude  des 
mecanismes  de  transport  et  de  la  structure  du  champ  scalaire  dans  le  cas  d'une  couche  limite  de  melange 
entre  un  jet  d'air  chaud  et  le  milieu  ambiant  au  repos,  etude  effectuee  par  Fiedler  2  1 'University  de 
Berlin  [2].  Une  analyse  du  champ  de  vitesse  et  de  temperature  permet  d'evaluer  dans  l'equation  du  bilan 
des  fluctuations  de  temperature  (  T' )  dans  le  champ  de  1' ecoulement  en  faisant  apparaitre  les  termes 
de  convection,  de  production,  de  diffusion  et  de  dissipation  comme  l'indique  la  figure  3  en  donnant  en 

M  .V 

fonction  du  rapport  1}  r  sens  de  1 ' ecoulement ,  ^  perpendiculaire  £  1 ' ycoulement)  les  valeurs 

des  processus  qui  apportent  un  gain  ou  qui  enregistrent  une  perte  :  la  production  est  le  composant  le  plus 
important,  il  possede  deux  maxima,  la  convection  et  la  diffusion  sont  presque  symetriques,  la  convection 
est  nulle  au  centre  de  la  couche  de  melange  oh  la  diffusion  atteint  un  maximum,  la  dissipation  n'est  pas 
mesuree  et  est  obtenue  par  difference.  Cette  etude  permet  d'apprecier  la  distribution,  les  termes 

V  T  y  "f"  *  correspondant  &  la  diffusivite  laterale  et  longitudinale  et  V'/  "T  ^  ,  Le  processus  de 

melange  a  egalement  ete  analyse  en  perturbant  1' ecoulement  par  des  ondes  periodiques  issues  d'un  haut 
parleur  place  au  bord  de  fuite  de  la  tuyere.  La  conclusion  des  auteurs  est  ; 

-  il  existe  une  forte  influence  des  perturbations  externes  de  faible  amplitude  sur  la  structure  de  l'€cou- 
lement,  et  tout  particulierement  sur  la  structure  scalaire,  se  traduisant  par  la  formulation  de  cellule 
d  structure  coherente  ; 

-  dans  le  cas  de  1' ecoulement  non  perturbe,  les  structures  coherentes  semblent  plutot  peu  importantes  ; 

-  les  regions  de  production  negatives  sont  presentees  dans  le  cas  de  1 'ecoulement  excite  ; 

-  le  moddle  du  gradient  de  diffusion  n'est  pas  applicable  ; 

-  1 ' excitation,  non  seulement  discipline  les  gros  tourbillons  (vortices),  mais  aussi  influence  la  struc¬ 
ture  de  la  turbulence  interne. 

Ce  travail  il lustre  bien  ce  que  l'on  peut  deduire  de  l'£tude  d'un  champ  de  scalaire  couple  avec 
l'etude  d'un  champ  de  vitesse  pour  modeliser  la  structure  de  l'ycoulement  et  en  deduire  les  processus  de 
transport . 

Le  troisieme  exemple  que  nous  avons  choisi  concerne  la  mise  en  £vidence  des  structures  coherentes  en 
combustion  par  des  techniques  photographiques  qui  traduisent  en  fait  le  champ  des  masses  volumiques.  Dans 
ce  domaine  un  travail  trds  interessant  a  ete  fait  par  Yule,  Chigier  et  Thompson  T3].  Dans  bien  des  cas  le 
role  des  structures  coherentes  est  plus  facilement  obtenu  dans  les  zones  de  combustion  par  suite  du 
caractdre  lumineux  des  flammes,  e'est  en  particulier  une  faqon  simple  de  visualiser  le  melange  combus- 
t ible/comburant .  Les  flammes  de  diffusion  turbulentes  offrent  une  grande  variety  de  structures  et  dans  la 
plupart  des  cas  des  anneaux  tourbillonnaires  apparaissent  analogues  *1  ceux  mis  en  yvidence  dans  les 
melanges  non  reactifs,  le  confinement  de  la  flamme  peut  changer  cette  structure,  la  longueur  de  la  zone  oh 
se  rencontrent  les  anneaux  tourbillonnaires  est  plus  grande  avec  combustion  que  sans  combustion.  On  note 
une  coalescence  de  ces  anneaux  dans  les  deux  situations  avec  et  sans  combustion.  Ces  anneaux  apparaissent 
egalement  dans  la  combustion  des  jets  de  combustible  pulvyrisy  et  1 'existence  de  ces  remous  cohyrents 
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depend  de  la  distribution  du  melange  comburant-combustible.  La  figure  4  donne  un  schema  de  remous  appa- 
raissant  dans  une  flamme  de  diffusion  et  separant  le  combustible  du  comburant.  La  figure  4  represente  par 
exemple  une  coupe  d'un  anneau  tourbillonnaire  ou  dans  un  systeme  bidimensionnel  etroit  une  ligne  de 
tourbillons  ou  une  coupe  d'un  tourbillon  helicoidal.  Ces  tourbillons  sont  analogues  et  s'enroulent  comme 
ceux  observes  dans  un  sillage  mais  ils  ne  sont  pas  alternatifs  mais  sym^triques  [4].  Cette  schematisation 
correspond  &  un  tourbillon  qui  n'en  a  pas  rencontre  d'autres,  ce  qui  provoque  des  instabilit£s[5J.  La  zone 
delimitant  le  combustible  de  l'air  a  la  forme  d'une  double  spirale  entourant  un  noyau  tourbillonnaire.  Une 
flamme  se  situe  S  la  lisiere  de  cette  double  spirale  separant  le  combustible  de  l'air  et  le  noyau  tourbil¬ 
lonnaire  contient  du  combustible  et  de  l'air  entraind  et  peut  constituer  une  zone  de  combustion  volumique 
£  richesse  variable  fonctionnant  comme  une  zone  de  recirculation  et  la  flamme  en  bordure  des  deux  zones 
peut  etre  consid^ree  comme  une  flamme  de  diffusion  laminaire. 

Cette  structure  peut  done  £tre  assez  bien  decrite  si  l'on  connait  1 'evolution  de  la  richesse,  le 
volume  de  cette  poche  augmente  avec  1 'accroissement  du  temps  de  sejour  dans  le  tourbillon.  Dans  la  situa¬ 
tion  de  la  figure  4  une  flamme  entoure  completement  le  tourbillon.  Ces  tourbillons  peuvent  rencontrer 
d'autres  tourbillons,  s'unir  pour  former  des  structures  plus  complexes  &  caractere  tridimensionnel ,  ou 
apres  coalescence  de  plusieurs  tourbillons  se  restructurer . 

Un  autre  type  de  structure  coherente  observe  dans  la  flamme  est  constitue  par  de  grands  tourbillons 
allonges  tridimensionnels  dont  la  structure  est  donnee  sur  la  figure  5  qui  met  en  evidence  une  incomplete 
coherence  circonferent ielle  autour  de  la  flamme  et  aussi  l'existence  d'une  distribution  de  tourbillons 
irreguliers  avec  la  presence  de  tourbillons  d'echelle  plus  petite  avec  des  poches  de  comburant,  de  combus¬ 
tible  et  de  produits  brGles.  Le  combustible  et  l'air  s ' engouf frent  dans  ce  tourbillon  et  se  melangent  plus 
rapidement  que  dans  la  situation  precedente.  Cette  structure  est  observee  dans  de  nombreuses  flammes  de 
diffusion.  Ces  trois  resultats  ont  ete  choisis  pour  montrer  ce  qui  peut  etre  deduit  de  1 'analyse  des 
champs  scalaires.  Les  fluctuations  de  pression  ont  ete  negligees  car  la  mesure  de  ce  parametre  est  beau- 
coup  plus  difficile. 

Decider  de  ce  qui  est  coherent  de  ce  qui  ne  l'est  pas,  comprendre  les  mecanismes  de  transport,  recon- 
naitre  les  modifications  apportees  &  la  structure  de  la  zone  turbulente  par  la  combustion,  voir  1 'action 
d'ondes  de  pression  sur  1 'organisation  des  phenomenes  et  en  particulier  sur  la  convection,  comprendre  les 
phenomenes  de  dissipation,  en  un  mot  apporter  des  informations  pour  modeliser  la  turbulence,  tels  sont  les 
resultats  que  l'on  peut  attendre  des  methodes  de  diagnostic  en  general  et  des  methodes  de  mesure  des 
scalaires  en  particulier. 

3  -  METHOVES  VE  MESURE  VE  LA  TEMPERATURE 


II  existe  de  nombreuses  methodes  pour  determiner  la  temperature  dans  un  ecoulement  turbulent  sans  ou 
avec  combustion,  mais  il  existe  deux  grandes  classes  :  les  methodes  optiques  et  les  methodes  par  sondes, 
e'est  pourquoi  nous  decrivons  une  methode  dans  chaque  classe  pour  voir  apparaitre  les  avantages  et  les 
inconvenients  de  chacune  d'elles. 

3.  /  -  Methods  optiquoA 

Ces  methodes  of frent  l'interet  de  ne  pas  perturber  1 ' ecoulement ,  ce  qui  est  un  gros  avantage  dans  le 
cas  d'ecoulement  turbulent,  mais  elle  integre  sur  une  partie  de  la  longueur  du  faisceau  les  phenomenes.  On 
peut  cependant  remedier  en  partie  &  cet  inconvenient  en  introduisant  localement  dans  le  gaz  des  traceurs 
qui  ont  la  temperature  mais  aussi  la  vitesse  locale  du  gaz.  La  technique  que  nous  allons  decrire  est  celle 
du  renversement  de  raies  mais  en  se  plagant  dans  le  domaine  infra-rouge,  ce  qui  permet  de  faire  des  mesures 
a  partir  de  600  K. 

3.1.1  -  MteuAeA_iMpinfaneeA_de_Za_temi?e{La£u/ie._de._ga.z_de_cqrribu6tiqn_e£  de._la 
^^§q4^ce._tRcAmZqui_paAplfAqmitA2c._Zn^a2^quge'l6i~ 

A)  -  P/Uncsipe.  de  la  methodz  de  meAu/io. 

Le  schema  optique  du  dispositif  qui  permet  de  discriminer  les  signaux  provenant  de  la  source  et  de  la 
flamme  est  presente  figure  6.  L' image  de  la  source  S  est  formee  dans  le  plan  C  par  la  lentille  L-d  .  La 
lentil le  L  £  en  reforme  une  image  dans  le  plan  de  la  flamme.  L' optique  qui  supporte  le  diaphragme 

d'ouverture,  conjugue  les  plans  images  precedents  sur  le  diaphragme  de  champ  di  .  La  lentille  Li,  permet  le 
transport  de  l'image  monochromat ique ,  isolee  par  le  monochroma teur /'V.  ,  sur  le  miroir  prismatique  P  qui 
la  divise  en  deux  parties  sensiblement  egales.  Grace  aux  lentilles  l_5  une  moitie  de  l'€nergie  arrive  sur 
un  recepteur  ,  1' autre  moiti£  sur  un  recepteur  R t  .  Un  couteau  plac€  dans  le  plan  C  peut  etre  d§place 
de  fagon  a  diviser  l'image  exactement  comme  le  miroir  prismatique.  L'un  des  rScepteurs,  par  exemple  R^, 
ne  regoit  alors  plus  aucune  energie  provenant  de  la  source,  alors  que  le  signal  d^livre  par  n'a  pas 
change . 

Er.  l'absence  de  flamme  et  avant  obturation  par  le  couteau  les  ddtecteurs  R  ^  et  delivrent  respecti- 
vement  les  signaux  proport ionnels  d  la  luminance  energetique  monochromat ique  de  la  source 

oD  est  la  temperature  de  luminance  de  la  source  3  la  longueur  d'onde  A  .  On  peut  aussi  faire  appa¬ 

raitre  son  emissivit£  £  (A}  en  ecrivant  : 

kJ5(a)  L(a/T5)  ,  fs(A)L(A,T5) 

oD  est  alors  sa  temperature  vraie.  Compte  tenu  des  differences  de  sensibilite  des  detecteurs  et  des 

etendues  geometriques  des  faisceaux  incidents,  les  coefficients  de  proport ionnal ite  et  £ ne  sont  pas 
identiques  mais  il  est  toujours  possible  de  les  £galer  en  jouant  sur  les  Stages  d'amplif ication. 
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En  presence  de  flamme  et  couteau  mis,  donne  un  signal  du  I  la  flamme  seule  alors  que  R^ 

_  „ : _ i  >  i  »  c-> _ s  i  «  ♦  " 


donne  un  signal  .  du  h  la  fiamme  et  a  la  source  . 

i*  *  ,  .♦  * 

-J  <  (.  A,x)  Ooc 


r  l  a,  rts)  *  -< 


-j*  (A, 


j  <  O,  *0  L  [A ,  Tf  (*)]  -e 


-J  o(  (A,  x  )  cCk 


(he 


(he 


est  coefficient  d'Emission  monochromatique  local  Egal  au  coefficient  d'absorption  d'apres 
la  loi  de  Kirchhoff  et  L  (T^,  ”^(*)^]  est  la  luminance  monochromatique  du  corps  noir  &  la  temperature 
locale  de  la  flamme  TJi'r)  • 

Dans  cette  formulation  on  a  suppose  1 'homogeneite  de  la  flamme  dans  une  coupe  transversale  du 
faisceau,  ou  ce  qui  revient  au  meme,  que  le  faisceau  est  infiniment  etroit.  Le  premier  terme  de  corres¬ 
pond  h  la  transmission  h  travers  la  flamme  de  l'Energie  provenant  de  la  source  selon  la  loi  de  Beer,  loi 
rigoureuse  dans  le  cas  strictement  monochromatique.  Le  deuxieme  terme  correspond  &  1' emission  de  la  flamme 
dans  une  tranche  d’epaisseur  djx.  et  h  son  absorption  par  les  tranches  suivantes  de  x  &  ♦  -i  ,  t-i  etant 
l'epaisseur  de  la  flamme.  Seul  ce  second  terme  apparait  dans  le  signal  .  Si  la  flamme  est  homogene 
1' integrate  est  simple  car  L(A>T<.  )  sort  de  l'integrale  et  1 'homogeneite  du  coefficient  d'emission 
n'est  meme  pas  necessaire.  On  obtient  les  relations  suivantes  : 

N),  =  L(A,TuS)(^_e)  +  ^  L(A,  Tf  )  £ 


oD  £  est  le  facteur  d'emission  defini  par 


£(*) 


-  J  o<  (  A  y  oc)  he 


- 


On  remarquera  que  £  depend  de  l'epaisseur  de  la  flamme.  La  connaissance  de  la  luminance  monochro- 
raatique  de  la  source  L(  permet  alors  de  calculer  h  chaque  instant  par  la  resolution  de  ce  systeme 

de  deux  equations,  les  inconnues  L(A,T^  }  et  €  (A) 

Pour  effectuer  cette  resolution,  dans  le  cas  oil  =  R 

on  pose  : 

-- 

°n  °b£ient  8  ^  .  K  €  L(A,Tu5) 


et  on  en  deduit  : 


l  c  ^ ,  Tf )  _  y& 

L  (  A,  TLS)  "  ys 


La  loi  de  Planck  pour  la  longueur  d'onde  A  permet  alors  de  remonter  a  la  temperature  de  la  flamme 
Tj;  .  Le  temps  de  reponse  n'est  limite  que  par  la  constante  de  temps  des  recepteurs  et  celle  de  l'enre- 
gistrement  des  signaux. 


Le  calcul  d'erreur  montre  que  la  precision  de  la  mesure  est  d'autant  plus  grande  : 

-  que  l'on  travaille  au  plus  pres  du  point  de  renversement  defini  par  l'egalite  des  temperatures  de 
la  flamme  et  de  la  source, 

-  et  que  le  facteur  d'emission  de  la  flamme  £  est  proche  de  1 'unite. 

II  faut  souligner  encore  une  fois  que  la  mesure  n'a  de  signification  claire  que  si  la  temperature  de 
la  flamme  est  uniforme  sur  tout  le  volume  gazeux  traverse  par  les  faisceaux  optiques  et  que  la  mesure 
interesse  tout  ce  volume.  En  general  la  flamme  presente  un  profil  de  temperature  avec  un  maximum  au  milieu 
du  jet.  Dans  ce  cas  la  mesure  conduit  soit  h  une  temperature  proche  de  la  temperature  maximale  si  le 
facteur  d'emission  est  faible,  soit  £  une  temperature  voisine  de  celle  en  bordure  du  jet  si  le  facteur 
d'emission  est  proche  de  1.  En  fait  lorsque  l'on  s'interesse  plus  particulierement  h  la  temperature  au 
centre  du  jet  on  montre  qu'il  existe  un  optimum  pour  le  facteur  d'emission,  optimum  qui  est  trEs  bien  mis 
en  evidence  par  1' etude  spectrographique  d'une  raie  d' emission-absorption  qui  prEsente  deux  maximums 
symetriques  par  rapport  au  centre  de  la  raie. 

Dans  le  cas  d'une  flamme  transparente  1' addition  d'atomes  etrangers  (du  sodium  par  exemple)  emettant 
des  radiations  definies,  est  necessaire  afin  de  la  rendre  monochromat iquement  lumineuse  ;  1 'hEtErogEnEitE 
des  temperatures  du  gaz  peut  n' avoir  alors  aucune  influence  sur  le  resultat  de  la  mesure  si  l'on  a  soin 
d'injecter  le  sel  de  sodium  localement.  II  faut  par  contre  prendre  garde  que  celui-ci  ait  atteint  l'Equi- 
libre  thermique. 

Bien  evidemment  la  methode  suppose  l'Equilibre  thermodynamique .  Dans  le  cas  contraire  la  loi  de 
Kirchhoff  doit  etre  modifiee  de  la  fa£on  suivante  : 

la  luminance  energetique  spectrale  d'un  gaz  quelconque  est  Egale  au  produit  de  son  facteur  d'absorp¬ 
tion  par  la  luminance  Energetique  spectrale  L  A  du  corps  noir  £  la  temperature  dEfinie  par  la 

loi  de  Boltzmann  Ecrite  pour  les  populations  et  responsables  de  1 'Emission  et  de  1 'absorption  de 

la  raie  de  longueur  d'onde  A 

_SiL  -  J .  J!±—  , 

*T  * 

oh  et  sont  les  poids  statistiques  respectifs  des  niveaux  1  et  2  de  l'atome  ou  de  la  molEcule 

Emetteur.  On  effectue  done  une  mesure  de  tempErature  de  population.  Dans  le  cas  courant  oh  il  existe  des 
Equilibres  partiels  correspondant  aux  transitions  Electroniques ,  de  vibrations  ou  de  rotations,  la  mesure 
conduit  h  la  tempErature  d'Equiliore  caractErist ique  du  type  de  transition  utilisE.  Par  exemple  si  l'on 
travaille  sur  l'une  des  raies  du  doublet  du  sodium  on  mesurera  une  tempErature  Electronique . 
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Enfin  la  methode  suppose  que  les  phenomenes  d'dmission  et  absorption  l'emportent  devant  tous  les 
phdnomSnes  de  diffusion  qui  sont  enregistrds  comme  une  perte  de  signal.  En  particulier,  il  est  possible 
de  travailler  sur  les  particules  en  suspension  dans  le  gaz  (flammes  Sclairantes)  d  condition  que  leur 
dimension  soit  trds  petite  devant  la  longueur  d'onde  utilisee. 

8)  -  T/umipoi-ction  de  la.  mitkode  en  initia-Muqe 


L'appareil  est  constitue  de  deux  parties  :  un  emetteur  et  un  recepteur  dont  l'alignement  optique  doit 
etre  aussi  parfait  que  possible.  Cet  appareil  utilise  l'absorption  et  l'emission  propre  des  flammes  dues 
essentiellement  d  la  vapeur  d’eau  et  au  gaz  carbonique.  L'intervalle  spectral  y  est  alors  limits,  non  plus 
par  un  monochromateur,  mais  par  un  filtre  interferentiel  qui  integre  sur  la  structure  fine  d'une  bande  de 
vibration-rotation.  On  commet  ainsi  une  erreur  systematique,  le  calcul  monochromatique  n'dtant  plus  rigou- 
reusement  exact.  Nous  allons  montrer  que  c'est  une  approximation  gengralement  suffisante.  L'homogdniite  en 
temperature  de  la  flamme  etant  admise,  1 ' integration  des  relations  donnant  y  et  y^  sur  l'intervalle 
spectral  A  A  =  A^  -  A^  donne  :  ^ 

yt  =  KjlL(\T()£(\)dA  y  y}  =  «J  e£(Mua.ts) 

^4  b  -4 

Les  variations  d'emissivite  du  gaz  €  CA)  en  fonction  de  la  longueur  d'onde  sont  generalement 
tres  fortes,  par  contre  la  luminance  du  corps  noir  pour  les  temperatures  Tr  et  Ts  et  l'emissivite  de 


tres  fortes,  par  contre  la  luminance  du  corps  noir  pour  les  temperatures  et  Ts  et  l'emissivite  de 

la  source  £ s  peuvent  Stre  considerees  comme  constantes  et  sortir  des  integrales  si  l'intervalle  spec¬ 
tral  est  suffisamment  etroit.  Dans  ces  conditions  la  resolution  du  systlme  peut  etre  conduite  comme  dans 
le  cas  monochromatique  : 


J*  =  L  (  Xm.  ,  T,  ^ 

^5  £s(VH(V,Ts) 
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oO  est  une  longueur  d'onde  moyenne  de  l'intervalle  .  L'emissivite  ainsi  obtenue  est  evidemment 

une  Imissivite  moyenne  £  ,  et  elle  ne  sera  importante  que  si  la  densit€  des  raies  spectrales  est  assez 

forte. 

Afin  d'estimer  1 ' importance  de  1 'erreur  systematique  due  2  cette  integration,  on  considere  le  cas  oQ 
la  variation  de  luminance  du  corps  noir  sur  aA  est  lineaire 

L  (  \  T  )  —  L(XWjTK  —  (  V.,  T  X  A  ~  ) 
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Puisque  f^A)  est  une  fonction  bornee  par 

f  ^ 


le  terme  correctif  peut  etre  major!  par  : 

Ik  1 


soit  en  utilisant  la  loi  de  Planck 


T<  ) 


^-"p-tV,  ) 


L'erreur  seta  negligeable  si  l'intervalle  spectral  est  suffisamment  Etroit  et  si  l'emissivite  moyenne 
de  la  flamme  n'est  pas  trop  faible.  Elle  sera  quasi-nulle  si  la  temperature  vraie  de  la  source  est  Igale  S 
celle  de  la  flamme. 


L'etude  spectrometrique  prealable  des  flammes  de  combustion  montre  que  l'emissivite,  par  ailleurs 
faible,  est  par  contre  tr!s  forte  A  la  longueur  d'onde  de  4,5/tmsur  les  bandes  chaudes  d'emission  du  gaz 
carbonique.  C'est  ce  que  montre  le  spectre  obtenu  sur  un  r!ac£eur  homoglne  (Fig.  7),  l'lpaisseur  de  gaz 
traversee  etant  de  quelques  centimetres.  A  cette  longueur  d'onde  on  utilise  un  filtre  interferentiel  dont 
la  largeur  relative  est  de  4%.  Dans  ces  conditions  l'erreur  systematique  sur  la  luminance  L  ) 

est  inflrieure  au  pour  cent. 

Le  temps  de  rlponse  de  1' ensemble  dltecteur-prlamplif icateur  est  de  l'ordre  de  une  micro-seconde , 
mais  il  ne  passe  pas  le  continu  et  il  est  nlcessaire  d ' interrompre  plriodiquement  le  flux  incident  pour 
recaler  le  zero.  On  utilise  la  rotation  d'une  roue  dent!e  qui  occulte  1' image  obtenue  dans  le  plan  d. 

(Fig.  6)  3  une  frequence  de  2500  Hz.  Les  deux  diaphragmes  de  champ  qui  delimitent  les  faisceaux  1  et  2 
sont  constituls  par  les  surfaces  sensibles  des  detecteurs,  le  diaphragme  d'ouverture  !tant  commun  aux  deux 
faisceaux.  Cette  separation  des  deux  diaphragmes  de  champ  a  pour  avantage  de  permettre  par  le  seul  d!pla- 
cement  axial  du  miroir  separateur  de  disjoindre  les  deux  images,  ce  qui  rend  plus  facile  l'occultation  de 
l'une  d'elles  par  le  couteau. 


La  photo  de  la  figure  8  montre  l'appareil  sur  un  banc  d' optique  et  au  premier  plan  le  calculateur 
analogique  qui  permet  de  calculer  en  temps  riel  la  luminance  et  le  facteur  d'emission  de  la  flamme.  Une 
difficult!  particulidre  pour  effectuer  ce  calcul  est  due  h  la  modulation  &  2500  Hz.  Les  signaux  se 
trouvent  haches  &  cette  frequence  avec  par  ailleurs  un  leger  dlphasage  entre  les  deux  dltecteurs,  ce  qui 
fait  que  le  rlsultat  du  calcul  n'est  valable  que  pendant  de  courts  intervalles  de  temps  de  &0jus  espacls 
de  300/js  environ.  On  proclde  done  par  Ichantillonnage  et  la  chaine  de  traitement  logique  a  l&rge  bande 
passante  raise  au  point  ne  retient  que  les  mesures  correctes  et  affiche  z!ro  peur  les  autres.  Cependant, 
pour  les  phenomlnes  fugitifs  dont  la  dur!e  est  inflrieure  d  quelques  centaines  de  microsecondes  la  modu¬ 
lation  est  inutile. 


C)  -  MeatiAe  doj>  fluctuations  de  tejnpeAaXuAc 

^^oupla.Q^_d<^  _^ZucMw^qn6  _d^_tm^^piAe:_2X_d,_^m^66Zvit^ 


Nous  avons  vu  que  la  mesure  n'a  de  signification  claire  que  si  la  temperature  de  la  flamme  est  uni¬ 
forme  sur  tout  le  volume  gazeux  traverse  par  les  faisceaux  optiques.  Cela  seroble  exclure  la  possibility 
d'effectuer  des  mesures  de  turbulence  thermique  lorsque  toutes  les  caracteristiques  du  milieu  (concentra¬ 
tions,  temperature)  varient  en  tout  point  et  3  tous  instants  de  fa^on  aleatoire.  Nous  allons  montrer  qu'il 
est  pourtant  possible  de  caracteriser  la  turbulence  dans  le  cas  od  elle  est  statistiquement  homogene  et 
lorsque  les  fluctuations  d'emissivite  sont  de  faible  amplitude  vis-d-vis  de  la  valeur  moyenne. 

Pour  le  demontrer  reprenons  la  formule  donnant  l'energie  reque  par  un  detecteur  en  provenance  d'une 
source  apres  absorption  par  la  flamme  et  en  provenance  de  la  flamme  elle-meme  et  ceci  pour  un  faisceau 
optique  monochromatique  tres  etroit  :  f  i 

.  I  «  cKoc  i  ~  J  o'* 

y  —  s  *  ’*  +  J  °<  L  €  cine 

-l  * 

Si  l'on  separe  valeur  moyenne  et  valeur  fluctuante  : 

t  .  /  .  # 

seules  les  fonctions  **  et  L,  dependent  de  oc  et  de  t  puisque  le  milieu  est  suppose  statistiquement 
homogene .  *  ^ 

On  en  deduit  la  meme  separation  oour  le  facteur  d’emission  avec  :  /  '  chx. 

I  ,*-*-£**  ,  ret)*-*1 

et  la  partie  fluctuante  du  signal  ^/tO  s'ecrit  :  *■£ 

,  ,  **  .J  («■  + «')  <u 

y<  Z  chx. 
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On  remarque  que, si  la  luminance  de  la  source  est  egale  &  la  luminance  moyenne  de  la  flamme,  les 
fluctuations  d'emissivite  disparaissent_au  premier  ordre.  Si  ces  fluctuations  sont  de  faible  amplitude  il 
est  possible  de  negliger  <x  '  Levant  et  il  ne  reste  que  les  fluctuations  liees  &  la  temperature. 

'  —  I  ,  -  °<  ( -  -x) 
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L’ interpretation  de  cette  integrale  sera  encore  facilitee  si  on  fait  en  sorte  que  le  coefficient  d 'absorp¬ 
tion  soit  petit  de  faqon  a  pouvoir  negliger  o<  (  .  /*« \  devant  1  : 


o<  [  C  chx. 
<4  f 


Pour  cela  on  se  place  sur  la  bande  d'emission  &  2,5^m  pour  laquelle  le  facteur  d'emission  est  tr£s  faible 
comme  le  montre  le  spectre  obtenu  sur  un  foyer  homogene  du  type  Longwell  (Fig.  7). 

Enfin,  si  les  fluctuations  de  temperature  sont  de  suf f isamment  faible  amplitude  il  est  possible  de 
lineariser  la  loi  de  luminance  L(^T  )  au  voisinage  de  la  temperature  moyenne  7^  ,  et  obtenir 

ainsi  1' integrale  des  fluctuations  de  temperature  sur  le  trajet  optique  : 

W 

SCO  t  2  i  2i  (A,  T  )  ,  _  j  T  (ocy  t)  olx. 
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avec  c  X  o<  z  c  le  facteur  d’emission  moyen. 

Dans  le  cas  oQ  le  milieu  n'est  pas  homogdne,  c'est-l-dire  lorsque  et  ne  sont  pas  independants 

de  oc  ,  le  decouplage  au  premier  ordre  des  fluctuations  n'est  plus  total  mais  les  fluctuations  d'emissivite 
sont  fortement  attenuees.  Le  terme  relatif  H  ces  fluctuations  s'ecrit  en  effet  : 


J  <*’(*,  t)  (  ty*)  -  b$ 


et  est  tr£s  inferieur  h  celui  que  l'on  aurait  on  absence  de  source  (  ^  *0),  D'ailleurs,  meme  lorsque  le 
milieu  est  homogene,  la  qualite  du  decouplage  est  fonction  de  la  precision  de  l'ajustement  de  la  temperature 
de  la  source  &  la  temperature  moyenne  de  la  flamme  qu'il  est  necessaire  de  determiner  par  une  mesure  ou  un 
calcul  prealable. 

Si  le  decouplage  est  parfait,  le  signal  pyrometrique  est  proport ionnel  &  1' integrale  sur  le  traiet 
optique  des  fluctuations  de  luminance  ou  de  temperature.  Quels  sont  les  effets  de  cette  integration  et 
notamment  quelle  relation  existe-t-il  entre  la  variance  (ou  1 'ecart-type)  du  signal  et  la  variance  de  la 
grandeur  mesuree  ?  On  montre,  en  utilisant  un  formalisme  du  &  Uberoi  et  Kovasznay  [7],  que  lorsque  la 
longueur  d' integration  est  grande  devant  l'echelle  integrale  elles  sont  proport ionnel les .  Le  rapport  de 
proportionnalite  est  trouve  egal  au  rapport  du  double  de  l'echelle  integrale  de  la  turbulence  A-  S  la 
longueur  d ' integration  il  .  C'est-3-dire  que  si  SL'  represente  les  fluctuations  mesurees  et  T*  les  fluc¬ 
tuations  locales,  leurs  variances  JL*  et  T'*  sont  liees  par  la  relation  suivante  : 
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Jl' 1  =  T  .  *  A  . 

a 

Cette  relation  peut  etre  interpretee  de  la  faqon  suivante  ;  si  on  d§coupe  le  trajet  optique  en 
interval les  de  longueur  2  ^ 

/Y\.  - 


•?  A 


on  peut 
ces  intervalles 


O  -r 1 

considerer  que  J(.  est  la  moyenne  de  /vu  grandeurs  independantes  attaches  &  chacun  de 

T<  -  V* . > 


La  variance  de  la  somme  est  alors  egale  &  la  somme  des  variances  et  puisqu'il  s'agit  de  variables  £ 
moyenne  nulle  on  retrouve  la  relation  ci-dessus. 

3.1.2  -  ^^t^rnina^on_de._i,_ech^lle._-Lnt^gn-al^_dc._la_p:iA.bule.nce. 

Pour  remonter  &  la  variance  locale  T  *  de  la  grandeur  mesuree  il  est  necessaire  de  connaitre 
l'echelle  integrale  J\,  .  Pour  l'obtenir  on  peut  penser  utiliser,  soit  1* analyse  spectrale,  soit  la  methode 
des  faisceaux  orthogonaux  mise  en  oeuvre  pour  la  premiere  fois  par  Fischer  et  Krause  l  8]  ,  mais  qui  pre¬ 
sente  1 ' inconvenient  de  necessiter  la  mise  en  oeuvre  de  deux  appareils  de  mesure. 

a]  AnaJS^e.  4pec#ia£e 

Par  analyse  spectrale  on  entend  le  spectre  obtenu  &  partir  du  signal  de  mesure.  Un  tel  spectre  corres¬ 
pond  par  transformee  de  Fourier  I  la  fonction  de  correlation  temporelle.  Par  contre  le  calcul  de  l’echelle 
integrale  ne  peut  se  faire  qu'H  partir  de  la  fonction  de  correlation  spatiale.  Le  passage  de  l'une  & 

1' autre  n'est  possible  que  si  l'on  admet  1 'hypothese  de  Taylor.  Cette  hypothese  consiste  H  supposer  que 
pour  les  grandes  vitesses  de  convection  la  turbulence  est  gelee.  II  y  a  alors  proportionnalite  entre  les 
deux  variables,  le  temps  et  la  coordonnee  spatiale  parallele  &  la  vitesse.  Un  dephasage  dans  le  temps 
correspond  done  3  une  translation  de  cette  coordonnee.  Le  spectre  que  l'on  obtient  alors  est  le  spectre 
dit  "unidimensionnel"  qui  est  la  transformee  de  Fourier  de  la  correlation  entre  deux  points  progressive- 
ment  eloignes  l'un  de  1' autre  definissant  un  vecteur  maintenu  parallele  H  une  direction  donnee,  ici  celle 
de  la  vitesse.  Mais  IS  encore  le  spectre  mesure  est  modifie  par  l'effet  d ' integration.  Dans  le  cas  oil  la 
longueur  d ' integration  est  grande  devant  -A-  ,  on  demontre  la  relation  suivante  entre  le  spectre  integre 
J”^  et  le  spectre  local  £  :  ^ 

JL  I  K  £(*)  U  « 

77  is 

Classiquement  on  calcule  l'echelle  integrale  H  partir  de  la  valeur  II  l'origine  du  spectre  unidimen¬ 
sionnel  :  »  —  /  o 

£  to)  Jj,  T 
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Le  calcul  de  E^(o]  £  partir  du  spectre  permet  d'ecrire  : 

E  m  :lfr(K)<AK  =  A.  SI 

S  >  71  J  ^  71 
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et  on  retrouve  la  relation  reliant  les  variances 


SI 


>1 


’ \ 
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L' analyse  spectrale,  meme  compte  tenu  de  1' hypothese  de  Taylor,  ne  permet  done  pas  d'obtenir  une  rela¬ 
tion  supplementaire  qui  aurait  permis  de  calculer  l'echelle  integrale  et  la  variance  locale.  Rdciproquement 
le  spectre  local  ne  peut  etre  determine  qu'&  un  facteur  pres.  Remarquons  au  passage  l'effet  de  1' integra¬ 
tion  sur  le  spectre  :  si  £  peut  etre  represente  par  une  loi  en  K  ^  -T^  sera  alors  en  ^  j  . 

L'effet  de  1 ' integration  se  fait  sentir  dans  tout  le  spectre  par  une  decroissance  plus  rapide. 

b)  F&idceaux  orthogonaux 

Cette  technique  consiste  &  mesurer  le  niveau  de  correlation  entre  deux  mesures  optiques  sur  deux 
trajets  orthogonaux  et  progressivement  Eloignes  l'un  de  l'autre  comme  le  montre  la  figure  9.  Les  points  JL 
et  delimitent  la  plus  courte  distance  cL  des  deux  droites.  Si  l'on  decoupe  les  deux  trajets  optiques 
en  intervalles  de  longueur  2  -A.  (en  prenant  A  et  B  comme  point  milieu  de  l'un  des  intervalles),  chacune 
des  mesures  et  Sl'i  est  la  moyenne  des  grandeurs  independantes  attachees  U  chacun  des  elements 

constituant  l'un  des  trajets  optiques.  On  conqoit  que  seuls  les  deux  intervalles  incluant  le  point  -A-  ou 
B  donneront  une  contribution  notable  d  la  correlation  entre  JV^  et  SI  %  puisque  la  distance  entre  tous 
les  points  des  autres  intervalles  sera  superieure  &  l'echelle  integrale A-  .  En  fait,  Fischer  et  Krause 
admettent  que  la  fonction  de  correlation  fl(cA\  que  l'on  obtient  point  par  point  en  eloignant  progressive¬ 
ment  les  faisceaux  l'un  de  l'autre,  est  proportionnelle  3  la  fonction  de  correlation  R  (d.)  entre  les 
points  ft  et  3  .  II  suffit  alors  de  normaliser  /3(A)  par  la  correlation  en  faisceaux  croises  {%(*>) 
pour  obtenir  P(dl)  ♦  puis  par  integration  l'echelle  integrale.  En  fait  1 ' interpretation  de  Fischer  et 
Krause  n'est  pas  correcte  car  elle  suppose  negligeable  les  correlations  entre  les  points  voisins  de  A  et 
3  »  tels  que  C  et  J)  par  exemple.  L'etude  theorique  de  cette  technique  a  done  ete  reprise  en  utili- 
sant  toujours  le  formal isme  dG  &  Uberoi  et  Kovasznay  valable  pour  un  milieu  statistiquement  homogene.  Ce 
calcul,  pour  deux  faisceaux  infiniment  etroits  de  longueur  Z  La  et  2  -^2.  dirigSs  l'un  selon  l’axe  1  et 
l'autre  selon  l'axe  2  conduit  &  la  relation  integrale  suivante  : 


fiU), 
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En 


particulier  la  correlation  en  faisceaux  croises,  pour  d-m  0,  est  donnSe  par 


/5(0)=_4z_  T^j\A 

A  .  ****.*'. 

oQ  est  une  echelle  int^grale  diff£rente  de  celle  utilisee  jusqu'ici  et  definie  par  : 

A  rAj  j  R(Z)tdt 

La  normalisation  de  /3W)  par  file)  permet  de  dlfinir  la  fonction  sans  dimension  <h(dL\ 

MU)  - 
Mo)~7'lXK 


one cion 


dont  1* inversion  donne 


R  (t)--ZT’ljL 


LW)Jd:6 


On  voit  qu'il  est  necessaire  pour  obtenir  R  de  tracer  P  en  fonction  du  parametre  au  lieu  de 
et  d'en  prendre  la  d€rivee.  La  valeur  de  cette  derivee  S  l'origine  permet  de  calculer  J\_  : 

Itycd*) 


'Z(JL') 


d  1  0 


— ^ 
t-k- 


Enfin  1  *  integration  de  cette  fonction  derivee  sur  cL  permet  d'obtenir  une  relation  entre  -A  et  A 
et  done  de  calculer  _A-  :  ^  .*o 


La  methode  des  faisceaux  orthogonaux  permet  done  bien  de  determiner  1 'echelle  int^grale  ;  elle 
permet,  mieux  encore,  d'obtenir  la  fonction  de  correlation  5  deux  points,  mais  a  pour  inconvenient  d'etre 
longue  (il  faut  effectuer  des  mesures  point  par  point)  et  assez  peu  precise  (la  fonction  R  s'obtient  par 
derivation  de  la  fonction  •  Aussi  on  a  recherche  et  on  propose  une  mesure  plus  simple  derivee  de 
celle-ci.  Elle  consiste  A  mesurer  la  correlation  normalisee  en  faisceaux  croises. 


c)  CoyizlaXion  noAmatl&ez_en  ficuAczaux  c/iql&H 

Il  est  toujours  difficile  d'effectuer  des  mesures  absolues  et  ceci  est  egalement  vrai  pour  les  corre¬ 
lations.  Ce  que  donne  le  traitement  numerique  e'est  la  correlation  normalisee  par  les  valeurs  quadratiques 
moyennes  des  mesures  et  .  Or  nous  avons  donne  une  expression  de  ces  valeurs  quadratiques 

moyennes  :  _  1$  .  V.  _  AA  _  Vs 

U'j)  =(t"A)  *l  (JA 


La  correlation  en  faisceaux  croises  donnee  par  : 

devient  apres  normalisation  : 

fA  = 
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Pour  aller  plus  loin  on  admet  que  les  deux  echelles  integrales -A_  et  A  sont  voisines.  C'est  le  cas 
lorsque  l'on  admet  pour  la  fonction  de  correlation  R  une  loi  exponent iel le  ;  on  a  alors  rigoureusement 
A  *  .A-  •  Pour  une  gauss ienne  on  trouve  : 


A  = 


A  -  o,  ?  A 


Or  ce  sont  deux  cas  extremes  correspondant  ^  un  nombre  de  Reynolds  de  turbulence  soit  tr£s  grand  (cas 
de  1 ' exponentielle) ,  soit  tr£s  faible  (cas  de  la  gaussienne)  [9]  •  On  admettra  que  pour  les  Reynolds  inter¬ 
mediates  cette  quasi  £galite  tient  encore.  Dans  ces  conditions  la  correlation  normalisee  en  faisceaux 
croises  jointe  S  la  connaissance  des  longueurs  d ' integration  2  XA  et  2  Ji ^  permet  de  calculer  l'€chelle 
integrale  :  n  / - 

x.AxJuu, 

Il  est  A  remarquer  que  la  methode  des  faisceaux  orthogonaux  ne  faisait  pas  intervenir  les  longueurs 
d' integration.  Les  deux  methodes  sont  done  independantes .  Dans  les  deux  cas,  les  correlations  sont 
d'autant  plus  faibles,  done  difficiles  3  mesurer,  que  ces  longueurs  sont  grandes  devant  1 'echelle  inte¬ 
grale. 


3.1.3  -  InjJfuence  dc_l[£p(L<AA£uA_deA  fiai&ceaux 

Jusqu'ici  on  a  suppose  les  faisceaux  infiniment  etroits,  ce  qui  n'est  pas  physiquement  realisable. 

En  ce  qui  concerne  le  decouplage  des  fluctuations  de  temperature  et  des  fluctuations  d ' emissivite ,  leur 
epaisseur  finie  n'en  restreint  pas  la  qualite  puisque  dans  le  cadre  de  1 ' approximation  geometrique  un 
faisceau  optique  est  la  superposition  de  rayons.  Par  contre  les  calculs  de  la  variance  et  des  correlations 
peuvent  en  etre  affectes  puisque  l’on  ajoute  une  integration  supplcmentaire  sur  les  dimensions  transver- 
sales  des  faisceaux. 
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L' importance  de  cet  effet  a  ete  evaluee  dans  le  cas  de  faisceaux  parall^lepipediques  de  section 
l^€.^  pour  celui  dirige  selon  l'axe  1 ,  et  Lf  pour  celui  dirige  selon  l'axe  2. 


On  obtient  : 


a_ 

x 


L ,  •/}  'a 

si'1  =  a'!  =  *  r\A. 


^  •*  — t 

oQ  A  est  une  troisieme  echelle  integrale  definie  par  : 

rL  -  ^  .  r  ±-  UL  d  & 

J-  "T1-  4 

=  * 


C  A 


Dans  le  cas  oQ  € 
croises  : 


%  W 

$  =  ^  ,  on  en  deduit  la  correlation  normalisee  en  faisceaux 

£  r  i 


P(»)-  ^  T T  -A. 

^7^  s/umn 


i  v-  aX 


A 


-) 


Si  l'on  admet  que est  du  meme  ordre  de  grandeur  que  -A.  et  -A-  ,  les  termes  correctifs  sont 

i'ordre  ( t/jJ)  ,  Cela  signifie  que  l'epaisseur  des  faisceaux  doit  etre  petite  devant  l'echelle  integrale. 


comme  on  pouvait  s'y  attendre. 


Pour  des  faisceaux  larges  de  2  mm  et  une  echelle  integrale  d'environ  1  cm,  la  precision  des  mesures 
est  encore  de  1%. 

3.1.4  -  Application  decette methode  a _VJ.tu.de.  do.  ia_tuA.bLii.cnce.  theAmique 

La  technique  de  mesure  des  fluctuations  de  temperature  et  les  methodes  des  faisceaux  orthogonaux  et 
des  faisceaux  croises  ont  ete  utilises  sur  le  jet  3  la  sortie  d'un  foyer  d'etude  de  la  dilution  malgre 
V  inconvenient  de  1 'heterogeneite  des  temperatures  moyennes.  On  procede  en  deux  temps  :  tout  d'abord 
mesure  des  fluctuations  de  temperature  integrees  en  emission  et  en  absorption  simultanees,  ensuite  deter¬ 
mination  de  l'echelle  integrale.  En  effet,  si  l'on  admet  que  les  caracterist iques  de  la  turbulence  sont 
independantes  de  la  grandeur  mesuree  (pourvu  que  ce  soit  un  scalaire)  les  techniques  des  faisceaux  ortho¬ 
gonaux  ou  croises  peuvent  etre  utilisees  aussi  bien  en  emission,  ce  qui  presente  l'avantage  de  simplifier 
1 ' appareillage  et  sa  mise  en  oeuvre  en  supprimant  les  emetteurs  des  deux  appareils  necessaires.  Cela 
permet  egalement  d'augmenter  considerablement  le  rapport  signal  sur  bruit  puisqu'il  n'est  plus  necessaire 
d'etre  monochroma tique . 

a)  Me&uAe_deb  _&luctuatiqnb  de_tempeAatuAe_integA^eb 

La  mesure  des  fluctuations  de  temperature  a  ete  utilisee  a  la  suite  de  la  mesure  de  temperature 
moyenne  en  ajustant  la  temperature  de  la  source  &  celle-ci,  il  est  necessaire  auparavant  de  remplacer  le 
filtre  interferentiel  S  4,5  jam  par  un  filtre  &  2,5 jjm  pour  lequel  le  facteur  d'emission  est  inferieur  S 
5%.  Le  principe  de  la  mesurfe  ne  necessite  plus  qu'une  seule  voie  de  mesure  et  la  voie  occultee  par  le 
couteau  peut  etre  utilisee  pour  la  determination  du  facteur  d'emission  moyen,  les  pointes  de  temperature 
pouvant  aller  jusqu'&  100  K  au-dessus  de  la  temperature  moyenne.  La  valeur  quadratique  moyenne  des 
fluctuations  integrees  sur  le  chemin  optique  est  egale  &  : 

/P  =  i  5  K 

b)  Technique  deb  fiaUceaux  orthogonaux 

On  utilise  done  deux  pyrometres  qui  int£grent  1' emission  infrarouge  du  jet  dans  tout  le  domaine 
spectral  jusqu'S  4,25/jym.  L’un  d'eux,  maintenu  fixe,  est  place  vert icalement  au-dessus  de  la  veine  ; 

1' autre  est  horizontal  et  peut  etre  translate  axialement  grace  S  un  chariot  mobile  commande  &  distance. 


La  figure  10  presente  les  courbes  de  correlation  spat iotemporelles  obtenues  &  partir  des  deux  mesures. 
Chaque  courbe  correspond  h  une  position  relative  des  deux  faisceaux,  &  etant  le  parametre  de  correlation 
temporelle  et  la  plus  courte  distance  entre  les  deux  faisceaux.  cL  est  modifie  par  pas  de  10  mm  sauf  & 
l'approche  de  1 ' intersection  oO  le  pas  est  resserre  £  2  mm.  Pour  ne  pas  surcharger  la  figure,  les  courbes 
correspondantes  n'y  figurent  pas.  A  une  legere  attenuation  pres  les  courbes  se  deduisent  les  unes  des 
autres  par  une  translation.  Ceci  s’explique  si  l'on  admet  1'hypothSse  de  Taylor  ;  ef fectivement ,  si  on 
porte  sur  un  diagramme  1’abscisse  'fc  / max  relative  au  maximum  des  courbes  de  correlation,  en  fonction  de 
la  distance  on  obtient  une  droite  et  la  vitesse  definie  par  la  pente  de  cette  droite  (84  m/s)  est  proche 
de  la  vitesse  attendue  pour  le  jet. 


Les  points  3  dephasage  temporel  nul  ( *  0)  permettent  de  tracer  la  courbe  P )  qui  presente  un 
maximum  pour  A  =  0  (Fig.  11).  D'apres  Fischer  et  Krause  si  on  normalise  cette  fonction  par  la  valeur 
P  correspondant  au  maximum,  on  obtient  la  fonction  de  correlation  RCdt)  entre  les  points  A  et  B  . 
Son  integration  permet  de  calculer  une  echelle  integrale.  On  trouve  2,5  cm.  En  fait,  1 ' interpretation  de 
Fischer  et  Krause  n'est  pas  correcte  car  elle  suppose  negligeables  les  correlations  entre  les  points 
voisins  de  A  et  B  ,  et  jin  a  vu  qu'il  est  necessaire  pour  obtenir  ft  id)  de  tracer  la  correlation  en 
fonction  du  parametre  <4.  au  lieu  de  oL  et  d'en  prendre  la  derivee.  Au  voisinage  de  l'origine  le  calcul 
de  cette  derivee  est  deiicat  ;  pour  l'effectuer  on  represente  la  fbnetion  par  un  developpement  de  Lagrange 
du  7eme  degre.  Loin  de  l'origine  la  derivee  est  determinee  graphiquement .  On  obtient  ainsi  la  courbe  de  la 
figure  12  qui  represente,  3  un  coefficient  pr£s,  la  fonction  de  correlation  cherchee.  On  voit  quj^elle  est 
beaucoup  plus  etroite  que  P(fll)  .  La  valeur  &  l’origine  permet  de  calculer  l'echelle  integrale  et 

1 ' integration  sur  A  permet  de  calculer  l'echelle  integrale  A.  .  On  les  trouve  tr£s  proches  l'une  de 
1' autre  :  ^ 

A  -  0,74  cm  et  A  -  0,72  cm 


II 


ce  qui  est  une  justification  supplementaire  de  1 ' approximation  qui  consiste  &  les  supposer  egales,  II  faut 
cependant  noter  que  ces  mesures  ne  sont  qu' approximat ives ,  En  particulier,  le  calcul  par  integration  est 
arrEtE  arbitrairement  au  premier  zero  de  la  courbe  de  correlation.  On  suppose  que  la  contribution  du  reste 
de  la  courbe,  alternativement  positive  et  negative,  est  negligeable. 

c)  Co  violation  en^at*ceaux_croi*§* 

Une  autre  methode  pour  atteindre  .A.  est  d'utiliser  la  correlation  normalisee  en  faisceaux  croises 
donnee  par  :  _ 

f<»).  A-T  ^  . », II, 

Cette  determination  fait  intervenir  les  longueurs  d' integration  2  ^  et  i  ^  l  .  Dans  le  cas  present 
l  X*  »  5  cm  et  tAl  -  30  cm.  On  en  deduit  une  echelle  integrate  JL  =  0,94  cm  sensiblement  plus  grande 
que  la  precedente  determination  ;  rnais  il  faut  noter  que  du  fait  de  1' existence  d'un  profil  de  temperature 
moyenne,  la  hauteur  du  jet  2  A  a  ne  coincide  pas  avec  l'Epaisseur  de  la  zone  d'emission  et  conduit  &  une 
surestimation  de  l'Echelle  integrale.  Pour  egaler  les  deux  determinations  il  faudrait  -  3  cm. 

L'ordre  de  grandeur  de  ces  mesures  est  trois  &  quatre  fois  plus  faible  que  celle  selon  1 ' interpretation  de 
Fischer  et  Krause. 

d)  Fluctuation*  de  temperature  locale 

Connaissant  l'echelle  integrale  il  est  maintenant  possible  de  determiner  la  valeur  quadratique 
moyenne  des  fluctuations  de  temperature  locale,  en  utilisant  la  relation  : 


—  % 

(-%>) 


i  A 

TZ7 


Pour  effectuer  le  calcul  on  retient  pour  l'Echelle  integrale  la  valeur  moyenne  entre  0,94  cm  et 
0,72  cm  soit  0,83  cm.  Le  facteur  multiplicatif  est  alors  egal  &  1,74  et  1* intensity  des  fluctuations  passe 
de  25  K  &  43  K.  L'intensite  de  la  turbulence  thermique  est  d'environ  7%. 

e)  Spectre 

En  principe  on  peut  obtenir  le  spectre  &  partir  de  la  fonction  de  correlation  par  transformation  de 
Fourier.  Mais  cette  fonction  de  correlation  a  ete  obtenue  point  par  point  et  sa  definition  est  insuffi- 
sante  notamment  pour  les  points  proches  de  l’origine  qui  fournissent  le  domaine  hautes  frequences  du 
spectre.  Aussi  s'est-on  contente  de  1' analyse  spectrale  directe  des  mesures,  les  spectres  obtenus  Etant 
alors  modifies  par  l'effet  d' integration.  Par  exemple,  les  deux  spectres  correspondant  &  la  mesure  en 
faisceaux  croises,  traces  en  coordonnees  semi-logarithmiques  jusqu'S  10  kHz  ont  une  allure  non  classique 
sensiblement  exponent ielle ,  de  frequence  caracteristique  environ  1300  Hz. 

rf)  Flamme _premelangee_turbulente 

Cette  methode  a  ete  appliquee  pour  determiner  la  fonction  densite  de  distribution  dans  le  cas  d'une 
flamme  turbulente  premelangee.  Une  telle  flamme  est  representee  sur  la  figure  13,  un  Ecoulement  de  gaz 
brules  chauds  est  parallEle  5  un  ecoulement  de  gaz  frais,  une  flamme  bidimensionnelle  se  developpe, 
l'abscisse  ce.  est  dans  le  sens  de  1 ' ecoulement  et  l’analyse  des  fluctuations  de  temperature  se  fait  3  une 
abscisse  ce  determinee  et  &  plusieurs  ordonnees  ^  situees  ^  la  traversee  du  front  de  flamme.  On  observe 
comme  le  montre  la  figure  13  deux  distributions  pour  deux  temperatures  semblant  indiquer  des  paquets  de 
gaz  frais  et  de  gaz  brflles. 

Nous  avons  insiste  sur  les  techniques  de  mesures  de  la  temperature  par  des  methodes  optiques  car 
elles  presentent  un  certain  nombre  d'avantages  mais  cette  mesure  ne  va  pas,  dans  1 ' interpretation,  sans  un 
certain  nombre  de  difficultes. 

3.2  -  MZthode*  par  *onde* 

Il  existe  de  tres  nombreuses  methodes  de  mesure  des  fluctuations  de  temperature  par  thermocouple, 
resistance  ou  thermistance,  comme  precedemment  nous  en  d^crivons  une  qui  montre  les  avantages  et  les 
inconvenients  de  cette  technique.  C'est  celle  developpe-0  par  Yule,  Taylor  et  Chigier  [9]  ,  il  s'agit  de 
mesure  de  la  temperature  par  thermocouple,  le  probldme  la  plus  important  dans  le  cas  d'une  mesure  de  la 
temperature  dans  un  Ecoulement  turbulent  etant  le  calcul  et  la  mesure  du  temps  de  reponse  du  thermocouple 
permettant  d'acceder  a  la  temperature  du  gaz  T  ,  &  partir  de  la  mesure  donnee  par  le  thermocouple  T4, 
et  la  constante  de  temps  f  ’ 

wr£ 

Le  schema  d'un  tel  thermocouple  est  representE  sur  la  figure  14  qui  montre  de  quelle  maniEre  on  peut 
concevoir  le  support  pour  ne  pas  trop  perturber  l'Ecoulement ,  La  temperature  de  la  jonction  ou  de  la  perle 
TJu  est  reliEe  £  la  tempErature  du  gaz  T<^.  en  considErant  trois  mEcanismes  :  le  transfer!  gaz-jonction, 
le  transfert  de  chaleur  &  travers  le  fil  rEsultant  d'une  dEperdition  par  conduction  et  le  transfert  de 
chaleur  entre  le  fil  et  le  gaz.  La  conduction  thermique  dans  le  gaz  peut  Etre  considerEe  comme  negligeable. 
Une  Etude  de  bilan  d'Energie  est  une  bonne  approximation  pour  avoir  une  idee  de  la  valeur  de  la  constante 
de  temps  et  des  parametres  qui  agissent  sur  celle-ci. 

La  disposition  de  ce  thermocouple  conduit  £  admettre  une  diffErence  de  constante  de  temps  pour  le  fil 
seul  et  pour  la  perle  placEe  autour  de  la  jonction.  En  particulier  le  temps  de  rEponse  dEpend  fortement  de 
la  dimension  de  la  perle  par  rapport  au  fil  constituant  le  thermocouple.  La  rEponse  du  thermocouple  est  en 
genEral  calculee  ou  determinEe  expErimentalement  soit  en  rEalisant  dans  le  gaz  un  Echelon  de  tempErature, 
soit  en  modulant  £  une  frEquence  constante  la  tempErature  du  gaz.  Cette  fa^on  de  faire  est  gEnErale  pour 
la  dEtermination  d'une  constante  de  temps  :  modulation  sinusoidale,  Echelon,  ou  technique  impulsionnel le . 

On  peut  en  effet  produire  un  Echelon  ou  une  impulsion  de  tempErature  dans  le  thermocouple  en  appor- 
tant  de  la  chaleur  dans  le  fil  par  un  Echelon  ou  une  impulsion  de  courant  Electrique. 
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Comme  nous  l'avons 
du  fil  plus  celle  de  la 
celle  de  la  perle.  Dans 
est  de  sorte  que  : 


signale  cette  constante  de  temps  correspond  en  fait  3  la  constante  de  temps  £ 
perle  Vm  et  dans  son  evaluation  tout  depend  du  rapport  de  la  dimension  du  fil  a 
le  cas  d'un  fil  support  de  longueur  infinie  la  constante  de  temps  qui  intervient 


V  VW 


dans  le  cas  d'un  thermocouple  de  longueur  finie  on  utilise  la  relation 


T.  i 


T 


T. 


Dr« 


t  7>t 


£.  _  A 


(7. 


oQ  intervient  un  coefficient  d'echange  &  entre  la_ temperature  du  fil  au  loin  et  la  temperature 

de  la  perle  et  une  constante  de  temps  de  la  perle  i .  La  methode  la  plus  exacte  consiste  3^mesurer 
par  une  methode  de  chauffage  electrique  impulsionnel  et  au  raoyen  d’un  calculateur  d'evaluer  et  fa  et 
done  d'obtenir  les  elements  correctifs  necessaires  pour  calculer  la  temperature  instantanee  du  gaz  a 
partir  de  la  temperature  mesuree  3  la  jonction  X  . 

T 

Dans  la  conception  de  cet  appareil  nous  avons  deux  types  de  calculateurs  : 


a)  l'un,  qui  admet  comme  donnees  les  parametres  geometriques  du  thermocouple  (fil  et  perle)  et  les  para- 
metres  mesur£s  pendant  1 'utilisation,  il  est  possible  de  cette  maniere  de  determiner  la  constante  de 
temps  de  la  perle  dans  les  meilleures  conditions,  c'est-3-dire  pratiquement  les  conditions  de  1 'expe¬ 
rience  et  de  calculer  ainsi  les  tenues  correctifs  dus  aux  pertes  thermiques  ; 

b)  un  autre  calculateur  qui,  a  partir  de  la  constante  de  temps  et  de  la  temperature  de  la  fonction,  deter¬ 
mine  la  temperature  instantanee  du  gaz. 


Les  sources  d'erreurs  avec  cette  technique  proviennent  de  plusieurs  causes  : 
aX  tout  d'abord  de  1 ' imprecision  dans  la  mesure  de  la  constante, 

b)  de  la  connaissance  plus  ou  moins  variable  pour  chaque  thermocouple  des  geometries  qui  different  d'un 
thermocouple  3  un  autre  et  meme  evoluent  en  cours  d'essai, 

c)  de  la  raeconnaissunce  des  perturbations  apportees  par  le  thermocouple  ou  la  sonde  3  1 ' €coulement , 
quoique,  3  l'heure  actuelle,  ces  inconvenients  aient  ete  reduits  par  l’emploi  dans  les  flammes  de  fils 
de  diametre  inferieur  3  5  pm  mais  qui  cependant  doivent  tenir  compte  de  1' impact  des  particules, 

d)  du  manque  de  donnees  sur  les  transferts  de  chaleur  en  transitoire,  cet  aspect  est  en  effet  important  et 
pratiquement  inconnu,  d’oQ  le  manque  de  precision  qui  en  decoule. 

Le  fait  d'utiliser  deux  thermocouples  voisins,  l'un  chauffe  electriquement  par  un  courant  alternatif 
&  haute  frequence  et  1 'autre  pas,  elimine  certaines  causes  d'erreurs  ;  en  effet  la  difference  de  la  force 
electromotrice  des  deux  thermocouples  permettrait  d'atteindre  les  coefficients  de  transfert  inetantanes  et 
qui  peuvent  etre  utilises  par  la  suite  pour  determiner  la  compensation  necessaire  3  appliquer  au  thermo¬ 
couple  non  chauffe. 


Par  cette  methode  de  la  sonde  on  determine  en  principe  les  fluctuations  de  temperature  en  un  point  et 
l'on  peut  avoir  ainsi  une  statistique  eulerienne. 

3.3  -  Vu>ciLt>t>iovi  6uA  £&&  methods  de  m&6u/ie  d(u  fluctuation!,  da  Impo-XatuAc 

Nous  avons  insist^  sur  la  mesure  de  ce  scalaire  pour  bien  montrer  les  limites  de  ces  techniques  et 
les  progres  qui  restent  3  faire. 

Les  methodes  de  mesure  de  la  temperature  par  voie  optique  sont  tres  interessantes  et  d'une  mise  en 
oeuvre  relativement  facile,  on  obtient  un  signal,  mais  il  faut  1 ' interpreter  correctement ,  ce  qui  n'est 
pas  toujours  aise,  nous  n' avons  decrit  que  la  methode  du  renversement  de  raies  dans  1' infrarouge  car  elle 
nous  parait  d'une  relative  facilite  mais  il  en  existe  d'autres  basees  par  exemple  : 

a)  sur  la  mesure  absolue  de  1' intensity  de  raies, 

b)  sur  le  rapport  de  1'intensite  de  deux  raies  (emission  de  raies  due  3  des  metaux  contenus  dans  la 
f lamme) , 

ou  encore, 

c)  sur  le  rapport  d'intensite  de  bandes  qui  sont  li£es,  dans  le  cas  de  molecule  diatomique,  au  mouvement 
interne  de  la  molecule  ;  en  spectrographie  Raman  par  exemple,  la  mesure  d'intensite  de  raies  Raman  ou 
du  rapport  d'intensites  des  raies  Stokes  et  Anti-Stokes  sera  etudiee  au  paragraphe  suivant  lorsque  la 
mesure  des  concentrations  sera  abordee. 

Les  methodes  de  mesure  3  court  temps  de  reponse  de  la  temperature  par  sonde  offrent,  comme  nous 
l'avons  souligne,  de  nombreuses  difficultes.  Nous  avons  choisi  de  d^crire  la  methode  du  thermocouple  parce 
qu'elle  est  ponctuelle.  Les  methodes  3  variation  de  resistance  sont  en  general  coupl§es  avec  les  techni¬ 
ques  du  fil  chaud  et  elles  servent  3  calculer  le  terme  de  compensation  necessaire  pour  une  determination 
correcte  de  la  vitesse,  elles  pr£sentent  1 ' inconvenient  de  ne  pas  etre  ponctuelles  et  d'int£grer  les 
phenomenes  de  transfert  instat ionnaires  sur une  certaine  longueur  de  la  sonde.  Ces  sondes  couplees  sont  tres 
utiles  car  elles  permettent  de  determiner  des  correlations  du  type  T/^/  -axTT'  f  T'y  ^  ''J> 

qui,  dans  un  probleme  bidimensionnel ,  determinent  outre  les  correlations  des  fluctuations  du  scalaire  la 
diffusion  de  ce  scalaire  dans  deux  directions  orthogonales  correspondant  3  des  fluctuations  de  vitesses  de 
composantes  ax'  et  -v'  .  La  figure  15  est  un  exemple  de  realisation  due  3  H.  Fiedler  et  P.  Mensing  F 10l  qui 
couple  1 'effet  de  la  vitesse  et  celui  de  la  temperature. 
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4  -  METHOPES  PE  MESURE  INST ANT ANEE  VES  CONCENTRATIONS  ET  VE  LA  MASSE  VOLUMIQUE 

Les  techniques  de  diagnostic  portant  sur  la  mesure  des  concentrations  ont  et£  consid€rablement  d£ve- 
loppees  ces  dernieres  annees,  quoique  prometteuses ,  il  faut  cependant  distinguer  dans  ces  m€thodes,  et  ce 
n'est  pas  tou jours  evident  S  la  lecture  des  rapports,  cel les  qui  donnent  des  concentrations  moyennes  et 
detectent  des  especes  minoritaires  comme  les  radicaux  (Of/  )»  celles  qui  conduisent  A  la  mesure  des 
concentrations  instantanees  et  qui  interessent  done  directement  les  ecoulements  turbulents  ainsi  que 
celles,  toujours  instantanees,  mais  qui  ne  distinguent  pas  la  nature  des  molecules  mais  sont  plutot  sensi- 
bles  aux  nombres  de  molecules  par  unite  de  volume,  c'est-A-dire  aux  fluctuations  de  la  masse  volumique. 
Toutes  ces  methodes  sont  optiques  et  utilisent  des  techniques  lasers. 

Nous  distinguerons  dans  ce  rapport  trois  methodes  qui  coramencent  S  donner  quelques  resultats  intlres- 
sants  : 

-  la  methode  de  la  diffusion  Raman  spontanee  ; 

-  la  methode  utilisant  la  fluorescence  ; 

-  la  methode  de  diffusion  Raman  antistokes  coherente  (DRASC)  ou  (CARS)  coherent  antistokes  Raman  scatte¬ 
ring  que  nous  allons  brievement  decrire  en  donnant  les  performances  actuelles,  cette  dernifcre  determi¬ 
nant  surtout  la  masse  volumique. 


4.1  -  M$thode  de  la  cU^iuion  Raman  Apontanea  IVRS)  [11] 

Lorsqu'un  faisceau  lumineux  traverse  un  gaz,  une  faible  fraction  du  flux  incident  est  diffuse  avec  un 
changement  de  frequence  dans  7T  steradians.  C'est  la  Diffusion  Raman  Spontanee  dans  laquelle  les  photons 
incidents  echangent  de  l'energie  par  valeurs  discretes  avec  le  systeme  moleculaire.  Le  spectre  du  flux 
diffuse  se  compose  de  raies  dont  les  positions  sont  representatives  des  divers  niveaux  energetiques  du 
systeme  moleculaire.  Generalement  on  s’interesse  aux  molecules  diatomiques  et  aux  transitions  de  vibration- 
rotation.  Ce  type  de  diffusion  est  bien  adapte  a  1' analyse  des  melanges  de  gaz  car  les  spectres  de  DRS 
vibrationnelle  n' interferent  que  rarement.  Examinons  le  cas  d'etats  vibration-rotation  de  la  molecule 
caracteriseeopar  les  nombres  quantiques  v  et  0  ,  un  processus  DRS  fait  passer  de  l'etat  (v  ,  3)  &  l'etat 

(  V T  *J  /  *3r  i  ),  la  difference  de  frequence  entre  le  flux  excitateur  et  le  flux  diffuse  etant  repre¬ 
sentative  de  cette  transition.  La  diffusion  Stokes  correspond,  pour  la  molecule,  b  un  gain  d'energie  et 
celle  Anti-Stokes  &  une  perte.  Soient  le  nombre  de  molecules  par  unite  de  volume,  Tv  et  Tr  les  tempe¬ 
ratures  de  vibration  et  de  rotation,  R.  la  puissance  incidente,  ( ^  \  ,  l  ^  .  les  sections 

efficaces  dif ferentielles  de  diffusion  pour  la  branche  Stokes  et  la  branche  Anti-Stokes,  Q0  la  constante 
rotationnelle ,  ->2_  1* angle  solide  de  collection,  la  frequence  de  vibration, ^ la  puissance  pour  une 
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Pour  realiser  la  DRS  on  utilise  un  montage  represente  par  le  schema  de  la  figure  16  comprenant  une 


^ROT  -  X!  (•?  3  + 


source  excitatrice  et  1  Ti / £  1 'analyse  par  un  spectrometre  de  la  lumiere  diffusee. 

D'apres  les  relations  precedentes  donnant  et  on  remarque  que  le  spectre  de  DRS  evolue  en 

fonction  de  ^  ,  de  TR  et  de  Tv  mais  on  peut  supposer  l'equilibre  thermodynamique  realise  et  suppose 
t  Tv  z  "7”  •  En  faisant  le  rapport  d'intensite  de  deux  raies  on  obtient  la  temperature,  on  peut,  par 

exemple,  en  prenant  V  =  I  faire  le  rapport  du  flux  Stokes  emis  dans  la  branche  G2.  fondamentale  (due  aux 
molecules  initialement  dans  l'etat  vibrationnel  fondamental)  au  flux  Stokes  emis  par  la  branche  CQ  due 
aux  molecules  qui  sont  initialement  dans  un  etat  vibrationnel  excite  donne,  c»  qui  correspond  A  la  rela¬ 
tion  :  -^cev 

At 

0,3)  — *(«3J]  i 

On  peut  egalement  effectuer  le  rapport  du  flux  integre  sur  l' ensemble  de  la  branche  Q  Stokes 
celui  integre  sur  1 'ensemble  de  la  branche  int£gree  Anti-Stokes,  ce  rapport  correspond  &  la  relation  ; 

ps  Lev, d — .(vtoq 

La  selection  spectrale  peut  etre  realisee  dans  les  deux  cas  par  filtres  interferent iels. 


£t  Ps  [ip,  3) -^>(-1,3)1  ^  ^ 
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La  temperature  etant  connue,  1' intensity,  integrie  dans  une  bande  spectrale  de  largeur  et  de  position 
quelconques,  est  proportionnelle  71  la  concentration  dans  le  gaz  de  la  molecule  choisie.  Suivant  que  la 
temperature  est  inferieure  ou  superieure  71  1500  K  on  prendra  soit  la  structure  rotationnelle ,  soit  la 
structure  vibrationnelle .  Done  1' intensite  integree  J  dans  une  bande  de  largeur  spectrale  A  A  est  egale 

I  =  cxitf  FOX.T) 

En  prenant  une  longueur  ^  *  1  mm,  =  0,6§  steradian,  pour  1' azote  &  ■  2  10^  molecules/cm^  et 
(  <k<r  /  <ISI  ^  =  5,4  lo3  (  raie  d'excitation  4880  A)  on  collecte  7  10“13,  soit  2  106  photons  pour  un  joule 
incident.  Mais  avec  une  optique  de  collection  &  10%,  un  detecteur  de  10%,  on  detecte  2  104  photoelectrons 
avec  une  erreur  relative  d' environ  1%. 

Les  performances  escomptees  sont  difficiles  71  atteindre  par  suite  des  flux  parasites  :  luminosity  du 
milieu,  fluorescence,  rayonnement  des  particules . . . 

La  source  excitatrice  peut  etre  continue  ou  pulsee,  malheureusement  il  n'est  pas  possible  d'obtenir 
les  fluctuations  avec  une  source  A  excitation  continue  (mesure  longue  pour  collecter  suf f isamment  de 
photons)  mais  avec  une  excitation  par  laser  pulse  il  est  possible  d'etudier  les  fluctuations  de  tempera¬ 
ture  et  de  concentration. 

Un  tel  montage  est  represente  sur  la  figure  17.  Trois  conditions  doivent  etre  r^alisles  : 

-  la  duree  de  1* impulsion  laser  doit  etre  courte  devant  les  temps  caracteristiques  de  la  turbulence, 

-  l'energie  doit  etre  suffisante  pour  permettre  une  mesure  precise, 

-  les  mesures  necessaires  71  la  determination  de  et  de  "T"  doivent  etre  effectules  simultaneraent . 

Actuellement  avec  un  laser  3  rubis  relaxe  de  30  Joules  et  pour  une  impulsion  laser  de  700  jjs,  le 
temps  de  mesure  est  de  l'ordre  de  10  *is,  la  precision  de  la  mesure  sur  la  concentration  de  .Ng,  mesuree 
dans  des  conditions  normales  par  la  mesure  de  la  raie  Stokes  est  de  l'ordre  de  2%  et  de  4%  71  1500  K.  Le 
volume  de  mesure  correspond  71  un  cylindre  de  1  mm  de  diametre  et  >(  *  1  mm  pour  la  longueur. 

L'extension  possible  de  cette  methode  consiste  71  realiser  des  impulsions  plus  longues  de  10  71  20  ms 
avec  des  temps  de  mesure  de  10  ps,  mais  cette  extension  se  fait  au  detriment  de  la  precision. 

4.2  -  Methods  utita>ayvt  ta  filuo* iz&cence.  [72]  [73] 

Cette  methode  est  analogue  &  la  DRS  mais  au  lieu  de  considerer  la  fluorescence  comme  un  effet  para¬ 
site  qui  vient  perturber  la  mesure  en  DRS,  on  la  considere  au  contraire  comme  le  phynomene  principal  et 
elle  sert  de  traceur  pour  les  autres  molecules  entrant  dans  les  produits  de  combustion.  Taran,  Pealat  et 
Baillyfll]  ont  montre  que  cette  technique  etait  valable  pour  determiner  la  temperature  et  la  concentration 
de  certaines  molecules  donnant  le  phenomene  de  fluorescence  .  La  technique  experimental  est  voisine  de  la 
DRS.  Le  spectrographe  peut  etre  remplace  par  un  filtre  qui  elimine  la  composante  Rayleigh  et  Mie.  Le 
photomult iplicateur  collecte  alors  toute  la  fluorescence  ainsi  que  les  raies  Raman  individuelles  mais  dont 
la  contribution  au  signal  total  est  negligeable.  Il  faut  cependant  que  les  gaz  contiennent  des  molecules 
f luorescentes ,  ceci  a  lieu  dans  les  zones  de  pyrolyse  des  hydrocarbures  en  particulier  dans  les  flammes 
air-kerosene  et  dans  le  domaine  de  temperature  compris  entre  600  et  1000  K,  au-delS  de  1000  K  ce  phenomene 
de  fluorescence  disparait  car  les  molecules  f luorescentes  sont  detruites.  L'avantage  en  faveur  de  la 
fluorescence  provient  de  son  signal  important,  ce  qui  facilite  les  mesures,  faut-il  encore  pouvoir  inter¬ 
preter  ces  mesures,  car  les  fluctuations  de  concentration  peuvent  etre  dues  71  plusieurs  causes  qui  inter- 
viennent  dans  le  bilan  de  l'espece  :  d'abord  un  effet  physique  du  au  transfert  de  masse  par  diffusion 
moleculaire  et  turbulent  et  ensuite  un  effet  chimique  de  production  et  de  destruction  difficile  71  evaluer 
en  milieu  turbulent.  Quoiqu'il  en  soit  les  performances  de  cette  methode  sont  actuel lement ,  avec  un  laser 
71  Argon  ionise  de  10  W  continu,  de  50  ps  en  ce  qui  concerne  le  temps  de  reponse,  c '  est-71-dire  que  l'on 
peut  mesurer  en  continu  des  frequences  allant  de  100  Hz  71  20  kHz,  la  resolution  spatiale  correspond  71  un 
cylindre  de  1  mm  de  diametre  et  1  cm  de  longueur.  Cette  methode  ne  marche  qu'entre  600  et  1000  K  et  le  gaz 
doit  contenir  des  hydrocarbures  qui  se  pyrolysent  71  cette  temperature.  Il  existe  dans  la  litterature  de 
nombreuses  publications  sur  1 ' utilisation  du  phenomene  de  fluorescence,  en  particulier  sur  la  mesure  des 
concentrations  de  radicaux  tels  que  CN,  CH,  OH  utilisant  des  lasers  pour  provoquer  une  fluorescence  satu- 
ree,  1 'etude  se  faisant  par  analyse  de  1 'emission  ou  de  1 'absorption,  mais  ces  methodes  ne  sont  pas  instan- 
tanees  ;  Bonczyk  et  Shirley  [16]  pretendent  que  leur  methode  est  instantanee  mais  ils  n'ont  publie  aucun 
resultat  en  ce  sens.  La  litterature  est  abondante  mais  ne  repond  pas  aux  problemes  poses  par  1 'etude  de  la 
turbulence . 

4.3  -  Methode  du  di^mlon  Raman  anti- Stoker  coheAente.  (PRASC)  [J4l  [75] 

Il  semble  que  1 ' utilisation  de  la  DRS  est  delicate  pour  ne  pas  dire  impossible  dans  le  cas  de  flammes 
instables  chargees  de  particules  brulees  ou  imbrGlees  qui  ont  dans  le  champ  une  trop  grande  luminosite  ; 
e'est  pourquoi  la  diffusion  Raman  anti-Stokes  coherente  (DRASC)  prend  de  1' importance  en  particulier  le 
rejet  de  la  lumiere  parasite  est  beaucoup  plus  facile  qu'en  DRS. 

Le  principe  est  le  suivant.  On  observe  un  signal  DRAS  lorsque  deux  ondes  planes  et  colineaires  de 
frequence  ^  p  et  sont  envoyees  sur  des  molecules  ayant  des  vibrations  actives  caj*  telles  que 

.  Une  onde  7i  la  frequence  anti-Stokes  -  2  est  diffus^e  dans  la 

di-ection  des  ondes  excitatrices  et  son  intensite  est  proportionnelle  au  carre  de  la  concentration  des 
mo  ecules  diffusantes.  Lorsqu'on  balaie  la  frequence  par  rapport  7*  le  signal  anti-Stokes  emis 

pr»  sente  done  une  augmentation  resonante  chaque  fois  que  r  c‘L>'fa  “ cu 5  est  egal  &  la  frequence  d'une 

transition  Raman  du  milieu.  Le  faisceau  diffuse  est  lui-meme  coherent  et  pr£sente  toutes  les  caract£ris- 
tiques  d'un  faisceau  laser  :  grande  intensite,  petit  angle  solide,  et  dans  la  meme  direction  que  les  ondes 
excitatrices,  ceci  rend  1 ' eliminat ion  de  la  lumiere  parasite  beaucoup  plus  aisee  nu'enDRS.La  fluorescence 
est  eliminee  d'autant  plus  facilement  que  le  signal  se  trouve  emis  du  c6t£  anti-Stokes  de  la  frequence 
excitatrice  alors  que  la  fluorescence  est  emise  du  cote  Stokes. 
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La  puissance  anti-Sotkes  collectee  est  proportionnelle  au  carre  d'une  susceptibilite  "X  du 

milieu  qui  traduit  la  reponse  des  molicules  3  l'excitation  coherente  ;  JC  est  la  somme  de  deux  contribu¬ 
tions  et  J<- :  pc  r  est  la  partie  resonante  et  represente  la  contribution  des  vibrations  actives 

en  Raman  et  de  frequence  proches  de  ;  pour  la  transition  Raman  £  y,  0  )  _ 

XRoc  [n(v,3>  - 

done  proportionnelle  au  nombre  de  molecules  par  unite  de  volume  de  l'espece  moleculaire  qui  poss£de  la 
transition  Raman  analysee  ;  X.  ^  ^  represente  la  tres  faible  contribution  des  autres  mouvements  de 
vibration  ainsi  que  de  la  non  linearite  de  distorsion  du  nuage  electronique  des  molecules,  e'est  un  para- 
metre  proportionnel  au  nombre  de  molecules  par  unite  de  volume  de  toutes  les  especes  moleculaires  presen- 
tes  au  voisinage  de  la  mesure.  En  faisant  varier  ^  -  ooc,  au  voisinage  des  frequences  de  vibration 
on  obtient  un  spectre  representat if  de  la  variation  l  ?C  \  ^  avec  ^ S  •  Un  exemple  de  montage  expe¬ 
rimental  est  donne  sur  la  figure  18,  on  distingue  les  deux  lasers  de'  frequence  fixe  et  de  frequence 

variable  uu  5  .  Mais  il  existe  des  difficultes  experimentales  inherentes  &  la  nature  non  lineaire  de  la 
DRASC  : 

-  L' interaction  n'est  efficace  au  point  de  mesure  que  s'il  y  a  superposition  spatiale  et  temporelle  des 
deux  faisceaux  excitateurs,  ce  qui  demande  un  alignement  soigne  de  ces  deux  faisceaux. 

-  Une  spectrocopie  Raman  I  haute  resolution  exige  une  grande  stabilite  en  frequence  des  deux  lasers. 

-  Des  monochromateurs  de  grand  pouvoir  de  rejection  deviennent  necessaires  pour  isoler  le  rayonnement 
anti-Stokes  3  detecter  des  rayonnements  excitateurs. 

-  La  puissance  anti-Stokes  variant  rapidement  avec  les  puissances  des  deux  lasers  il  faut  affranchir  la 
mesure  de  leurs  fluctuations  en  puissance. 

Des  mesures  de  grandeurs  instantanees  sont  possibles  avec  la  technique  DRASC  mais  pour  1' instant  on 
obtient  une  mesure  des  fluctuations  de  la  masse  volumique  avec  les  performances  suivantes  :  duree  de  la 
mesure  10  ns  H  la  cadence  de  10  Hz  avec  une  precision  dans  la  mesure  de  la  masse  volumique  P  de  3%  ;  or. 
obtient  ainsi  pour  la  mesure  de  P  la  valeur  moyenne  et  l'ecart  type. 

Toutes  ces  techniques  progresseront  si  les  techniques  des  lasers  progressent  egalement  [16]  [17]. 

4.4  -  Ji6uata>cition  at  diagnostic  dz  V zc.oalmo.nt  pan.  dzi>  techniques  dz  cUj&iuion, 

emiAAiovi  eX  abAotiptcon 

§P*d£_globa£e_du_SLayqnn<ime.nt_di66u6&  [IS J 

Des  methodes  globales  ont  ete  utilisees  pour  determiner  les  variations  de  la  masse  volumique  & 
partir  de  la  lumiere  diffusee  ;  un  tel  montage  est  donne  sur  la  figure  19  :  un  faisceau  laser  est  place 
dans  ce  cas  perpendiculairement  S  un  jet  turbulent  qui  se  melange  avec  l'air  ambiant.  La  lumiere  diffusee 
est  analysee  par  un  photomul tipi icateur  que  l'on  repere  par  un  angle  ©  ,  l'origine  de  1’ angle  etant  le 

faisceau  laser.  De  cette  lumiere  diffusee  par  le  melange  turbulent  on  en  deduit  la  loi  angulaire  de 
l’intensite  (indicatrice  de  diffusion)  et  on  analyse  le  spectre  de  la  lumiere  diffusee,  spectre  au  sens 
propre  du  terme  ainsi  que  les  spectres  des  fluctuations.  Les  fluctuations  de  masse  volumique  resultent 
ainsi  de  trois  fluctuations  (pression,  temperature  et  concentration) 

Cette  methode  globale  parait  interessante  mais  1 ' interpretation  reste  delicate  comme  en  temoigne  les 
travaux  de  Lelievre  et  Picard. 

Pour  reconstituer  l'effet  tridimensionnel  du  phenomene,  R.  Goulard  [19]  a  propose  une  technique 
generale  en  emission  et  absorption  pour  effectuer  des  mesures  dans  un  espace  tridimensionnel.  Il  s'agit 
d’une  methode  £  balayage  mult iangulaire  telle  que  celle  representee  sur  la  figure  20  oh  le  profil  en 
fonction  d'une  variable  spatiale  peut  etre  deduit  de  la  variation  de  1 'angle  et  de  1 'angle  Q  .  Ce 

champ  de  recherche  peut  encore  etre  ameliore  en  rendant  la  mesure  instantan£e. 

Ce  dernier  aspect  conduit  vers  une  solution  du  probleme  qui  consiste  &  obtenir  les  variations  d'un 
parametre  dans  1' espace  et  dans  le  temps. 

4.5  -  {/^uatc&dtcon  da  la  vaAA.(itLov\  dz  P  pa/i  doj>  mGthodeA  photo gKapkigue. 6 

Nous  avons  indique  au  debut  de  ce  rapport  1' importance  que  revet  la  connaissance  des  structures  et 
l'interet  que  l'on  attache  &  la  topologie  des  ecoulements  turbulents.  Les  methodes  utilises  actuellement 
pour  mieux  decrire  les  structures  posent  des  problemes  difficiles  5  resoudre  car  les  informations  enregis- 
trees  doivent  etre  convenablement  interceptees .  Les  techniques  photographiques  ou  cinematographiques 
sont  : 

-  1 ' utilisation  de  traceurs  (solides  ou  f luides) , 

-  la  strioscopie, 

-  1 ' ombroscopie , 

-  l'holographie . 

La  strioscopie  et  la  methode  des  ombres  donnent  les  fluctuations  de  masse  volumique  de  la  turbulence, 
mais  elles  integrent  ces  fluctuations  sur  tout  le  trajet  du  rayon  lumineux.  Les  techniques  du  traceur  et 
l'holographie  sont  interessantes  car  elles  identifient  la  position  de  la  particule  fluide  mais  elles  sont 
d' interpretation  complexe  lorsque  l'on  desire  identifier  la  structure  de  la  turbulence.  La  methode  holo- 
graphique  contient  les  informations  necessaires,  mais  la  tridimensionnalite  stochastique  du  proc£d6  rend 
1' interpretation  de  ces  informations  difficiles.  Nous  d€crirons  une  methode  qui  nous  parait  interessante 
pour  l'etude  de  l'ecoulement  turbulent,  elle  consiste  3  utiliser  pour  la  visualisation,  la  luminescence 
de  substances  organiques  comme  la  2-3  Butanedione  [20].  Elle  permet  ainsi  d'introduire  un  gaz  phosphores¬ 
cent  pour  visualiser  l'ecoulement.  On  peut  de  cette  mani£re  etudier  les  processus  de  diffusion  de  la 
turbulence  ainsi  que  les  phenom£nes  de  melange  turbulent. 
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Cette  technique  consiste  done  &  detecter  le  transfert,  par  voie  optique,  de  l'energie  d'une  molecule 
&  une  autre  par  collisions.  Ces  collisisons  coincident  avec  le  melange  des  especes  &  l'echelle  moleculaire. 


La  molecule  est  excitee  par  un  rayonnement  ayant  une  certaine  distribution  spectrale,  la  lumiere 
incidente  est  absorbee  et  reemise  sous  forme  de  fluorescence  pour  un  etat  singulet  excite  et  sous  forme  de 
phosphorescence  si  1' emission  provient  d'un  triplet  excite.  Sur  le  plan  pratique  ceci  se  traduit  par  un 
mecanisme  represente  sur  la  figure  21  :  la  lumiere  incidente  est  composee  de  plusieurs  longueurs  d'onde  ; 
elle  est  filtree  pour  n'emettre  que  dans  le  bleu,  le  rayonnement  est  absorbe  par  la  butanedione  excitee  et 
la  molecule  eraet  dans  le  vert,  e'est  la  phosphorescence,  dans  le  bleu  e'est  de  la  fluorescence  qui  s'ajoute 
au  bleu  de  la  lumiere  incidente.  On  peut  par  la  suite  enregistrer  sur  une  plaque  sensible  le  rayonnement 
vert  de  la  phosphorescence  qui  est  tres  intense. 

Les  avantages  de  cette  technique  sont  les  suivants  : 

-  II  est  possible  de  distinguer  le  melange  du  aux  remous  crees  par  la  turbulence  du  melange 
du  3  la  diffusion  moleculaire. 

-  On  peut  faire  une  mesure  directe  de  concentration  du  traceur. 

-  On  peut  isoler  une  partie  de  1 ' ecoulement  et  l'analyser. 

-  On  possede  egalement  une  vue  d* ensemble  des  phenomenes. 

-  On  enregistre  la  structure  du  melange  figee  a  un  instant  donne. 

II  existe  cependant  des  inconvenients  3  cette  technique,  en  particulier  elle  n'est  valable  qu'3 
basse  temperature  car  il  faut  trouver  des  substances  pour  lesquelles  1' emission  de  phosphorescence  ou  de 
fluorescence  demeure  valable  dans  la  gamine  des  temperatures  etudiees  : 

-  Si  l'on  desire  obtenir  des  photographies  instantanees  la  vitesse  ne  peut  etre  trop  rapide  et  de  ce  fait 
on  mesure  dif f icilement  le  champ  des  vitesses. 

-  L'examen  est  bidimensionnel  et  le  phenomene  est  tridimensionnel . 

-  L'analyse  qualitative  des  films  est  relativement  facile  mais  tres  difficile  lorsqu'on  veut  obtenir  des 
resultats  quantitatifs. 


5  -  MESURE  VES  FLUCTUATIONS  VE  PRESSIOU 

La  mesure  dans  un  ecoulement  d'une  pression  moyenne  d' arret  ou  statique  est  encore  du  domaine  du 
possible,  mais  mesurer  avec  precision  une  pression  fluctuante  dans  un  ecoulement  reste  encore  du  domaine 
du  futur. 


D'un  point  de  vue  theorique  la  pression  apparait  tout  d'abord  dans  1 'equation  d'etat,  de  la  forme, 
pour  un  milieu  &  composition  variable  : 

'V  -  r  *  T  ^ 

etant  la  fraction  massique  de  l'espece  j  . 

Elle  apparait  egalement  dans  1' equation  du  mouvement  et  dans  1' equation  du  bilan  d'energie.  Les 
problemes  de  fermeture  font  apparaitre  des  termes  de  correlation  avec  la  pression  et  il  est  souvent  tres 
difficile  d'apprecier  1' importance  de  ces  termes  qui  sont  de  la  forme  : 

Q  '  /  ,  1) 


qui  apparaissent  dans  le  bilan  d'energie  cinetique  de  turbulence.  Il  faut  comparer  par  exemple  le  terme 
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oQ  i  represente  les  fluctuations  de  l'Snergie  de  dissipation,  les  fluctuations  du  flux  massique 
de  diffusion  dans  la  direction  $  ,  w les  fluctuations  du  taux  massique  de  production  chimique  j  * 
la  chaleur  de  formation  de  l'espece  ^  les  fluctuations  du  flux  de  chaleur  dans  la  direction/^. 

Sur  le  plan  theorique  et  sur  le  plan  de  la  modelisation  des  phenomenes  turbulents,  la  mesure  des 
fluctuations  de  pression  revet  done  une  tres  grande  importance. 


Il  n' existe  actuellement  que  des  methodes  de  mesure  par  sonde  et  des  mesures  indirectes. 

Me4  uAe  _  po/L  _A 

On  dispose  bien  de  sondes  de  pression  de  petites  dimensions  (de  l'ordre  du  mm  de  diametre)  de  temps 
de  reponse  de  10  kHz  mais  qui  perturbent  1 'ecoulement ,  e'est  loin  d'etre  1' ideal  pour  faire  une  mesure  de 
pression  correcte  dans  un  ecoulement  turbulent^  Elies  sont  en  general  couplees  d  des  mesures  de  vitesse 

instantanees  et  des  correlations  du  type  ( TTtt  )  ont  et£  obtenues.  Ceci  a  lieu  dans  des  ecoulements 
isothermes  car  ces  sondes  sont  encore  tres  sensibles  aux  gradients  de  temperature.  Les  acousticiens  utili- 
sent  de  telles  sondes  pour  analyser  le  champ  acoustique  du  bruit  engendre  par  les  ecoulements. 


17 

MethodeA  _  CndiAecteA 

Elies  sont  de  deux  sortes  : 

a)  Les  fluctuations  de  pression  locales  determinees  a  partir  des  mesures  de  pressions  instantanees  A  la 
paroi,  une  analyse  theorique  plus  ou  moins  complexe  permet  de  passer  des  mesures  des  fluctuations  a  la 
paroi  aux  valeurs  de  la  pression  locale  instantanee. 

b)  Des  mesures  de  masse  volumique  et  de  temperature  en  un  meme  point  d'oG  l'on  deduit  les  fluctuations  de 
pression  en  ce  point. 

Comme  on  peut  s'en  rendre  compte  les  mesures  de  pression  dans  les  ecoulements  turbulents  sont  encore 

tres  imprecises.  C'est  un  probleme  difficile  qui  est  loin  d'avoir  requ  des  solutions  satisfaisantes . 

6  -  CONCLUSION  SUR  LA  MESURE  VES  SCALAIRES  EM  ECOULEMENT  TURBULENT 

Ce  rapport  comprend  cinq  parties  : 

a)  La  premiere  partie  est  consacree  aux  objectifs  de  ce  groupe  de  travail  qui  ont  largement  evolue  au 
cours  du  temps.  On  est  parti  d'une  idee  assez  vague  du  Professeur  A.  Ferri  qui  portait  sur  les  mesures 
des  fluctuations  de  pression  et  sur  1 ' interpretation  de  ces  mesures  dans  un  ecoulement  de  gaz  chaud  et 
dans  les  flammes. 

Ce  theme  se  rapportant  A  un  probleme  tres  complique  n'a  pas  requ  1  *  approbation  de  tous  les  membres 
du  groupe  de  travail,  on  a  done  change  d'objectif  et  defini  un  nouveau  theme  relatif  cl  la  mesure  des 
scalaires  dans  un  ecoulement  turbulent  chaud  ou  dans  un  ecoulement  en  combustion  et  portant  sur  la 
mesure  des  temperatures,  des  concentrations  ou  de  la  masse  volumique  et  de  la  pression. 

b)  La  deuxieme  partie  detaille  comment,  H  partie  de  ces  mesures,  il  sera  possible  de  mieux  connaitre  la 
structure  turbulente  des  zones  de  melange  et  des  flammes  et  de  discemer  en  particulier  la  coherence  de 
ces  structures. 

On  montre  l'interet  qu'il  y  a  a  utiliser  des  techniques  photographiques  capables  de  visualiser  ces 
structures  et  le  role  que  peut  iouer  dans  la  description  de  1 'ecoulement  les  spectres  de  turbulence 
obtenus  a  partir  des  mesures  fluctuantes  de  temperatures  et  de  concentrations  puisqu'il  n'est  pas,  pour 
1* instant,  possible  d'exploiter  les  mesures  des  fluctuations  de  pression.  Ces  resultats  sont  importants 
puisqu'ils  seront  &  la  base  de  tout  progres  dans  la  modelisation  des  ecoulements  turbulents  et  des 
flammes  turbulentes  suivant  qu'il  s'agit  de  flammes  premelangees  ou  de  diffusion,  de  flammes  ouvertes 
ou  confinees. 

c)  La  troisieme  partie  est  consacree  A  la  mesure  instantanee  de  la  temperature.  Deux  techniques  sont 
decrites  :  les  methodes  optiques  et  les  methodes  par  sondes.  Nous  avons  insiste  sur  les  difficultes 

d' interpretation  des  methodes  optiques  car  elles  integrent  les  phenomenes  sur  une  portion  du  faisceau, 
ces  interpretations  se  font  avec  de  nombreuses  hypotheses.  Les  methodes  par  sonde  necessitent  la 
connaissance  de  donnees  permettant  d'evaluer  &  partir  de  la  mesure  par  sonde,  les  fluctuations  reelles 
de  la  temperature  du  gaz.  Ces  methodes  sont  interessantes  si  elles  sont  couplees  avec  des  sondes  de 
mesure  de  la  vitesse. 

d)  La  quatrieme  partie  traite  des  mesures  instantanees  des  concentrations.  Ces  techniques  utilisant  les 
performances  des  lasers  sont  en  pleine  evolution.  Trois  methodes  sont  plus  particulierement  decrites  : 
la  diffusion  Raman  spontanee  (DRS) ,  les  methodes  utilisant  le  phenomene  perturbateur  de  la  fluorescence 
et  la  methode  de  diffusion  Raman  anti-Stokes  coherente  (DRASC) .  De  tres  nombreux  travaux  ont  montre  que 
ces  trois  methodes  sont  interessantes  pour  determiner  la  concentration  moyenne,  mais  peu  de  ces 
methodes  sont  adaptees  &  la  mesure  instantanee  de  la  concentration.  Seule  la  DRS  donne  dans  la  version 
impulsionnelle  les  concentrations  instantanees  et  la  fluorescence,  mais  dans  un  domaine  de  temperature 
etroit,  la  DRASC  pour  1' instant  ne  determine  que  les  fluctuations  de  masse  volumique.  Dans  ce  meme 
paragraphe  nous  avons  discute  des  methodes  globales  de  diffusion  et  des  techniques  appropriees  pour 
traiter  le  probleme  tridimensionnel .  Nous  avions  prevu  initialement  un  paragraphe  sur  le  couplage  de 
ces  scalaires  avec  la  mesure  de  la  vitesse  instantanee,  nous  1 'avons  supprime  car  celui-ci  est  discute 
dans  chaque  partie  et  nous  avons  inclus  les  techniques  de  visualisations  qui  reposent  en  fait  sur  des 
mesures  de  masse  volumique.  Nous  avons  donne  quelques  indications  sur  la  technique  des  gaz  traceurs 
phosphorescents  et  sur  les  avantages  et  les  inconvenients  de  cette  technique. 

e)  La  cinquieme  partie  est  consacree  3  la  mesure  instantanee  de  la  pression.  Cette  mesure  est  tres  diffi¬ 
cile  "in  situ".  II  faut  passer,  pour  l'atteindre,  A  une  modelisation  a  priori  de  la  zone  de  mesure. 
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TURBULENT  TRANSPORT  PHENOMENA 


by  Marcel  BARRERE 

Working  group  09  report 

Office  National  d‘ Etudes  et  de  Recherches  Aerospatiales  (ON ERA) 
92320  Chat  i I  Ion  (France) 


1  -  INTRODUCTION 

This  PEP  Working  Group  has  been  initiated  by  Professor  A.  Ferri  and  we  would  like,  to  put  into 
situation  the  activity  of  this  Group,  to  recall  Prof.  Ferri' s  comments  on  this  problem  of  transport 
phenomena  in  a  turbulent  medium  : 

"The  relation  between  transport  properties  and  unsteadiness  of  the  flows  is  not  understood  ;  the  progress 
performed  is  not  sufficient  -  usually  the  new  approaches  complicate  the  analysis  without  improving  capability. 
Several  approaches  have  been  suggested  in  substitution  for  predicting  flow  properties  of  the  Prandtl  type 
eddy  viscosity  : 

a)  Use  of  average  fluctuation  quantities  and  transport  properties  related  to  kinetic  energy.  No 
rational  is  available  for  the  relation.  Several  independent  relations  have  been  suggested.  In  addition,  the 
relation  does  not  take  into  account  the  importance  of  large  eddies  with  respect  to  small  size  eddies  which  is 
very  important  and  is  the  basis  of  the  Prandtl  eddy  viscosity  model. 

b)  Use  of  time-dependent  properties.  (Intermittency) .  Such  an  approach  is  more  physical  ;  however, 
the  methods  to  determine  transport  properties  from  the  gross  physical  properties  of  the  flow  are  not  yet 
available.  It  is  known  that  several  parameters  affect  the  transport  properties  :  the  scale  intensity  and 
distribution  of  the  turbulence  at  the  initial  station  ;  the  generation  of  turbulence  due  co  shear  and 
energy  changes  ;  the  decay  of  turbulence,  and  the  change  of  frequency  and  scale  of  the  turbulence  ;  the 
generation  of  noise  due  to  turbulence  that  produce  pressure  fluctuations  and  diffuses  by  a  wave 
mechanism  and  is  a  source  of  intermittency. 

Little  quantitative  information  presently  exists  on  such  parameters  and  no  information  exists  on 
turbulence  in  three-dimensional  flow.  The  panel  recommends  that  a  substantial  effort  be  organized  to  develop 
better  knowledge  on  such  quant ites.  The  shortcomings  in  physical  knowledge  make  the  results  of  numerical 
methods  inaccurate  and  often  their  usefulness  and  capability  are  dangerously  overstated.  The  problem  is  still 
more  complex  for  two  phase  flow.w 

So  the  Working  Group's  purpose  was  to  improve  our  understanding  of  the  flow  structure,  with  a  view 
to  arrive  at  a  prediction  of  transport  phenomena  in  a  turbulent  flow.  We  brought  to  light  the  conclusion 
that  to  reach  this  stage  of  understanding  it  was  very  important  to  possess  better  analytical  and 
measuring  tools  and  that,  in  particular,  we  place  in  priority  the  measurement  of  scalar  parameters 
and  of  their  spatio-temporal  fluctuations.  We  already  have  at  our  disposal  safer  techniques  for 
measuring  fluctuations  of  the  velocity  vector,  mainly  laser  techniques. 

This  Working  Group  selected  four  scalar  parameters  :  concentration,  temperature,  pressure  and 
density,  which  seem  the  most  significant.  Prof.  Ferri 's  wish  had  been  to  emphasize  pressure  fluctuation 
measurements,  but  as  it  appears  that  this  parameter  is  difficult  to  measure  and  that  there  does  not  exist, 
for  the  moment,  any  direct  means  for  measuring  it,  the  Working  Group  preferred  to  extend  this  suggestion 
to  the  measurement  of  scalars  in  general.  These  measurements,  of  a  more  general  nature,  should  allow  a 
better  understanding  of  turbulent  flows  with  or  without  chemical  reactions  and  in  a  medium  whose  density 
p  might  fluctuate  (compressible  or  incompressible  flow). 

The  knowledge  of  physical  and  chemical  phenomena  in  a  turbulent  medium  has  direct  applications  on 
methods  of  quantitative  calculation  of  flows  in  combustors  and  chemical  reactors  ;  it  is  the  case  in 
particular  of  the  combustion  chambers  of  main  and  reheat  combustors,  of  ramjet  combustors,  of  liquid  and 
solid  propellant  rocket  combustors,  of  chemical  laser  chambers,  of  reactions  in  the  wake  of  turbojets 
and  rockets,  of  plain  or  tubular  reactors  encountered  in  chemical  engineering,  of  the  evolution 
of  the  atmospheric  and  stratosphere  composition. 

Before  dealing  with  the  contents  of  this  report,  we  would  like  to  emphasize  two  difficulties  that 
imposed  the  orientation  chosen  : 

a  -  Chemical  reactions  take  place  if  reacting  species  are  mixing  at  a  very  small  scale  :  turbulence 
intervenes  only  by  the  way  in  which  it  brings  the  species  close  to  each  other,  and  may  have  a  beneficial 
effect  if  it  brings  rapidly  hot  burnt  gases  into  contact  with  premixed  fresh  gases  ;  this  effect  is 
double,  and  concerns  both  the  temperature  and  the  concentration  parameters.  We  thus  understand  why,  in 
this  cases,  turbulence  increases  the  chemical  production  rate.  But  this  mixing  effect  by  turbulence  may, 
in  some  cases,  reduce  the  production  rate  ;  it  is  the  situation  of  a  diffusion  flame  for  which  oxidizer 
and  fuel  are  injected  separately  into  the  combustor  ;  the  turbulence  structure  may  then  be  such  that  some 
oxidizer  and  fuel  molecules  do  not  meet,  which  limits  the  chemical  production  rate.  We  may  thus  understand 
by  which  mechanism  turbulence  modifies  the  chemical  reaction  rate,  but  we  understand  less  well  how  the 
che.nical  reaction,  by  the  energy  it  provides  locally  through  very  rapid  increase  of  volume  and  temperature 
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(decrease  of  density),  can  modify  turbulence  ;  that  is  why  temperature  and  concentrations  measurements  in 
time  and  space  are  necessary. 

b  -  However,  there  remains  that  the  problem  of  measuring  pressure  and  its  fluctuations  in  turbulent 
flow  is  still  without  a  solution.  Mean  pressure  measurements  in  such  a  flow  are  often  uncertain,  as  it  is 
necessary  to  avoid,  eliminate  or  calibrate  the  probe  sensitivity  to  pressure  fluctuations  as  well  as  its 
parasitic  sensitivity  to  pressure  fluctuation  direction  and  intensity  ;  in  particular  we  should  determine 
in  which  conditions  the  flow  around  the  probe  may  be  considered  as  quasi  steady.  Usually  we  assume  that  a 
calibration  performed  in  steady  flow  can  be  used.  Verifications  in  different  experimental  conditions 
should  give  an  idea  of  the  accuracy  of  the  pressure  fluctuation  measurement  p’  ;  in  particular,  these 
checking  experiments  should  include  measurement  of  p'  behind  a  grid  generating  an  isotropic  turbulence  ; 
measurements  of  p'  can  be  deduced  from  velocity's  correlations.  This  measure  of  pressure  fluctuations  meets, 
as  will  be  seen  later,  with  many  difficulties. 

After  these  preliminary  remarks,  the  report  comprises  four  parts  : 

i  -  a  general  introduction  concerning  the  turbulent  flow  structure  ; 

ii  -  temperature  measurements,  mainly  methods  that  seem  best  adapted  to  turbulence  studies  ; 

iii  -  concentration  measurements,  with  particular  emphasis  on  new  laser  techniques  ;  density 

measurements  and  visualization  processes  ; 

iv  -  lastly  pressure  measurements . 

We  mention  in  the  conclusion  the  orientations  that  could  be  given  on  the  scientific  view  point,  and 
the  recommendations  to  be  proposed  to  PEP. 

2  -  TURBULENT  FLOW  STRUCTURE 

To  better  model  turbulent  flows  we  must  be  capable,  from  vectorial  and  scalar  measurements,  to  better 
describe  and  analyze  them.  The  first  difficulty  arises  from  the  fact  that  there  is  not  one  turbulent  flow 
but  a  great  variety  of  flows  that  may  comprise  coherent  structures  and  random  structures.  One  of  the  first 
tasks  of  the  analysis  is  to  bring  to  light  these  coherent  structures,  and  a  major  part  of  the  measuring 
activity  concerns  this  research. 

To  emphasize  the  kind  of  solutions  that  can  be  offered  by  diagnostic  methods,  we  shall  take  a  few 
examples  underscoring  first  the  difficulties  pertaining  to  the  measurement  of  vectorial  or  scalar 
parameters  : 

-  the  measurement  is  made  at  a  fixed  point,  and  it  is  often  difficult  to  follow  tracers  placed  within 
the  flow  ; 

-  the  measurement,  if  made  by  an  optical  process,  integrates  over  the  whole  beam  path  ;  but  if  we 
introduce  a  probe,  this  may  perturb  the  flow  ; 

-  it  is  also  mandatory  that  the  measuring  system  as  a  whole  allows  a  sufficiently  short  measuring  time 
in  order  to  follow  the  least  fluctuation  ; 

*  it  is ‘of ten  difficult  to  place  at  the  same  point  several  measuring  devices  so  as  to  precisely 
reconstitute  the  flow  structure  in  both  time  and  space,  as  there  is  not  enough  information  available. 

So  we  find  ourselves  faced  with  a  complex  scalar  and  vectorial  field,  resulting  from  the  superposition 
and  interaction  of  elementary  fields. 

Let  us  consider  for  instance,  and  this  will  be  the  first  example,  the  structure  of  a  diffusion  flame, 
made  by  a  jet  of  cold  oxidizer  (air)  and  a  parallel  jet  of  fuel  (hydrogen)  concentric  with  the  former.  By 
means-  of  optical,  or  hot  wire,  methods,  we  intend  to  study  : 

a)  the  spatial  distribution  field  of  the  oxidizer  flow  velocity  : 

comprising  two  parts  :  a  mean  value  u  and  a  fluctuating  value  u'  ; 

b)  the  density  field,  which  fluctuates  because  of  the  turbulence  of  the  two  mixing  flows,  oxidizing 
gas  on  the  one  hand,  injected  fuel  gas  on  the  other,  the  latter  provoking  a  density  variation  : 

^ 1  Ic*- yyA)+ e; 

c)  the  space-dependent  density  field  due  to  the  mixing  of  combustion  products  with  fresh  oxidizer 
and  fuel  species  : 

%  =  %  A)  y*) 

d)  lastly,  the  space-dependent  density  field  dee  to  chemical  reactions  of  combustion,  which 
liberate  energy  : 

.4) 4<*f  («,  yy*) 

These  elementary  fields  can  be  studied  separately,  and  their  structure  can  be  described  in  detail  from 
the  correlation  functions  and  energy  spectra,  expressed  as  a  function  of  wave  number  [ 1 1 . 
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These  functions  can  express  the  mechanism  of  transport  and  dissipation  if  we  consider  the  superposition 
and  interaction  of  these  elementary  fields.  One  operating  procedure  consists  in  studying  separately  the 
structure  of  each  elementary  field  and  after  the  complex  structure  of  the  resulting  field.  It  is  important 
first  to  know  the  field  ,  as  the  concentration  and  temperature  scalar  quantities  are  transported  by 


the  turbulence  due  to  velocity  fluctuations.  In  this  particular  case,  the  analysis  shows  that  turbulence  is 
isotropic  and  follows  the  Kolmogorof  law  of  viscous  dissipation. 

The  field  of  cold  mixing  of  the  hydrogen  jet  in  an  air  flow  can  be  studied  thanks  to  an  important 

result  provided  by  the  analysis  of  energy  spectra,  such  as  those  represented  on  figure  1,  as  well  as  the 

field  in  the  combustion  zone.  We  note  in  particular,  in  this  figure,  the  difference  of  spectra  with  and 

without  combustion.  Let  us  remark  first  that  the  small  structures  corresponding  to  large  values  of  the  wave 
number  k  have  slopes  different  from  those  of  Kolmogorof  :  -3  instead  of  -5/3  ;  the  slope  for  viscous 
dissipation  is  steeper  for  scalars  than  that  corresponding  to  velocity. 


The  field  due  to  chemical  reaction  has  a  structure  very  different  from  that  of  the  other  fields 

(  >  5  ,  T  )  .  The  spectrum  presents  a  maximum  around  1  mm  for  K"1 .  This  may  indicate  cells 

functioning  as  micro-combustors,  and  an  increased  stability  fo  chemical  reaction  zones  of  this  size.  Thus, 
this  concept  excludes  a  mechanism  of  combustion  along  a  flame  front. The  analysis  of  turbulent  energy  budget, 
defined  as  the  integral  relative  to  all  the  frequencies  of  the  energy  spectrum,  provides  a  complement  of 
information  on  the  exchange  mechanism  in  flames.  The  law  of  energy  distribution  in  the  flame  according  to 
wave  number  and  the  decrease  of  total  energy  may  result  from  the  decrease  of  the  weight  of  the  Fourier 


coefficients  for  the  lowest  frequencies  (K  <  0.3  mm"1),  and  it  is  the  Fourier  coefficients  situated  around 
K  =  0.9  ram”1  that  have  the  greatest  weight. 


Figure  2  (ref.  1)  also  shows  that  the  cold  hydrogen-air  mixing  £  c  t K  )  )  and  the  fresh  gas  -  burnt 

gas  mixing  have  identical  spectra,  as  well  as  the  global  spectrum  £  CK  )  ,  with  two  slopes  :  one 
corresponding  to  K"1*5  for  large  structures  and  one  corresponding  to  K"3  for  small  structure,  showing  an 
identical  mechanism  for  the  hydrogen  -  air  mixing  and  the  fresh  gas  -  burntgas  mixing. 


We  think  that  results  of  this  type  are  very  useful  for  studying  flow  structure,  mixing  processes  and 
combustion  processes. 


The  second  example  we  chose,  that  shows  particular  interest  of  these  scalar  measurements,  concerns  the 
study  of  transport  mechanisms  and  of  the  scalar  field  structure  in  the  case  of  the  boundary  layer  of  the 
mixing  between  a  jet  of  hot  air  and  the  surrounding  medium  at  rest,  a  study  carried  out  by  Fiedler  at  the 
Berlin  University  [2j.  An  analysis  of  the  velocity  and  temperatures  field  makes  it  possible  to  evaluate  the 
equation  of  temperature  fluctuations  (  T'2  )  in  the  flow  field,  by  making  apparent  the  terms  of  convection, 

production,  diffusion  and  dissipation,  as  shown  on  figure  3  which  gives,  as  a  function  of  the  ratio  H  —  ^  ~ 

'  X  -  oc4 

(x  in  the  flow  direction,  y  perpendicular),  the  values  of  the  process  that  brings  a  gain  or  a  loss  :  the 
values  of  the  production  is  the  most  important  component  and  posses  two  maximums,  convection  of  the  mixing 
layer  where  diffusion  reaches  a  maximum,  dissipation  is  not  measured  and  is  obtained  by  difference.  This 
study  makes  it  possible  to  appreciate  the  distribution,  the  terms  v'T'  and  u'T'  corresponding  to  lateral 
and  longitudinal  diffusivity  and  v'T'2. 

The  mixing  process  has  also  been  analyzed  by  perturbing  the  flow  with  periodic  waves  issued  from  a 
loudspeaker  placed  at  the  nozzle  trailing  edge.  The  author's  conclusion  is  : 

-  Here  exists  a  strong  influence  of  external  perturbations  of  small  amplitude  on  the  flow  structure, 
especially  on  the  scalar  structure,  expressed  by  the  formation  of  cells  of  coherent  structure  ; 

-  in  the  case  of  unperturbed  flow,  coherent  structures  seem  to  be  rather  unimportant  ; 

-  regions  with  negative  production  are  present  in  the  case  of  excited  flow  ; 

-  the  diffusion  gradient  model  is  not  applicable  ; 

-  the  excitation  not  only  disciplines  the  large  vortices  but  also  influences  the  structure  of  internal’ 
turbulence. 

This  work  illustrates  well  what  can  be  deduced  from  the  study  of  a  scalar  field  coupled  with  that  of  a 
velocity  field  with  a  view  to  modelling  the  flow  structure  and  deducing  the  transport  processes. 

The  third  example  we  chose  concerns  the  observation  of  coherent  structures  in  combustion  by  photographic 
techniques,  which  in  fact  reveal  the  density  field.  In  this  respect,  a  very  interesting  work  has  been 
carried  out  by  Yule,  Chigier  and  Thompson  £3j.  In  many  cases,  the  role  of  coherent  structures  is  more 
easily  obtained  in  combustion  zones  in  view  of  the  luminous  character  of  flames  ;  it  is  in  particular  a 
simple  way  for  visualizing  fuel-oxidizer  mixing.  Turbulent  diffusion  flames  present  a  great  structural 
variety,  and  in  most  cases  there  appear  vortex  rings  similar  to  those  observed  in  non-reactive  mixings  ; 
the  flame  confinement  may  change  this  structure,  and  the  length  of  the  vortex  ring  zone  is  greater  with 
combustion  than  without.  We  note  a  coalescence  of  these  rings  in  both  situations,  with  and  without 
combustion.  These  rings  also  appear  in  the  combustion  of  jets  of  pulverized  fuel,  and  the  existence  of  these 
coherent  eddies  depends  on  the  distribution  of  the  oxidizer-fuel  mixing.  Figure  A  gives  the  schematics  of  an 
eddy  vortex  appearing  in  a  diffusion  flame  and  separating  fuel  from  oxidizer  ;  it  may  represent  the  cross- 
section  of  a  vortex  ring  or,  in  a  narrow  two-dimensional  system,  a  vortex  line,  or  again  the  cross-section 
of  a  helical  vortex.  These  vortices  are  similar  and  roll  up  as  those  observed  in  a  wake  ;  but  they  are  not 
alternative,  but  symmetrical  [4}.  This  schem?.tizat ion  corresponds  to  a  vortex  that  has  not  encountered 
others,  as  such  an  encounter  provokes  instabilities  [5j.  The  zone  separating  fuel  from  air  has  the  shape  of 
a  double  spiral  surrounding  a  vortex  core.  A  flame  takes  place  at  the  boundary  of  this  double  spiral 
separating  fuel  from  air,  and  the  vortex  core  contains  some  fule  and  entrained  air,  and  may  constitute  a  zone 
of  volume  combustion  with  variable  equivalence  ratio  functioning  as  a  recirculation  zone,  and  the  flame  at 
the  boundary  between  the  two  zones  can  be  considered  as  a  laminar  diffusion  flame. 


This  structure  can  be  rather  well  described  if  we  know  the  evolution  of  equivalence  ratio  while  this 
pocket  expands  with  the  increase  of  residence  time  within  the  vortex.  In  the  situation  of  figure  4,  a  flame 
surrounds  completely  the  vortex.  These  vortices  can  meet  other  vortices,  unite  to  form  more  complex 
structures  with  a  three-dimensional  character,  or  form  a  new  structure  after  coalescence  of  several 
vortices. 

Another  type  of  coherent  structure  observed  in  a  flame  is  made  of  large  three-dimensional  elongated 
vortices,  whose  structure  is  shown  on  figure  5  which  describes  an  incomplete  circumferential  coherence 
around  the  flame,  and  also  the  existence  of  a  distribution  of  irregular  vortices,  with  the  presence  of 
smaller  vortices  with  cells  of  oxidizer,  fuel  and  burnt  products.  Fuel  and  air  are  swallowed  by  this  vortex, 
and  mix  more  rapidly  them  in  the  previous  situation.  This  structure  is  observed  in  many  different  flames. 

These  three  examples  have  been  chosen  with  a  view  to  show  what  can  be  deduced  from  the  analysis  of 
scalar  field.  Pressure  fluctuations  have  been  neglected,  as  the  measurement  of  this  parameter  is  much  more 
difficult  than  that  of  others. 

To  decide  what  is  coherent  and  what  is  not,  to  understand  transport  mechanism,  to  recognize  the 
modifications  brought  to  the  turbulent  zone  structure  by  combustion,  to  discover  the  action  of  pressure 
waves  on  the  organization  of  various  phenomena  and  particularly  on  combustion,  to  understand  dissipation 
phenomena,  in  short  bring  forth  information  for  modeling  turbulence,  such  the  results  can  be  expected  from 
diagnostic  methods  in  general,  and  from  methods  for  measuring  scalar  parameters  in  particular. 

3  -  TEMPERATURE  MEASURING  METHODS 


There  exist  a  number  of  methods  for  measuring  temperature  in  a  turbulent  flow  with  or  without  combustion, 
and  these  can  be  divided  into  two  major  classes,  optical  methods  and  probing  methods  ;  we  shall  describe  one 
method  of  each  class  so  as  to  bring  to  light  the  advantages  and  drawbacks  of  each. 

3. 1  -  Optical  methods 

These  methods  present  the  interest  of  not  perturbing  the  flow,  which  is  a  great  advantage  in  case  of 
turbulent  flow,  but  they  integrate  the  phenomena  over  part  of  the  beam  length.  We  can  however  partly  remedy 
this  drawback  by  introducing  locally  into  the  gas  some  tracers  which  have  the  same  temperature,  and  also  local 
velocity,  as  the  gas  itself.  The  technique  we  shall  describe  is  that  of  line  reversal,  but  in  the  infrared 
range  which  makes  measurements  possible  from  600  K  upwards. 

3.1.1  -  Instantaneous  measurement  of  combustion  gas  temperature  and  of  thermal  turbulence  by  infrared 
pyrometry  m  — - - 

A  -  Principle^of _the__measurement 

The  optical  diagram  of  the  device  permitting  the  discrimination  of  signals  issued  from  the  source 
and  the  flame  is  presented  on  figure  6.  The  image  of  the  source  S  is  formed  in  the  plane  C  by  the  lens  L^. 

The  lens  forms  another  image  in  the  flame  plane.  The  optical  system  L^>  which  carries  the  aperture  dia¬ 
phragm  ,  conjugates  the  above  image  planes  on  the  field  diaphragm  d.  The  lens  L,  ensures  the  transport  of 

the  monochromatic  image,  isolated  by  the  monochromator  M,  on  the  prismatic  mirror  P,  which  divides  it  into 
two  almost  equal  parts.  Thanks  to  lenses  L,.,  half  of  the  energy  reaches  the  receiver  Rj ,  the  other  half 

the  receiver  R^.  A  knife  placed  in  plane  C  can  be  displaced  so  as  to  divide  the  image  exactly  at  the 

prismatic  mirror.  So  one  of  the  receivers,  e.g.  R2,  does  not  receive  any  energy  from  the  source,  while  the 
signal  delivered  by  R^  is  unchanged. 
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In  the  absence  of  a  flame  and  before  shutting  off  by  the  knife,  the  detectors  Rj  and  R2  deliver  signals 
and  Y2q,  respectively  proportional  to  the  monochromatic  energy  luminance  from  the  source  L  (.  ^ . 
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where  T",  5  is  the  temperature  luminance  of  the  source  at  wavelength  A  .  We  can  also  bring  about  its 
emissivity  €  t  X)  by  writing  : 

Y.o*  k«*s(MLU,Ts)  .  *20  --  Kj  £s(M  Ux,Ts) 

where  1 5  is  then  its  true  temperature.  In  view  of  the  differences  between  detector  sensitivities  and 
geometric  extend  of  incident  beams,  the  proportionality  coefficients  and  Kj  are  not  identical,  but  it  is 

always  possible  to  make  them  equal  by  acting  in  the  amplifying  stages. 


In  the  presence  of  the  flame  and  with  the  knife  in  place,  R2  gives  a  signal  Y2  due  to  the  flame  alone, 


while  R  gives  a  signal  Y  due  to  the  flame  and  the  source 
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is  the  local  coefficient  of  monochromatic  emission,  equal  to  the  absorption  coefficient  according 
to  the  kirchhof  law,  and  monochromatic  luminance  of  the  black  body  at  the  flame  local 

temperature 


In  this  formulation  we  assume  the  flame  homogeneous  within  a  cross-section  of  the  beam  or,  which  is 
equivalent,  the  beam  infinitely  narrow.  The  first  term  of  Y  corresponds  to  the  transmission  through  the 
flame  of  the  energy  issued  from  the  source  according  to  the  Beer  law,  an  exact  law  in  the  strictly  mono¬ 
chromatic  case.  The  second  term  corresponds  to  the  flame  emission  within  a  thickness  slice  doc  and  to  its 
absorption  by  the  following  slices  from  oc  to  A  ,  2  JL  being  the  flame  thickness.  Only  this  second  term 
appears  in  the  signal.  If  the  flame  is  homogeneous  as  regards  temperature,  integration  is  simple,  as 

L  )  comes  out  of  the  integral  and  the  homogeneity  of  relation  : 

^  L  +  K>,  L(A,TFVf 

*  --  LU,T,)e 

u 
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where  £  is  the  emission  factor,  defined  by  :  ~  ^  ^ 


We  shall  remark  that  £  depends  on  the  flame  thickness.  Knowing  the  source  monochromatic  luminance 
L  (\  TtS  )  makes  it  then  possible  to  calculate  at  any  time,  by  resolution  of  this  two-equation  system, 

the  unknownd  L  ,-T^  )  and  E  (.M  . 


To  perform  this  resolution,  in  the  case  where  Kj  =  ^ 


we  obtain  : 


'ji  =  Ke  Ux,tJ 


K,  we  set  : 


and  we  deduce  : 

L(X, 

The  planck  law  for  wavelength  X  then  makes  it  possible  to  go  back  to  the  flame  temperature  . 

The  response  time  is  only  limited  by  those  of  the  receivers  and  of  the  signal  recorders. 

An  error  calculation  shows  that  the  measuring  precision  is  all  the  higher  as  : 

-  we  work  closer  to  the  reversal  point,  defined  by  the  equality  of  the  flame  and  source  temperatures, 

-  the  flame  emission  factor  £  is  closer  to  unity. 

We  must  emphasize  once  more  that  the  measurement  has  only  a  clear  meaning  if  the  flame  temperature  is 
uniform  over  the  whole  gas  volume  crossed  by  the  optical  beams,  and  that  this  measurement  concerns  the  whole 
volume.  Usually  the  flame  presents  a  temperature  profile,  with  a  maximum  on  the  jet  axis.  In  this  case,  the 
measurement  yields  either  a  temperature  close  to  the  maximum  temperature  if  the  emission  factor  is  small,  or 
a  temperature  close  to  that  on  the  jet  edge  if  the  emission  factor  is  close  to  unity.  In  fact,  when  we  are 
interested  more  particularly  with  the  jet  centre  temperature,  we  show  that  there  exists  an  optimum  for  the 
emission  factor,  an  optimum  which  is  very  well  brought  to  light  by  the  spectrographic  study  of  an  emission- 
absorption  line  which  presents  two  maximums,  symmetrical  about  the  line  centre. 

In  case  of  a  transparent  flame,  the  addition  of  foreign  atoms  (e.g.  sodium)  emitting  on  well  defined 
lines  is  necessary  in  order  to  make  it  monochroraatically  luminous  :  the  gas  temperature  heterogeneity  has 
then  no  more  influence  on  the  measurement  result  if  we  took  care  in  injecting  the  sodium  salt  locally.  But 
we  have  to  make  sure  that  it  has  reached  its  thermal  equilibrium. 

Obviously  the  method  assumes  thermodynamic  equilibrium.  If  it  is  not  the  case,  the  Kirchhoff law  must  be 
modified  as  follows  :  the  spectral  energy  luminance  L  y  of  any  gas  is  equal  to  the  product  of  its  absorption 

factor  by  the  spectral  energy  luminance  L  >  of  the  black  body,  defined  by  the  Boltzmann  law  written 

for  the  populations  and  ^  responsible  for  the  emission  and  absorption  of  the  line  of  wavelength  X  : 

V  AT  ’  ' 

where  and  are  the  respective  statistic  weights  of  levels  1  and  2  of  the  emitting  atom  or 

molecule.  Thus,  we  perform  a  measurement  of  population  temperature.  In  the  normal  case  where  there  exist 
partial  equilibria  corresponding  to  electronic,  vibration  and  rotation  transitions,  the  measurement  provides 
the  equilibrium  temperature  characteristic  of  the  type  of  transition  used.  If,  for  instance,  we  work  on  one 
of  the  lines  of  the  sodium  doublet,  we  shall  measure  an  electronic  temperature. 

Lastly,  the  method  assumes  that  the  emission  and  absorption  phenomena  take  over  all  diffusion  phenomena, 
which  are  recorded  as  a  signal  loss.  In  particular  it  is  possible  to  work  on  particles  in  suspension  within 
the  gas  (bright  flames)  on  condition  their  seize  is  very  small  besides  the  wavelength  used. 

B  -  Transposit ion_of_the_method_to_the_inf rared_range 
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The  instrument  is  made  of  two  parts,  an  emitter  and  a  receiver,  whose  optical  alignment  should  be  as 
perfect  as  possible.  This  instrument  uses  the  own  absorption  and  emission  of  the  flame,  essentially  due  to 
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water  vapor  and  carbon  dioxide.  The  spectral  interval  is  then  no  more  limited  by  a  monochromator,  but  by  an 
interferential  filter  which  integrates  over  the  fine  structure  of  a  vibration-rotation  band.  We  thus  make  a 
systematic  error,  the  monochromatic  calculation  being  no  longer  strictly  valid.  We  shail  show  that  this 
approximation  is  usually  •  sufficient .  The  temperature  homogeneity  of  the  flame  being  assumed,  the  integration 
of  the  relations  giving  Yj  and  Y^  on  the  spectral  interval  A  A  -  gives 
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The  variations  of  the  gas  emissivity  £(M  with  wavelength  are  usually  very  large  ;  on  the  other 
hand,  the  black  body  luminance  for  the  temperature  and  Ts  ,  and  the  source  emissivity  £3  can  be 

considered  as  constant  so  that  they  can  be  taken  out  of  the  integrals  if  the  spectral  interval  is  narrow 
enough.  In  these  conditions,  the  resolution  of  the  system  can  be  performed  as  in  the  monochromatic  case  : 

y*  -  L  (  ^  y  rF  )  /  _  _±_  j  e  C  A  )  A  \  1 

(  ^/vn.  )  L  (k<m>  ^  J 

where  is  the  mean  wavelength  of  the  A  A  interval.  The  emissivity  obtained  this  way  is  obviously  a 

mean  value  £  ,  and  it  will  be  important  only  if  the  spectral  line  density  is  high  enough. 


In  order  to  estimate  the  importance  of  the  systematic  error  due  to  this  integration,  we  consider  the 
case  where  the  variation  of  the  black  body  luminance  over  A  A  is  linear  : 

Ux,  T)  ,  L(k*,T)  +  t 


We  then  calculate 
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As  E (  XN  is  a  function  bounded  by  1 ,  the  corrective  term  has  the  absolute  maximum  value 
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or,  using  the  Planck  law 
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The  error  will  be  negligible  if  the  spectral  interval  is  narrow  enough  and  if  the  mean  emissivity  of 
the  flame  is  not  too  low.  It  will  be  almost  zero  if  the  true  temperature  of  the  source  is  equal  to  that  of 
the  flame. 


The  preliminary  spectrometric  study  if  combustion  flames  shows  that  the  emissivity,  very  weak  otherwise, 
is  very  strong  at  the  4.5  m  wavelength  on  the  hot  emission  bands  of  carbon  dioxide.  That  is  what  we  can  see 
on  the  spectrum  obtained  on  a  stirred  reactor  (fig.  7),  the  thickness  of  the  gas  crossed  being  a  few  centi¬ 
metres.  At  this  wavelength  we  use  an  interferential  filter  of  4%  relative  width.  In  these  conditions,  the 
systematic  error  on  the  luminance  |_  C  A  T"  \  is  lower  than  1%. 

$  J 

The  response  time  of  the  detector-preamplifier  unit  is  of  the  order  of  one  microsecond,  but  it  does  not 
work  on  a  continuous  signal  and  it  is  necessary  to  periodically  interrupt  the  incident  flux  to  reset  the 
zero.  We  use  the  rotation  of  a  cog  wheel  occulting  the  image  obtained  on  plane  <A.  (fig. 6)  at  a  frequency 
of  2500  Hz.  The  two  field  diaphragms  limiting  beams  1  and  2  are  made  of  the  detectors  sensitive  surfaces, 
the  aperture  diaphragm  being  common  to  both  beams.  This  separation  of  the  two  fields  diaphragms  presents  the 
advantage  of  allowing,  by  only  the  axial  displacement  of  the  beam  splitter,  the  disjunction  of  the  two 
images,  which  makes  it  easier  to  occult  one  of  them  with  the  knife. 

The  photograph  of  figure  8  shows  the  instrument  on  an  optical  bench  with,  on  the  foreground,  the 
analogue  computer  which  calculates  on-line  the  flame  luminance  and  emission  factor.  A  particular  difficulty 
in  this  calculation  arises  from  the  2500  Hz  modulation.  The  signals  are  chopped  at  this  frequency,  with  a 
slight  outphasing  between  the  two  detectors,  so  that  the  calculation  result  is  only  valid  during  short  time 
intervals  of  80  us,  separated  by  about  300  us.  So  we  operate  by  sampling,  and  the  broad-band  logic  proces¬ 
sing  unit  developed  for  the  instrument  retains  only  the  correct  measurements  and  displays  zero  for  the  others. 
However,  modulation  is  unnecessary  for  transient  phenomena  lasting  less  than  a  few  hundred  microseconds. 

C  -  Temgerature_f luctuation_measurement 

i  -  Uncoupling  of  temperature  and  emissivity  fluctuations 

We  have  seen  that  the  measurement  has  a  clear  meaning  only  if  the  flame  temperature  is  uniform  over  the 
whole  gas  volume  crossed  by  the  optical  beams.  This  seems  to  preclude  the  possibility  of  performing  thermal 
turbulence  measurements  when  all  the  characteristics  of  the  medium  (concentration,  temperature)  vary  at  all 
points  and  at  all  times  in  a  random  fashion.  We  shall  show  that  it  is  possible  to  characterize  turbulence  in 
case  it  is  statistically  homogeneous  and  when  the  emissivity  fluctuations  are  of  low  amplitude  besides  their 
mean  value. 
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To  prove  this,  let  us  take  again  the  formula  giving  the  energy  received  by  a  detector  from  a  source 

narrow  monochromatic  optical  beam  : 


after  absorption  by  the  flame  and  from  t^e  flame  itself,  this  fo£  a  very 
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We  remark  that  if  the  source  luminance  is  equal  to  the  flame  mean  luminance,  the  emissivity  fluctuations 
disappear  at  the_  first  order.  If  these  fluctuations  are  of  low  amplitude,  it  is  possible  to  neglect  •<  ' 
compared  with  e<  ,  and  there  only  ren^in  the  fluctuations  related  to  temperature  : 
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The  interpretation  of  this  integral  will  be  made  still  easier  if  we  arrange  it  so  that  the  absorption 
coefficient  be  small  ;  we  are  then  able  to  neglect  ^  Lif-  oc  }  besides  unity  : 
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To  this  end,  we  place  ourselves  on  the  2.5  pm  emission  band  for  which  the  emission  factor  is  very  small, 
as  shown  by  the  spectrum  obtained  on  a  stirred  combustor  of  Longwell  type  (fig. 7). 

Lastly,  if  the  temperature  fluctuations  are  of  sufficiently  low  amplitude  it  is  possible  to  linearize  the 
luminance  law  L ( A, T  )  around  the  mean  temperature  T.  ,  and  thus  to  obtain  the  integral  of  temperature 
fluctuations  over  the  optical  path  :  * 
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the  mean  emission  factor. 
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In  case  the  medium  is  not  homogeneous,  i.e.  when  and  £  are  not  independent  of  x,  the 

fluctuation  uncoupling  at  the  first  order  is  no  more  total,  but  the  emissivity  fluctuations  are  much 
attenuated.  The  term  relative  to  thes<*  fluctuations  writes  : 

J  )(  b  loc)  -  bs)  dx 
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as  is  much  lower  than  that  we  would  have  had  without  a  source  (  ®  0)  .  Anyway,  even  when  the  medium  is 

homogeneous  the  quality  of  the  coupling  is  a  function  of  the  adjustment  of  the  source  temperature  to  the 
flame  mean  temperature,  the  latter  having  to  be  determined  experimentally  or  by  a  preliminary  calculation. 

ii  -  Effects  of  fluctuation  integration 

If  the  uncoupling  is  perfect,  the  pyrometric  signal  is  proportional  to  the  integral,  along  the  optical 
path,  of  the  luminance  or  temperature  fluctuations.  What  are  the  effects  of  this  integration,  and  in 
part icular  which  relation  is  there  between  the  signal  variance  (or  standard  deviation)  and  that  of  the 
measured  parameter  ?  We  show  that  they  are  proportional  using  a  formalism  due  to  Uberoi  and  Kovasznay  C7^  , 
when  the  integration  length  is  large  compared  with  the  integral  scale.  The  proportionality  ratio  is  found  to 

be  equal  to  the  ratio  of  twice  the  turbulence  integral  scale  to  the  integration  length  l^L _ ,  i.e.  that 

if  -fl  represents  the  measured  fluctuations  and  T'  the  local  fluctuations,  their  variances  £L/l  and  T'* 
are  related  in  the  folllwing  way  : 
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This  relation  can  be  interpreted  as  follows  :  if  we  divide  the  optical  path  into  n  intervals  of  length 
2  A  :  n  -  U/JJL>  we  can  consider  that  -CL*  is  the  average  of  n  independent  values  T*.  attached  to  each  of 
these  intervals  : 
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The  variance  of  this  sum  is  then  equal  to  the  sum  of  the  variances  and,  as  then  are  variable  with  zero 
average,  we  fin  again  the  above  relation. 
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3.1.2  -  Determination  of  the  turbulence  integral  scale 

To  obtain  the  local  variance  T  *  of  the  measured  parameter,  it  is  necessary  to  know  the  integral 
scale  -A.  .  To  this  end  we  may  use  either  the  spectral  analysis  or  the  orthogonal  beam  method  developed  for 
the  first  time  by  Fisher  and  Krause  [8^  >  but  which  presents  the  disavadvantage  of  requiring  the  implemen¬ 
tation  of  two  measuring  units. 


a)  Sgectra_analysis 


By  spectral  analysis  we  mean  that  of  the  spectrum  obtained  from  the  measured  signal.  Such  a  spectrum 
corresponds,  by  Fourier  transform,  to  the  temporal  correlation  function.  But  the  integral  scale  calculation 
can  only  be  performed  from  the  spatial  correlation  function.  The  passage  from  one  to  the  other  is  only 
possible  if  we  admit  the  Taylor  assumption.  This  assumption  consists  in  supposing  that,  for  high  convection 
velocities,  turbulence  is  frozen.  There  is  then  proportionality  between  the  two  variables,  time  and  a 
spatial  coordinate  parallel  to  velocity.  The  chance  of  phasing  in  time  then  corresponds  to  a  translation  of 
this  coordinate.  The  spectrum  that  we  obtain  this  way  is  the  so-called  "one-dimensional"  spectrum,  which  is 
the  Fourier  transform  of  the  correlation  between  two  points  gradually  separated  from  each  other,  defining 
a  vector  maintained  parallel  to  a  given  direction,  here  that  of  velocity.  But  there  again  the  measured 
spectrum  is  modified  by  the  integration  effect.  In  case  the  integration  length  is  great  besides  ,  we 

demonstrate  the  following  relation  between  the  integrated  spectrum  and  the  local  spectrum  £ ^  : 
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Classically  we  calculate  the  integral  scale  from  the  value  at  origin  of  the  one-dimensional  spectrum 
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The  calculation  of  E^( 0  from  the  spectrum  _f3j  allows  writing 
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and  we  find  again  the  relation  linking  the  variances  : 
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Thus  spectral  analysis,  even  with  the  Taylor  assumption,  does  not  yield  a  complementary  relation  for 
calculating  integral  scale  and  local  variance.  Reciprocally,  the  local  spectrum  can  only  be  determined  but 
for  a  factor.  Let  us  remark  in  passing  the  effect  of  integration  on  the  spectrum  :  if  £4  can  be 


represented  by  a  law  in  \( 
decrease  over  the  whole  spectrum 


is  then  in  K  -(***  <j>he  integration  effect  results  in  a  faster 


b)  Orthogonal _beams 


This  technique  consists  in  measuring  the  correlation  level  between  two  optical  measurements  on  two 
orthogonal  paths  gradually  set  apart  from  each  other,  as  shown  on  figure  9.  Points  A  and  B  limit  the 
shortest  distance  d  of  the  two  straight  lines.  If  we  divide  the  two  optical  paths  into  intervals  of  length 
2  -Au  (taking  A  and  B  as  median  points  of  one  of  these  intervals),  each  of  the  measurements  SZ'^  and  SI'/ 
is  the  average  of  the  independent  parameters  attached  to  each  element  making  up  one  of  the  optical  paths.  We 
understand  that  only  the  two  intervals  including  points  A  or  B  will  bring  a  notable  contribution  to 
correlation  between  and  i  ,  as  the  distance  between  all  pints  of  the  other  intervals  will  be 

greater  than  the  integral  scale  -A.  .  In  fact,  Fischer  and  Krause  ^8]  admit  that  the  correlation  function 
that  we  obtain  by  discrete  values  by  gradually  increasing  d  is  proportional  to  the  correlation 
function  between  points  A  and  B.  It  then  suffices  to  normalize  fi  (oi)  by  the  crossed  beam 

correlation  /i(o)  to  obtain  RW},  then,  by  integration,  the  integral  scale.  The  Fischer-Krause  interpretation 
is  not  actually  correct,  as  it  assumes  as  negligible  any  correlation  between  points  close  to  A  and  B  such 
as  C  and  D.  That  is  why  the  theoretical  study  of  this  technique  has  been  started  again,  still  using  the 
Uberoi  -  Kovasznay  formalism,  valid  for  a  homogeneous  medium.  This  calculation,  for  two  infinitely  thin  beams 
of  length  2  JLa  and  2  along  axes  1  and  2  respectively,  leads  to  the  following  integral  relation 

acd)  _  f  run  m  . 
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In  particular  the  crossed  beam  correlation,  for  d  ■  0,  is  given  by  : 
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where  -A-  is  an  integral  scale  different  from  that  used  up  to  now,  and  defined  by  : 

a'  =  4^-  J-R  <u 
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The  normalization  of  by  allows  one  to  define  the  dimensionless  function  y  : 


whose  inversion  gives  : 
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We  see  that  it  is  necessary,  to  obtain  R  ,  to  trace  V  as  a  function  of  parameter  instead  of  a., 

and  to  take  its  derivative.  The  value  of  this  derivative  at  the  origin  permits  one  to  calculate  ; 
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Lastly,  the  integration  overtfl  of  this  derived  function  yields  a  relation  between  -A.  and  J .  ,  thus 

the  value  of  A  • 


Thus,  the  orthogonal  beam  method  actually  yields  the  integral  scale  ;  better  still,  it  yields  the  two- 
point  correlation  function,  but  it  has  the  disadvantage  of  being  long  (measurements  have  to  be  made  point) 
and  not  very  precise  (function  R  obtained  by  derivation  of  function  Pp  )  ;  so  we  look  for,  and  now 
propose,  a  simpler  measurement  derived  from  the  former  :  it  consists  in  measuring  the  crossed  beam  normalized 
correlation. 

e)  Crossed  beam  normalized  correlation 

It  is  always  difficult  to  perform  absolute  measurements,  and  this  still  holds  for  correlations.  What  is 
given  by  numerical  processing  is  the  correlation  normalized  by  the  r.m.s.  values  of  the  —52^  and  -52^0 
measurements  ;  and  we  already  gave  an  expression  of  these  r.mS.  values  : 
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The  crossed  beam  correlation,  given  by  : 


becomes,  after  normalization  : 
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To  go  further,  we  assume  that  the  two  integral  scales  _A_  and  X  are  close  to  each  other  ;  it  is  the 
case  when  we  admit  an  exponential  law  for  the  correlation  function  R  :  we  then  have  strictly  J\_  =  -A.  .  For 
a  Gaussian,  we  find  : 


A  =  A-  A  ^  o,?  A. 
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These  are  two  extreme  cases  corresponding  to  a  turbulence  Reynolds  number  either  very  great  (case  of 
exponential)  or  very  small  (case  of  Gaussian)  £9]  .  We  shall  admit  that,  for  intermediary  Reynolds  numbers, 
this  quasi  equality  still  holds.  In  these  conditions,  the  crossed  beam  normalized  correlation,  joined  to  the 
knowledge  of  integration  lengths  2  and  2^1  ,  permits  the  calculation  of  the  integral  scale  : 
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It  should  be  remarked  that  the  orthogonal  beam  method  did  not  involve  integration  lengths.  The  two 
methods  are  thus  independent.  In  both  cases  the  correlations  are  all  the  weaker,  the  more  difficult  to 
measure,  as  these  lengths  are  greater  compared  with  regard  to  the  integral  scale. 

3.1.3  -  Influence  of  beam  width 

Up  to  now  we  assumed  the  beams  to  be  infinitely  thin,  which  is  not  physically  realistic.  As  regards 
uncoupling  of  temperature  and  emissivity  fluctuations,  a  finite  width  does  not  restrict  its  equality  as, 
within  the  framework  of  geometric  approximation,  an  optical  beam  is  a  superposition  of  rays.  On  the  other 
hand,  variance  and  correlation  calculations  can  be  affected,  as  we  then  add  a  complementary  integration  over 
the  beam  crosswise  dimensions. 

The  importance  of  this  effect  has  been  evaluated  in  the  case  of  parallelepiped  beams,  of  cross  section 
1 1  4 ^  4^  for  that  along  axis  1  and  e $  for  that  along  axis  2. 


We  obtain  : 
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where  J\_  is  the  third  integral  scale,  defined  by  : 
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correlation  : 

3  [XX  ~  J?)  ' 

A. 

and  -A.  ,  the  corrective  terms  are 

the  integral  scale,  as  could  be 

For  2-mm-v/ide  beams  and  an  integral  scale  of  about  1  cm,  the  measuring  precision  is  still  1%. 

3.14  -  Application  of  the  method  to  the  study  of  thermal  turbulence 

The  temperature  fluctuation  measuring  technique  and  the  orthogonal  beam  and  crossed  beam  methods  have 
been  applied  on  the  jet  at  the  exhaust  of  a  combustor  for  dilution  studies,  in  spite  of  the  disadvantage 
of  an  heterogeneity  of  mean  temperature  in  this  set-up.  We  processed  in  two  stages  :  first,  measurement 
of  temperature  fluctuations  integrated  on  simultaneous  emission  and  absorption,  then  determination  of  the 
integral  scale.  The  reason  is  that  if  we  assume  that  the  turbulence  characteristics  are  independent  of  the 
measured  parameter  (as  long  as  it  is  a  scalar) ,  the  technique  of  orthogonal  or  crossed  beams  can  also  be 
used  in  emission,  which  presents  the  advantage  of  simplifying  the  instrumentation  and  its  implementation 
by  doing  away  with  the  emitters  of  the  two  necessary  devices.  This  also  yields  a  considerable  increase  of 
the  signal/noise  ratio,  as  it  is  no  more  necessary  to  be  monochromatic. 

a)  Measurement_of_integrated_temgerature_f luctuations 

The  temperature  fluctuation  measurement  has  been  made  after  that  of  the  mean  temperature  by  adjusting 
the  source  temperature  ;  it  is  then  necessary  to  replace  beforehand  the  4.5  pm  interference  filter  by  a 
2.5  ym  filter,  for  which  the  emission  factor  is  lower  than  5%.  The  principle  of  the  measurement  does  not 
require  then  any  more  than  a  single  measuring  channel,  and  the  channel,  interrupted  by  the  knife  edge,  can 
be  used  for  the  determination  of  the  mean  emission  factor  ;  the  temperature  peaks  being  able  to  go  up  to 
100  K  above  the  mean  temperature.  The  r.m.s.  value  of  the  fluctuations  integrated  over  the  optical  path  is 
equal  to  : 

/JF  ^5  K 

b)  0£thogonal_beam_ technique 

So  we  use  two  pyrometers  that  integrate  the  jet  infrared  emission  over  the  whole  spectral  range  up  to 
4.25  ym.  One  of  them,  maintained  fixed,  is  placed  vertically  on  top  of  the  test  section  ;  the  other  is 
horizontal,  and  can  be  translated  axially  thanks  to  a  remotely  controlled  trolley. 

Figure  10  represents  the  spatio-temporal  correlation  curves  obtained  from  these  two  measurements.  Each 
curve  corresponds  to  a  relative  position  of  the  two  beams,  A  being  the  parameter  of  temporal  correlation, 
and  cL  the  shortest  distance  between  the  two  beams.  This  distance  cL  is  modified  by  10-mm-steps,  apart  from 
the  neighbourhood  of  intersection  where  the  pitch  is  reduced  to  2  mm.  In  order  not  to  clutter  the  figure, 
the  corresponding  curves  were  not  traced.  Apart  from  a  slight  attenuation,  the  curves  can  be  deduced  from 
one  another  by  translation.  This  can  be  explained  if  we  admit  the  Taylor  assumption  ;  actually,  if  we  have 
on  a  diagram  the  abscissa  A  ^ax  relative  to  the  maximums  of  the  correlation  curves  as  a  function  of 
distance,  we  obtain  a  straight  line  and  the  velocity  defined  by  the  slope  of  the  line  (84  m/s)  is  close  to 
the  expected  jet  velocity. 


In  case 
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we  deduce  the  crossed  beam  normalized 
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If  we  assume  that  -A-  is  of  the  same  order  of  magnitude  as  A. 
of  order /A ^  .  This  means  that  the  beam  width  must  be  small  besides 


The  points  with  zero  temporal  outphasing  (t  =  0)  permit  the  tracing  of  the  curve  P(ct),  which  presents 
a  maximum  for  d  *  0  (fig.  11).  According  to  Fischer  and  Krause  8  ,  if  we  normalize  this  function  by  the 

value  P(o)  corresponding  to  the  maximum,  we  obtain  the  correlation  function  Rcol}  between  points  A  and  B. 
Its  integration  allows  one  to  calculate  an  integral  scale  ;  we  find  2.5  cm.  In  fact,  the  Fischer-Krause 
interpretation  is  not  correct,  as  it  supposes  negligible  the  correlations  between  points  close  to  A  and  B, 
and  we  saw  that  it  is  necessary,  to  obtain  R(<A,  to  trace  the  correlation  as  a  function  of  parameter  dJ* 
instead  of  (k .,  and  to  take  its  derivative.  Around  the  origin,  this  derivation  is  delicate  ;  to  perform  it, 
we  represent  the  function  by  a  Lagrange  expansion  of  the  7th  degree.  Far  from  the  origin,  the  derivative  is 
determined  graphically.  Thus,  we  obtain  the  curve  of  figure  12  which  represents,  but  for  a  constant,  the 
required  correlation  function.  We  see  that^it  is  much  narrower  than  f(d).  The  value  at  the  origin  makes  it 
possible  to  calculate  the  integral  scaj^e  -X.  ,  and  integration  over  <A.  yields  the  integral  scale  -A.  ;  we 

find  them  very  close  to  each  other  :  -A-  =  0.74  cm  and  _A-  *  0.72  cm,  which  is  a  new  justification  of  the 
approximation  that  consists  in  assuming  them  equal. 

It  should  however  be  noted  that  these  measurements  are  only  approximate.  In  particular,  the  calculation 
by  integration  is  arbitrarily  interrupted  at  the  first  zero  of  the  correlation  curve  ;  we  suppose  that  the 
contribution  of  the  rest  of  the  curve,  alternately  positive  and  negative,  is  negligible. 

c)  Crossed  beam  correlation 


Another  method  to  find  -A-  is  to  use  the  crossed  beam  normalized  correlation,  given  by  : 


P(o)  = 


IT 
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This  determination  involves  the  integration  lengths  2*4  and  2 .  In  the  present  case  5  cm 

and  =30  cm,  we  deduce  an  integral  length  :  _A_  =  0.94  cm,  grea:er  than  the  previous  determination 


but  we  should  note  that,  because  of  the  existence  of  a  profile  of  meai  temperature,  the  jet  height 
does  not  coincide  with  the  width  of  the  emission  zone,  and  leads  to  a:  overestimation  of  integral  scale.  To 
equalize  the  two  determinations  we  should  have  :  =  3  cm.  The  order  of  magnitude  of  these  measurements  is 

three  to  four  times  smaller  than  that  according  to  the  Fischer-Krause  interpretation. 


d)  Fluctuat ion_of _local_temgerature 


Knowing  the  integral  scale,  it  is  now  possible  to  determine  the  r.m.s.  value  of  the  local  temperature 
fluctuations,  using  the  relation  : 


To  perform  this  calculation,  we  retain  for  the  integral  scale  the  mean  value  between  0.94  cm  and  0.72  cm, 
i.e.  0.83  cm.  The  multiplying  factor  is  then  equal  to  1.74  and  the  fluctuation  amplitude  passes  from  25  K  to 
43  K.  The  thermal  turbulence  intensity  is  about  7%. 


e)  Spectrum 

In  principle,  we  can  obtain  the  spectrum  from  the  Fourier  transform  of  the  correlation  function.  But 
this  correlation  function  has  been  obtained  point  by  point  and  its  definition  is  insuficient,  especially 
for  the  points  near  the  origin  which  provide  the  high  frequency  range  of  the  spectrum.  That  is  why  we  were 
satisfied  with  the  direct  spectral  analysis  of  the  measurements,  the  spectra  obtained  being  then  modified 
by  integration  effect.  The  two  spectra  corresponding  for  example  to  the  crossed  beam  measurement,  traced  in 
semi-logarithmic  coordinates  up  to  10  KHz  have  a  non-classical  shape,  almost  exponential,  with  a 
characteristic  frequency  of  about  1300  Hz. 

f )  Turbulent_gremixed_f lame 

This  method  has  been  applied  to  determine  the  distribution_density_f unction  in  the  case  of  a  premixed 
turbulent  flame.  Such  a  flame  is  represented  on  figure  13  :  a  flow  of  hot  burnt  gases  is  parallel  to  a  flow 
of  fresh  gases,  and  a  two-dimensional  flame  develops  :  abscissa  x  is  in  the  flow  direction,  and  the 
temperature  fluctuation  analysis  is  carried  out  at  a  given  abscissa  x  and  at  several  oridinates  y  across  the 
flame  front.  We  observe,  as  shown  on  figure  13,  two  distributions  for  two  temperatures,  apparently  due  to 
cells  of  fresh  and  burnt  gases  |  2 1 3  • 

We  insisted  on  temperature  measuring  techniques  by  optical  methods  as  they  present  a  number  of 
advantages,  but  this  measurement  is  not  devoid  of  some  difficulties  as  regards  its  interpretation. 


3.2  -  Methods  using  probes 


There  exists  a  large  number  of  methods  for  measuring  temperature  fluctuations  by  thermocouples, 
resistors  or  thermistors  ;  as  before  we  shall  describe  one,  showing  the  advantages  and  drawbacks  of  this 
technique.  It  is  that  developed  by  Yule,  Taylor  and  Chigier  pi  ,  based  on  the  use  of  thermocouples,  for 
which  the  most  important  problem,  in  the  case  of  turbulent  flow,  is  to  calculate  the  thermocouple  response 
time,  as  the  gas  temperature  has  to  be  determined  from  the  temperature  measured  by  the  thermo¬ 
couple,  and  the  time  constant  'C  .  ’ 
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The  schematic  of  such  a  thermocouple  is  represented  on  figure  14,  which  shows  how  we  can  design  the 
support  to  minimize  the  support  perturbation.  The  junction  temperature  is  related  to  the  gas 

temperature  l  by  considering  three  mechanisms  :  the  gas-junction  transfer,  the  heat  transfer  through 
the  wire  resulting  from  a  loss  by  conduction,  and  the  convective  and  radiative  heat  transfer  from  the  wire. 
A  study  of  the  energy  budget  is  a  good  approximation  to  give  an  idea  of  the  values  of  the  time  constant  and 
of  the  important  parameters. 

The  arrangement  of  this  thermocouple  induces  one  to  account  for  a  difference  of  time  constant  for  the 
wire  alone  and  for  the  bead  surrounding  the  junction.  In  particular  the  response  time  depends  much  on  the 
size  of  the  bead  relative  to  the  termocouple  wire.  The  thermocouple  response  is  usually  calculated  or 
determined  experimentally,  either  by  realizing  in  the  gas  a  temperature  step  or  by  modulating  at  a  constant 
frequency  the  gas  temperature.  This  operating  mode  is  general  for  the  determination  of  a  time  constant  : 
sine  modulation,  step  or  pulse  technique. 


We  can  create  a  temperature  step  or  pulse  on  the  thermocouple  by  producing  heat  in  the  wire  by  a  step 
or  pulse  of  electric  current. 


As  already  mentioned,  this  time  constant  corresponds  in  fact  to  the  time  constant  ^ 
that  of  the  bead  ,  and  in  its  evaluation  everything  depends  on  the  wire  size  relative  to 
In  the  case  of  a  supporting  wire  of  infinite  length,  the  time  constant  that  is  involved  is 

T  -  T.  v  T,  vje 

3  ip  f  » 

In  the  case  of  a  thermocouple  of  finite  length,  we  use  the  relation  : 


of  the  wire  plus 
that  of  the  bead. 
T .  ,  so  that  : 


T  :  T  tT 
%  f  f  -Jt 
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where  A  is  an  exchange  coefficient  between  the  wire  temperature  at  a  great  distance  and  the  bead 

temperature,  and  the  bead  time  constant.  The  most  exact  method  consists  in  measuring  by  a  method 

of  electric  pulse  heating  and,  by  means  of  a  computer,  evaluating  T*,  and  A,  and  thus  obtaining  the 
corrective  elements  necessary  for  calculating  the  instantaneous  temperature  of  the  gas  from  the  temperature 
"7>>  measured  at  the  junction. 


In  the  design  of  this  device  there  are  two  types  of  computations  : 
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a)  one,  admitting  as  data  the  geometric  parameters  of  the  thermocouple  (wire  and  bead)  and  the  para¬ 
meters  measured  during  use  ;  it  then  can  determine  the  bead  time  constant  in  the  best  conditions,  practically 
those  of  the  experiment,  and  the  corrective  terms  due  to  thermal  losses  ; 

b)  the  other  which,  from  the  time  constant  and  the  measured  temperature,  yields  the  gas  instantaneous 
temperature. 

The  sources  of  error  with  this  technique  arise  from  several  causes  : 

a)  first,  the  inaccuracy  in  the  time  constant  measurement  ; 

b)  the  variability  of  the  knowledge  of  the  geometry  of  the  thermocouples,  that  vary  from  one  to  another 
and  even  during  the  test  ; 

c)  the  uncertainty  of  the  perturbations  brought  by  the  thermocouple  6r  the  probe  to  the  flow,  although 
this  drawback  has  now  been  reduced  through  the  use,  in  flames,  of  wires  with  a  diameter  of  less  than  5  m  ; 
particle  impacts  though  have  to  be  accounted  for  ; 

d)  the  lack  of  data  on  heat  transfer  in  transient  regime  ;  this  important  aspect  is  practically  unknown, 
whence  the  limited  overall  accuracy. 

The  fact  of  using  two  neighbouring  thermocouples,  one  electrically  heated  by  a  high  frequency  current 
and  the  other  not,  eliminates  some  causes  of  errors,  as  the  differences  of  e.m.f.  of  the  two  thermocouples 
could  be  measured  and  be  used  later  to  determine  the  compensation  to  be  applied  to  the  unheated  thermocouple. 

By  this  probing  method  we  determine,  in  principle,  the  temperature  fluctuations  of  one  point,  and  we  can 
thus  obtain  a  Eulerian  statistic. 

3.3  -  Discussions  of  the  temperature  fluctuation  measuring  methods 

We  insisted  on  the  measurement  of  this  scalar  parameter  to  emphasize  the  limits  of  these  techniques  and 
the  progress  that  remains  to  be  made. 

The  temperature  measuring  methods  by  optical  means  are  very  interesting  and  relatively  easy  to 
implement  :  we  indeed  obtain  a  signal  but  we  have  to  interpret  it,  which  is  not  always  easy  ;  we  described 
the  line  reversal  method  in  infrared,  as  it  seems  to  be  relatively  easy,  but  there  exist  other  methods, 
based  for  instance  : 

i  -  on  the  absolute  measurement  of  line  intensities, 
ii  -  on  the  ratio  between  two  line  intensities  (emission  of  lines  due  to  metals  contained  in  the  flame), 

iii  -  on  the  ratio  of  band  intensities  which  one  related,  in  the  case  of  a  diatomic  molecule,  to  the 
internal  motion  of  the  molecule  ;  in  Raman  spectroscopy,  for  instance,  the  measurement  of  the  Raman  line 
intensity  of  the  ratio  cf  intensity  of  the  Stokes  and  anti-Stokes  lines  will  be  examined  in  the  next  section, 
when  concentration  measurement  will  be  studied. 

The  short  response  time  methods  for  measuring  temperature  with  a  probe  present,  as  already  mentioned, 
many  difficulties.  We  chose  to  describe  the  thermocouple  method  because  it  is  a  point  measurement.  The 
resistance  variation  methods  are  usually  coupled  with  hot  wire  techniques,  and  are  used  to  calculate  the 
compensation  term  necessary  for  a  correct  determination  of  velocity  ;  they  present  the  disadvantage  of  not 

being  point  measurements  and  of  integrating  the  unsteady  transfer  phenomena  over  the  probe  length.  These  _ 

coupled  probes  are  very  useful,  as  they  make  it  possible  to  determine  correlations  of  the  type  T'^  ,  -u.'  T\ 
*^T/  »  ,  /V/*  t  aa.***'  9  which,  in  a  two-dimensional  problem,  determine,  apart  from  the  correlations  of 

the  scalar  fluctuations,  the  diffusion  of  this  scalar  in  two  orthogonal  directions  corresponding  to  velocity 
fluctuations  of  the  components  u'  and  vf.  Figure  15  is  an  example  of  realization  due  to  H.  Fiedler  and 
P.  Mensing  [103,  which  relates  the  velocity  effect  to  that  of  temperature. 

4  -  METHODS  OF  INSTANTANEOUS  CONCENTRATION  AND  DENSITY  MEASUREMENT 

The  diagnostic  methods  concerning  combustion  measurements  were  considerably  developed  these  last  few 
years  :  though  they  are  promising,  we  should  discriminate  among  these  methods,  and  this  is  not  always  obvious 
in  the  reports,  those  that  provide  mean  concentrations  and  detect  minor  constituents  such  as  radicals  (OH), 
those  leading  to  instantaneous  concentration  measurements,  hence  of  direct  interest  for  turbulent  flow 
studies,  ar  1  those,  still  instantaneous,  which  do  not  discriminate  the  nature  of  the  molecules  but  rather  are 
sensitive  cd  the  number  of  molecules  per  unit  volume,  i.e.  to  density  fluctuations.  They  all  are  optical 
methods,  and  made  use  of  laser  techniques. 

We  shall  distinguish  in  this  report  three  methods  that  begin  to  yield  some  interesting  results  : 

-  the  spontaneous  Raman  scattering  method, 

-  the  method  making  use  of  fluorescence, 

-  the  coherent  Anti-Stokes  Raman  Scattering  (CARS)  method,  which  we  shall  briefly  describe  and  of 
which  we  shall  give  the  present  performance  ;  it  specially  determines  density. 

4. I  -  Spontaneous  Raman  Scattering  methods  (SRS)  Jill 

When  a  light  beam  crosses  a  gas,  a  small  fraction  of  the  incident  flux  is  scattered,  at  a  different 
freqyency,  over  4tt  steradians.  It  is  the  spontaneous  Raman  Scattering,  in  which  the  incident  photons 
exchange  energy  by  discrete  values  with  the  molecular  system.  The  scattered  flux  spectrum  comprises  lines 
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whose  positions  are  representative  of  the  various  energy  levels  of  the  molecular  system.  Usually  we  are 
interested  with  diatomic  molecules  and  with  vibration-rotation  transitions.  This  type  of  scattering  is  well 
adapted  to  the  analysis  of  gas  mixtures,  as  the  SRS  vibration  spectra  interfere  only  rarely.  Let  us  examine 
the  case  of  vibration-rotation  states  of  the  molecule,  characterized  by  the  quantum  numbers  V  and  3  ;  an  SRS 
process  entails  the  passage  from  a  state  (v,3)  to  the  state  3*^  ),  the  frequency  difference  between 

the  exciting  flux  and  the  scattered  flux  being  representative  of  this  transition.  The  Stokes  scattering 
corresponds,  for  the  molecule,  to  an  energy  gain,  and  the  anti-Stokes  one  to  a  loss.  Let  N  be  the  number  of 
molecules  per  unit  volume,  Tv  and  Tr  the  vibration  and  rotation  temperatures,  the  incident  power, 
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To  implement  SRS,  we  use  an  optical  device,  schematically  represented  on  figure  16,  comprising  an 
exciting  source  and,  at  a  right  angle,  spectrometer  analyzing  the  scattered  light. 


From  the  above  relations  giving  and  we  remark  that  the  SRS  spectrum  varies  as  a  function  of 

W  ,  Tfl  and  Tv  ;  but  we  can  assume  that  the  thermodynamic  equilibrium  is  reached,  and  write  T. 

By  taking  the  ratio  between  the  intensities  of  the  two  lines,  we  obtain  the  temperature  ;  we  may,  for 
instance,  taking  V  =  ! ,  calculate  the  ratio  between  Stokes  flux  emitted  by  the  fundamental  02L  branch 
(due  to  the  molecules  initially  in  the  fundamental  vibration  state)  and  the  Stokes  flux  emitted  by  the 
branch  which  are  initially  in  a  given  excited  vibration  stage,  which  corresponds  to  the  relation  : 

_  |  1  ^ 

^  %  KO) - 3)J  *  ±_  c  ^ 

— ,cn>  1  i 

We  may  also  calculate  the  ratio  between  the  flux  integrated  over  the  whole  Stokes  Q.  branch  and  that 
integrated  over  the  integrated  anti-Stokes  branch  ;  this  ratio  corresponds  to  the  relation  : 

z,,3 

Spectral  selection  can,  in  both  cases,  be  realized  by  interference  filters. 

Temperature  being  known,  the  intensity,  integrated  over  a  spectral  band  of  arbitrary  width  and  position, 
is  proportional  to  the  concentration  in  the  gas  of  the  molecule  chosen.  Whether  the  temperature  is  under  or 
above  1500  K,  we  shall  take  either  the  rotational  or  vibrational  structure.  Thus,  the  intensity  I  integrated 
over  a  band  of  spectral  width  A  X  is  : 

1  =  a  ) 

If  we  take  a  length  ^  =  1  ram,  =  0.66  steradian,  for  nitrogen  with  N  =  2.199  molecules/cm3,  and 

(  o' /  o Z,  )  a  5.4. io3  (excitation  line  4880  A),  we  collect  a  proportion  of  7.10-13,  i.e.  2 . 1 06  photons  per 
incident  joule.  For  a  10%  collecting  optics,  a  10%  detector,  we  detect  2.104  photoelectrons  with  a  relative 
error  of  1%. 

The  expected  performance  is  difficult  to  reach  because  of  spurious  fluxes  :  surrounding  luminosity, 
fluorescence,  particle  radiations,  etc. 

The  exciting  source  can  be  continuous  or  pulsed  ;  unfortunately  it  is  not  possible  to  obtain  the 
fluctuations  with  a  continuous  exciting  source  (the  measuring  operation  must  be  long  if  we  are  to  collect 
enough  photons)  :  but  it  is  possible  to  study  temperature  and  concentration  fluctuations  with  an  excitation 
with  pulsed  laser. 

The  corresponding  diagram  is  presented  on  figure  17.  These  conditions  must  be  fulfilled  : 

-  the  laser  pulses  should  be  short  besides  the  turbulence  characteristic  times, 

-  the  energy  must  be  sufficient  to  allow  a  precise  measurement, 

-  the  measurements  necessary  for  the  determination  of  N  and  T  must  be  performed  simultaneously. 

At  the  present  time,  with  a  30-joule  free-running  ruby  laser  and  for  a  700-vis  laser  pulse,  the  measuring 
time  is  of  the  order  of  10  us,  the  measuring  accuracy  on  N2  concentration  by  measurement  of  the  Stokes  line 
is  of  the  order  of  2%  in  normal  conditions,  and  of  4%  at  1500  K.  The  measuring  volume  corresponds  to  a 
cylinder  1  mm  long  and  1  mm  in  diameter. 
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The  possible  extension  if  this  method  consists  in  realizing  longer  pulses,  of  10  to  20  ms,  with 
measuring  times  of  10  ys,  but  this  extension  is  detrimental  to  accuracy. 

4.2  -  Methods  based  on  fluorescence  (J2j  [.IXl 

This  method  is  similar  to  SRS,  but  instead  of  considering  fluorescence  as  a  parasitic  effect  perturbing 
the  SRS  measurement,  we  consider  it,  on  the  contrary,  as  the  principal  phenomenon  and  use  it  as  a  tracer  for 
the  other  molecules  of  the  combustion  products.  Taran,  Pealat  and  Bailly  [l  lj  showed  that  this  technique  is 
valid  for  determining  the  temperature  and  concentration  of  some  molecules  generating  the  fluorescence 
phenomenon.  The  experimental  technique  is  similar  to  that  of  SRS.  The  spectrograph  can  be  replaced  by  a 
filter  eliminating  the  Rayleigh-Mie  component.  The  photo  multiplier  then  collects  all  the  fluorescence,  as 
well  as  the  individual  Raman  lines,  but  the  contribution  of  the  latter  on  the  total  signal  is  negligible. 
However  it  is  mandatory  that  the  gas  contains  fluorescent  molecules  ;  it  is  the  case  in  pyrolysis  zones  of 
hydrocarbons,  in  particular  in  air-kerosene  flames  and  in  the  600-1000  K  temperature  range  ;  above  1000  K 
this  fluorescence  phenomenon  disappears,  as  the  fluorescent  molecules  are  destroyed. 

The  advantage  of  fluorescence  is  due  to  its  intense  signal,  which  makes  measurements  easier,  at  least 
as  long  as  we  can  interpret  them,  as  concentration  fluctuations  can  be  due  to  several  causes  intervening  in 
the  species  budget  :  first,  a  physical  effect  of  mass  transfer  by  molecular  and  turbulent  diffusitn,  and  also 
a  chemical  effect  of  production  and  destruction,  difficult  to  evaluate  in  a  turbulent  environment.  All 
things  considered,  the  performance  of  this  method  is,  with  a  10-W  continuous  ionized  argon  laser,  a  50  ys 
response  time,  meaning  that  we  can  measure  continuously  frequencies  between  100  Hz  and  20  kHz,  the  spatial 
revolving  power  corresponding  to  a  1 -mm-dia-cylinder  1  cm  long. 

This  method  works  only  between  600  and  1000  K,  and  the  gas  must  contain  hydrocarbons  that  pyrolyze  in 
this  temperature  range.  There  exist  in  the  literature  many  publications  on  the  use  of  the  fluorescence 
phenomenon,  in  particular  on  the  concentration  of  radicals  such  as  CN,  CH,  OH  using  lasers  to  create  a 
saturated  fluorescence,  the  study  being  performed  by  emission  or  absorption  analysis  ;  but  these  methods  are 
not  instantaneous  :  Bonczyk  and  Shirley  £l6j  claim  that  their  method  is  instantaneous,  but  published  no 
result  on  this  point.  The  literature  is  abundant,  but  does  not  apply  to  the  problems  raised  by  the  study  of 
turbulence  * 


4.3  -  Coherent  Anti-Stokes  Scattering  methos  (CARS) 

It  seems  that  the  use  of  SRS  is  delicate,  not  to  say  impossible,  in  the  case  of  unstable  flames  loaded 
with  burnt  or  unburnt  particles,  which  have  in  the  field  too  great  a  luminosity  ;  that  is  why  CARS  is 
important  :  in  particular,  rejection  of  parasitic  light  is  much  easier  than  with  SRS. 


The  principle  is  as  follows  :  we  observe  a  CARS  signal  when  two  plane  and  colinear  waves,  of  frequencies 
C O and  ^5  ,  are  directed  towards  molecules  having  active  vibrations  uu  ^  such  that  ~  -  <^5. 

A  wave  at  the  anti-Stokes  frequency  ^  s  is  scattered  in  the  direction  of  the  exciting  waves, 

and  its  intensity  is  proportional  to  the  square  of  the  concentration  of  the  scattering  molecules.  When  we 
scan  the  frequency  co^  relative  to  ,  the  anti-Stokes  signal  emitted  presents  a  resonant  increase 

whenever  r  oj  c,  is  equal  to  the  frequency  of  a  Raman  transition  of  the  medium.  The  collected 

anti-Stokes  light  constitutes  the  CARS  spectrum  of  the  medium.  The  scattered  beam  is  itself  coherent,  and 
presents  all  the  characteristics  of  a  laser  beam  :  great  intensity,  small  solid  angle,  and  this  in  the  same 
direction  as  the  exciting  waves  :  this  makes  the  elimination  of  parasitic  light  much  easier  than  with  SRS. 
Fluorescence  is  eliminated  all  the  more  easily  as  the  signal  happens  to  be  emitted  on  the  anti-Stokes  side  of 
the  exciting  frequency  ,  while  fluorescence  is  emitted  on  the  Stokes  side. 

The  collected  anti-Stokes  power  is  proportional  to  the  square  of  a  susceptibility  3^  of  the 

medium,  which  expresses  the  response  of  the  molecule  to  the  coherent  excitation  ;  JC  is  the  sum  of  two 
contributions,  and  yL  ^  :  X  p,  is  the  resonant  part,  and  represents  the  contribution  of  the  active 

Raman  vibrations,  of  frequency  v  close  to  uu>  ^  -  <-*->  ^  ;  for  the  Raman  transition  *3  / )  • 

Xp,  <  [  W(O)  -  wtv',3')] 


thus  proportional  to  the  number  of  molecules  per  unit  volume  of  the  molecular  species  possessing  the  Raman 
transition  analyzed  ;  represents  the  very  small  contribution  of  the  other  vibration  motions,  as  well 

as  the  distorsion  non-linearity  of  the  molecule  electronic  cloud  ;  it  is  a  parameter  proportional  to  the 
number  of  molecules  per  unit  volume  of  all  the  molecular  species  present  around  the  measuring  point.  By 

U-)  ^  near  the  vibration  frequencies,  we  obtain  a  spectrum  representative  of  the  variation  of 


varying  uJ  ^ 

|  X  \  with  oU, 


t 


An  example  of  experimental  mounting  is  given  on  figure  18  ;  we  can  see  the  two  lasers,  one  with  a  fixed 
frequency  ,  the  other  with  a  variable  frequency  laj  $  .  But  there  are  experimental  difficulties 

inherent  to  the  nonlinear  nature  of  CARS  : 


-  Interaction  is  efficient  at  the  measuring  point  only  if  there  is  a  spatial  and  temporal  superposition 
of  the  two  excinting  beams,  which  requires  a  careful  alignment  of  the  two  beams. 


-  High  resolution  Raman  spectroscopy  requires  a  great  frequency  stability  of  both  lasers. 

-  High  rejection  power  monochromators  are  necessary  to  isolate  the  anti-Stokes  radiation  to  be 
discriminated  from  exciting  radiations. 


-  Anti-Stokes  power  varying  rapidly  with  the  two  lasers  powers,  the  measurement  has  to  be  made  free  of 
their  power  fluctuations. 


Measurements  of  instantaneous  parameters  are  possible  with  the  CARS  technique,  but  for  the  moment  we 
obtain  density  fluctuation  measurements  with  only  the  following  performance  :  measurement  duration  10  ns 
at  a  10-Hz  rate,  with  a  3%  accuracy  of  the  density  measurement  :  we  thus  obtain  for  the  mean  value  and 
the  standard  deviation. 
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All  these  techniques  will  progress  according  to  progress  on  lasers  (J6j,  [l 7 j  . 

4.4  -  Flow  visualization  and  diagnostic  by  scattering,  emission  and  absorption  techniques 

Global  study  of  the  scattered  radiation  p8j 

Global  methods  have  been  used  to  determine  density  variations  from  scattered  light.  Such  a  set-up 
is  represented  on  figure  19  :  a  laser  beam  is  placed,  in  this  case,  perpendicular  to  a  turbulent  jet  mixing 
with  the  surrounding  air.  The  scattered  light  is  analyzed  by  a  photomultiplier  oriented  at  an  angle  d  with 
the  laser  beam.  From  the  light  scattered  by  the  turbulent  mixture  we  deduce  the  angular  law  of  intensity 
(scattering  indicator)  and  we  analyze  the  scattered  light  by  is  wavelength  spectrum  and  its  fluctuation 
spectra.  Density  fluctuations  results  from  three  fluctuations,  of  pressure,  temperature  and  concentration  : 

This  method  seems  interesting  but  its  interpretation  remains  difficult,  as  witnessed  by  the  work  of 
Lelievre  and  Picard  £l8J. 

To  recover  the  three-dimensional  effect  of  the  phenomenon,  R.  Goulard  (j9j  proposed  a  general  emission 
and  absorption  technique  to  perform  measurements  in  a  three-dimensional  space.  It  is  a  multi-angular  scanning 
technique  such  as  that  represented  on  figure  20,  where  the  profile  as  a  function  of  a  spatial  variable  can 
be  deduced  from  the  variations  of  angles  and  Q  .  This  field  of  research  can  scill  be  improved  by 
making  the  measurement  instantaneous. 

This  last  aspect  leads  towards  the  solution  of  the  problem,  which  consists  in  obtaining  the  time  and 
space  variations  of  a  parameter. 

4.5  -  Density  variation  visualization  by  photographic  methods 

We  noted  at  the  beginnning  of  this  report  the  importance  of  the  knowledge  of  structures  and  the  interest 
of  turbulent  flow  topology.  The  methods  used  at  present  for  better  describing  the  structures  raise  difficult 
problems,  as  the  recorded  information  has  to  be  appropriately  interpreted.  Still  or  motion  picture  techniques 
are  : 


-  use  of  tracers  (solid  or  fluid) , 

-  Schlieren  method, 

-  shadowgraphy , 

-  holography. 

Schlieren  and  shadowgraphy  methods  provide  the  density  fluctuations  of  turbulence,  but  they  integrate 
these  fluctuations  over  the  whole  beam  paths.  Tracing  techniques  and  holography  are  interesting,  as  they 
identify  the  fluid  particle  position,  but  their  interpretation  is  complex  if  we  expect  to  identify  the 
turbulence  structure.  The  holographic  method  collects  the  required  information,  but  the  stochastic  three- 
dimensionality  of  the  process  makes  difficult  the  interpretation  of  this  information.  We  shall  describe  a 
method  that  seems  interesting  for  studying  turbulent  flow  :  it  consists  in  using,  for  the  visualization, 
the  luminescence  of  organic  compounds  such  as  2-3  butanedione  Q20j .  It  also  makes  it  possible  to  introduce  a 
phosphorescent  gas  to  visualize  the  flow.  In  this  way,  one  may  study  turbulence  diffusion  processes  as  well 
as  turbulent  mixing  phenomena. 

This  technique  consists  in  detecting  the  energy  transfer,  by  optical  means,  from  one  molecule  to  another 
by  collisions.  These  collisions  coincide  with  the  mixing  of  species  at  molecular  scale. 

The  molecule  is  excited  by  a  radiation  with  a  known  spectral  distribution  ;  the  incident  light  is 
absorbed  and  re-emitted  in  the  form  of  fluorescence  for  an  excited  singlet  state,  and  in  the  form  of 
phosphorescence  if  the  emission  is  provided  by  an  excited  triplet.  Practically,  this  is  expressed  by  a 
mechanism  represented  on  figure  21  :  the  incident  light,  made  of  several  wavelength,  is  filtered  so  as  to 
emit  only  in  blue  ;  the  radiation  is  absorbed  by  the  excited  butanedione,  and  the  molecule  emits  in  the  green 
range  as  phosphorescence,  and  in  the  blue  range  as  fluorescence,  which  combines  with  the  blue  of  the  incident 
light  ;  we  then  record  on  a  sensitive  plate  only  the  green  radiation  of  phosphorescence  which  is  very  intense. 

The  advantages  of  this  technique  are  the  following  : 

-  it  is  possible  to  discriminate  mixing  due  to  eddies  created  by  turbulence  from  mixing  due  to 
molecular  diffusion  ; 

-  we  can  take  a  direct  measurement  of  the  tracer  concentration  ; 

-  we  can  isolate  part  of  the  flow,  and  analyze  it  ; 

-  we  record  the  mixture  structure,  as  frozen  at  a  given  time. 

There  are  however  drawbacks  in  this  technique  ;  in  particular  : 

-  it  is  only  valid  at  moderate  temperatures,  as  we  have  to  find  compounds  for  which  the  phosphor¬ 
escence  of  fluorescence  emission  remains  visible  in  the  temperature  range  under  study  ; 

-  if  we  wish  to  obtain  instantaneous  photographs  the  flow  velocity  should  not  be  too  high,  so  we 
have  difficulty  in  measuring  velocity  fields  ; 

-  the  information  is  two-dimensional  while  the  phenomena  are  three-dimensional  ; 

-  the  qualitative  analysis  of  the  photographs  is  relatively  easy,  but  quantitative  results  are  very 
difficult  to  obtain. 
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PRESSURE  FLUCTUATION  MEASUREMENTS 


Measuring  in  a  flow  a  mean  stagnation  or  static  pressure  is  within  our  possibility,  but  measuring 
precisely  a  fluctuating  pressure  within  a  flow  remains  a  futuristic  wish. 

Theoretically,  pressure  appears  first  in  the  equation  of  state,  of  the  form,  for  a  medium  of  variable 
composition  :  _  nj  • 

-t = p  r  t  *7- 

being  the  mass  fraction  of  species  j  . 

It  also  appears  in  the  motion  equation  and  in  the  energy  budget  equation.  Closure  problems  involve 
terms  of  correlation  with  pressure,  and  it  is  often  very  difficult  to  appreciate  the  importance  if  these 
terms,  which  are  of  the  form  : 
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which  appear  in  the  turbulence  kinetic  energy  budget.  We  have,  for  instance,  to  compare  the  term  -dp  E. — I 
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where  ^  represents  the  fluctuations  of  dissipation  energy,  J;ft  the  fluctuations  of  mass  diffusion  flux 
in  the  (b  direction,  the  fluctuations  of  mass  chemical  production  rate,  '  the  formation  heat 

of  species  j  ,  °\!fo  t^ie  fluctuations  of  heat  flux  in  the  ft  direction. 

On  the  theoretical  viewpoint  and  as  regards  turbulent  phenomena  modelling,  pressure  fluctuation 
measurements  are  thus  of  prime  importance.  But  at  present  there  only  exist  measurements  by  probe  and  indirect 
measurements. 


-  Measurements  by  probe 

We  possess  pressure  probes  of  small  size  (of  the  order  of  1  mm  in  diameter),  with  a  10  kHz  bandwidth, 
but  they  still  perturb  the  flow  ;  it  is  far  from  ideal  for  taking  correct  pressure  measurements  in  a  turbulent 
flow.  They  are  usually  coupled  with  instantaneous  velocity  measurements  and  correlations  of  the  type 
O'  c  ^  A  / ?  \) )  have  been  obtained.  This  happens  in  isothermal  flows,  as  these  probes  are  still  very  sensitive 
to  temperature  gradients.  Acousticians  use  such  probes  to  analyze  the  acoustic  field  generated  by  flows. 

-  Indirect  methods 

They  are  of  two  kinds  : 

a)  Local  pressure  fluctuations,  determined  from  instantaneous  pressure  measurements  at  the  wall  ;  a  more 
or  less  complex  analysis  makes  it  possible  to  pass  from  the  wall  fluctuation  measurements  to  the  values  of 
instantaneous  local  pressure  ; 

b)  density  and  temperature  measurements  at  the  same  point,  from  which  we  deduce  pressure  fluctuations 
at  that  point. 

As  we  can  see,  pressure  measurements  in  turbulent  flows  are  still  very  inaccurate.  It  is  a  difficult 
problem,  still  far  from  having  received  satisfactory  solutions. 

6  -  CONCLUSION  ON  MEASUREMENTS  IN  TURBULENT  MEDIUM 

This  report  comprises  five  parts  : 

a)  The  first  part  is  devoted  to  the  objectives  of  this  Working  Group,  which  have  widely  changed  in  time. 
We  started  from  a  rather  vague  idea  of  Prof.  A.  Ferri,  which  concerned  pressure  fluctuation  measurements  and 
their  interpretation  in  a  flow  of  hot  gases  and  in  flames. 

This  theme,  concerning  a  very  complicated  problem,  has  not  received  the  approval  of  all  W.G.  Members  ; 
so  we  changes  our  objective  and  defined  a  new  theme,  concerning  scalar  measurements  in  a  hot  turbulent  flow 
and  in  combustion  flow,  as  regards  measurements  of  temperatures,  concentrations  or  densities  and  pressures. 

b)  The  second  part  describes  in  detail  how,  from  these  measurements,  it  will  be  possible  to  better 
understand  the  turbulent  structure  of  mixing  zones  and  flames,  and  in  particular  to  detect  coherence  in 
these  structures.  We  show  the  interest  of  using  photographic  techniques  capable  of  visualizing  these 
structures.  We  have  also  indicated  the  roles  that  turbulent  spectral  measurements  can  play  in  the  description 
of  the  flow.  It  is  not  possible,  for  the  moment,  to  exploit  pressure  fluctuation  measurements.  These  results 
are  important,  as  they  will  be  the  basis  of  any  progress  in  the  modelling  of  turbulent  flows  and  turbulent 
flames,  whether  premixed  or  diffusion,  or  open  or  confined  flames. 
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c)  The  third  part  is  devoted  to  the  instantaneous  measurement  of  temperature.  Two  techniques  are 
described  :  optical  methods  and  probing  methods.  We  insisted  on  the  difficulty  of  interpreting  optical 
methods,  as  they  integrate  phenomena  over  part  of  the  beam  ;  these  interpretations  require  many  assumptions. 

The  probing  methods  require  the  knowledge  of  data  allowing  the  evaluation,  from  probe  measurements,  of  the 
real  fluctuations  of  the  gas  temperature.  These  methods  are  interesting  if  they  are  coupled  with  velocity 
measuring  probes. 

d)  The  fourth  part  deals  with  instantaneous  measurements  of  concentration.  These  techniques,  which 
take  advantage  of  the  laser  performance,  are  at  present  in  full  evolution.  Three  methods  are  particularly 
described  :  Spontaneous  Raman  Scattering  (SRS) ,  methods  making  use  of  the  fluorescence  perturbating  phenomenon, 
and  Coherent  Anti-Stokes  Raman  Scattering  (CARS).  A  large  number  of  studies  showed  that  these  three  methods 
are  interesting,  but  few  of  these  methods  are  Adapted  to  the  instantaneous  concentration  measurement.  Only 

SRS  yields,  in  its  pulse  version,  instantaneous  concentrations  ;  fluorescence  also,  but  within  a  narrow 
temperature  range  ;  CARS,  for  the  moment,  determines  only  density  fluctuations.  In  this  same  section,  we 
discussed  global  scattering  methods,  and  appropriate  techniques  for  treating  the  three-dimensional  problem. 

We  had  initially  planned  a  section  on  the  coupling  of  these  scalars  with  the  measurements  of 
instantaneous  velocity  ;  we  suppressed  it,  as  it  has  been  discussed  in  each  part,  and  we  included 
visualization  techniques,  which  in  fact  rest  on  density  measurements.  We  gave  a  few  indications  on  the 
technique  of  phosphorescent  tracing  gases,  and  on  the  advantages  and  drawbacks  of  this  technique. 

e)  The  fifth  part  is  devoted  to  the  instantaneous  measurement  of  pressure.  This  measurement  is 
difficult  in  situ.  We  must  pass,  to  reach  it,  through  an  a  priori  modelling  of  the  measuring  zone. 


BIBLIOGRAPHY:  See  page  18 
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Fig.  1  -  Spectre  des  fluctuations  de  masse 
volumique  dG  aux  fluctuations  de  concentration 
dans  le  jet  froid  (1  et  2)  et  dans  la  flamme 
de  diffusion  turbulente  (3,4).  Pres  de  la 
tuyere  pour  I  et  3  (0  a  200  fois  la  largeur  de 
la  tuyere)  et  a  200  a  400  fois  la  largeur  pour 
2  et 4  (d'apres  Souiouchev  et  Sokolenko  [l]) 

-  Density  fluctuation  spectra  due  to 
concentration  fluctuations  in  a  cold  jet 
(curves  1  and  2)  and  in  a  turbulent  diffusion 
flame  (curves  3  and  4) .  Measurements  are  close 
to  the  nozzle  (0  to  200  times  the  nozzle 
width)  for  curves  1  and  3  and  200  to  400  times 
this  width  for  curves  2  and  4.  (after 
Souiouchev  and  Sokolenko  1 1]  ) 
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Fig.  4  -  Coupe  d  un  remous  tourbillon 
qui  se  forme  dans  la  zone  de  transi¬ 
tion  d'une  flamme  de  diffusion  [3] 

-  Cross  section  of  an  eddy- 
vortex  forming  in  the  transition 
zone  of  a  diffusion  flame  [3] 
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Fig.  5  -  Coupe  d'un  remous  tourbillon  dans  une 
flamme  de  diffusion  [3} 

-  Cross  section  of  an  eddy-vortex  in  a 
diffusion  flame  £3] 
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Fig.  6  -  Principe  de  la  methode  de  mesure  de  la  temperature  d'un  gaz  par  emission  et  absorption 
simultanees  [6j 

-  Principle  of  the  method  for  measuring  gas  temperatures  by  simultaneous 

emission  and  absorption  [6] 
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13  -  Determination  de  la  densite  de  probabilite  par  une  mesure 


de  temperature 


1 2  lj 


-  Determination  of  probability  density  by  temperature  measurement 


[21] 
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Fig.  14  -  Schema  de  la  sonde  a  thermocouple  a  fil  fin  [9] 
-  Schematics  of  thin  wire  thermocouple  probe  [9] 


Thermo  -  sonde 
Thermo-probe 


Fig.  15  -  Sonde  de  mesure  de  la  vitesse  et  de  la 
temperature  [loj 

-  Velocity  and  temperature  measuring  probe  [lo] 
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Fig.  16  -  Schema  simplifie  d'un  montage  de  DRS  [l5] 
-  Simplified  diagram  of  SRS  [l5l 
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Fig.  17  -  Montage  experimental  de  DRS  en 
impulsions  [15} 

-  Experimental  arrangement  of 
pulsed  SRS  [15] 


Fig.  18  -  Dispositif  experimental  de  DRASC  [l 4}  . 
LP  -  laser  fournissanr  la  frequence  w^. 

LS  -  laser  accordable  fournissant  la  frequence  wg 

M  -  miroir  ,  D  -  lame  dichroique 
SP  -  lame  semi  ref lechissante  prelevant  une  par- 
tie  des  faisceaux  excitateurs  pour  les  envoyer 
dans  une  cuve  de  reference  ,  L  -  lentilles 
F  -  filtres  laissant  passer  le  signal  anti- 
Stokes  et  bloquant  les  frequences  w  et  wg 

PM  -  photomultiplicateurs. 

-  Experimental  arrangement  of  DRASC  [l4] 
LP  -  laser  providing  the  fixed  frequency  w 

LS  -  tunable/laser  provinding  the  frequency  wg 

M  -  mirror  .  D  -  beam  splitter 

SP  -  beam  splitter  picking  up  part  of  exciting 

beams  to  direct  them  to  reference  cell 

L  -  lenses  ,  F  -  filters  passing  the  anti- 

Stokes  signal  and  blocking  frequencies  w  and 

wg  ,  PM  -  photomultipliers. 


Fig.  19  -  Dispositif  experimental 
d'analyse  de  la  lumiere  diffusee 
[18] 

-  Experimental  device  for 
analyzing  scattered  light  £l8] 
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Fig.  20  -  Absorption  par  faisceaux  en  eventail  (Acta-Scanner)  [l9] 
-  Fan-shaped  beams  in  absorption  (Acta-Scanner)  [l9j 
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Fig.  21  -  Excitation  directe  du  butanedione  £2oJ 


-  Direct  excitation  of  budanedione  molecule 
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14.  Abstract 

-  In  1976,  the  Propulsion  and  Energetics  Panel  set  up  its  Working  Group  09  on  ‘Turbulent 
Transport  Phenomena’.  The  objective  was  to  improve  knowledge  of  scalar  fluctuation  in 
turbulence  with  or  without  chemical  reaction  in  variable  density  flows  and  to  develop  techniques 
for  measurement  of  scalar  fluctuations,  especially  static  pressure. 

This  report  contains  the  results  of  the  Working  Group.  After  the  description  and  the  discussion 
of  the  temperature  measurement  methods  (optically  and  using  probes)  the  measurement 
methods  of  non  constant  concentration  and  density  are  treated.  After  regarding  the  possibilities 
of  fluctuating  pressure  measurement  some  conclusions  are  drawn. 

This  Advisory  Report  was  prepared  at  the  request  of  the  Propulsion  and  Energetics  Panel. 

En  1976,  le  Panel  Energetique  et  Propulsion  crea  le  Groupe  d’Etudes  09  sur  ‘Les  Phenomenes  de 
Transport  Turbulent’.  Ce  Groupe  avait  pour  objectif  l’amelioration  des  connaissances  sur  les 
fluctuations  scalaires  de  la  turbulence  avec  ou  sans  reaction  chimique  dans  les  ecoulements  a 
densite  variable  et  la  mise  au  point  de  techniques  de  mesure  de  fluctuations  scalaires,  en 
particulier  de  pression  statique. 

Ce  rapport  contient  les  resultats  des  activites  de  ce  Groupe  d’Etudes.  Apres  une  description  et 
un  examen  des  methodes  de  mesures  de  temperatures  (a  l’aide  de  disposilifs  optiques  et  de 
sondes),  sont  etudiees  les  methodes  de  mesure  de  concentration  et  de  densite  non  constante.  Un 
tour  d’horizon  des  possibility  de  mesures  de  pression  fluctuante  est  suivi  de  la  presentation  de 
certaines  conclusions. 
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